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SUMMARY 


This  study  was  initiated  to  provide  documentation 
and  guidelines  for  the  application  of  finite  element  analysis 
(FEA)  to  perform  mechanical  and  thermal  analyses  of  microelec¬ 
tronic  packages.  The  documentation  includes  an  investigation  of 
the  use  of  mechanical  engineering  analyses  as  a  viable  means  of 
assessing  the  reliability  of  microelectronics  devices.  The 
guidelines  show  how  to  apply  FEA  to  electronic  devices. 

In  the  finite  element  technique,  a  complex  structure 
is  modeled  using  an  assemblage  of  finite  elements  (rods,  beams, 
plates,  or  solids).  During  the  finite  element  analysis  values 
of  a  field  quantity  (for  example,  displacement)  are  estimated 
throughout  the  element  by  the  use  of  i  nterpol ati on  functions. 
The  stiffness  matrix  for  each  element  is  then  obtained  by 
applying  one  of  several  possible  energy  principles.  The  FEA 
computer  program  then  obtains  the  stiffness  matrix  for  the 
entire  structure  by  properly  assembling  the  stiffness  matrices 
of  all  the  elements  using  energy  principles  and/or  constitutive 
laws.  Inputting  nodal  forces  and  nodal  boundary  conditions 
results  in  an  output  of  the  field  variable  at  each  node. 

The  technical  approach  to  this  study  was  to:  1) 

perform  a  technical  assessment  to  identify  microelectronic 
package  problems  and  computer  programs ,  2)  evaluate  the  computer 
programs  and  perform  sample  calculations,  3)  perform  finite 
element  analysis  on  a  critical  problem,  and  4)  instrument  and 
test  an  empty  hybrid  package  for  comparison  to  FEA  and  closed- 
form  analyses  to  assist  in  preparing  guidelines  for  analysis. 

This  study  has  shown  that  although  FEA  is  primarily 
a  tool  for  designers  of  large  complex  structures,  a  great 
potential  exists  if  applied  to  microelectronic  packages.  By 
using  FEA  for  early  warnings  of  potential  failures,  the 
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mechanical  engineer  can  play  a  significant  role  In  assuring  or 
improving  the  reliability  of  microelectronic  systems.  The  most 
frequent  microelectronic  package  problems  were  hermetic ity, 
broken  or  corroded  internal  wires,  and  broken  wirebonds.  The 
survey  showed  that  most  companies  do  not  use  FEA  for  microelec¬ 
tronic  packages  and  that  considerable  cost  savings  can  be 
realized  by  substituting  FEA  for  MIL-STD-883  screening  tests. 

The  computer  program  evaluation  resulted  in  the 
following  recommended  programs  for  microelectronics  package 
analysis:  ABAQUS ,  ANSYS,  NISA,  and  STARDYNE.  STARDYNE  was 
found  to  be  the  best  program  for  analyzing  the  majority  of 
microelectronic  problems. 

Recommended  future  work  related  to  this  report 
includes;  understanding  the  sealing  of  large  custom  packages, 
inspection  criteria  for  hermetic  chip  carrier  solder  joints, 
selective  MIL-STD-883  qualification  by  analysis,  work  in 
nonlinear  analysis,  and  a  pre  and  postprocessor  computer  program 
study  linked  to  CAD/CAM  activities. 

The  reader  is  referred  to  the  text  of  the  report  for 
a  more  detailed  discussion  of  the  study.  Each  section  is 
concluded  with  a  detailed  summary,  while  a  condensed  section  at 
the  end  of  the  report  presents  all  of  the  conclusions  and 
recommendati ons  . 

It  is  refreshing  that  RADC  took  the  initiative  to 
sponsor  a  study  of  this  nature.  It  is  only  by  the  knowledge 
gained  in  this  type  of  study  that  we  can  advance  the  state-of- 
the-art  for  mechnical  design  and  analysis. 


1.0  INTRODUCTION 

The  microelectronic  package  is  one  of  the  key 
elements  in  electronic  ecuipment  reliability  i  h  cost.  However, 
little  information  is  available  for  assessing  the  overall 
acceptabi 1 i ty  of  the  various  mi croel ectronic  packages  in  the 
early  design  stages.  The  intent  of  this  study  is  to  evaluate 
existing  finite  element  computer  programs  to  be  used  in 
assessing  the  reliability  of  mi croel ectronic  devices. 

/ 

2  .0  PROGRAM  DESCRIPTION 

2  . 1  Background 

Finite  element  analysis  (FEA)  can  be  compared  to  a 
puzzle  where  the  structure  or  picture  is  broken  down  into  small 
non-uniform  parts.  With  a  puzzle,  one  must  assemble  the 
individual  pieces  to  reconstruct  the  picture.  In  finite  element 
analysis  the  structure  is  assembled  using  the  rules  of  struc¬ 
tural  mechanics:  eouilibrium  of  forces  and  continuity  of 
di  s  pi acements  . 

FEA  has  been  applied  to  large  complex  structures  as 
a  mechanical  and  structural  analytical  tool.  These  same  FEA 
techni cues  can  be  applied  to  the  evaluation  of  microelectronic 
packages  for  structural  and  thermal  integrity  when  exposed  to 
military  specifications  such  as  MIL-STD-883  and  MIL-STD-810  . 


2 .2  Program  Objectives 

The  objective  of  the  program  was  to  evaluate  new  and 
innovative  analytical  methods  to  ensure  the  mechanical  integrity 
of  microelectronics  packages  used  for  Hybrid  Integrated  Circuits 
{ H I C )  ,  Microwave  Integrated  Circuits  (MIC),  and  Stripline 
Microwave  Circuits  (SMC).  The  mechanical  computer  analysis 
could  be  applicable  to  HIC,  MIC,  and  SMC,  as  well  as  any 
ml c roel ec tron i c  package.  The  data  generated  by  computer 
analysis  will  be  compared  to  the  data  generated  by  physical 
environmental  testing  of  a  particular  hybrid  package.  The 
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programs  used  for  the  analysis  have  adaptability  to  the  various 
thermal  and  mechanical  properties  of  the  package  types. 

8elow  is  a  concise  summary  of  the  study  objectives: 
o  Investigate  the  feasibility  of  exposing  failure 

modes  by  use  of  computer  analyses  of  package 
types  rather  than  physical  testing, 
o  Evaluate  feasibility  of  reducing  cost  by 

performing  FEA  as  an  acceptable  alternative  to 
MIL-STD-883,  Group  0  package  evaluation, 
o  Conduct  a  survey  of  commercially  available 

finite  element  and  thermal  analysis  programs, 
o  Generate  a  set  of  guidelines  documenting  the 

practical  use  of  some  of  these  programs, 
o  Compile  a  list  of  mi c roel ect ron i c  problems  that 

can  be  evaluated  using  mechanical  analyses, 
o  Perform  a  solution  of  a  critical  problem  by  the 
finite  element  method,  compared  to  the  closed- 
form  method,  with  correlation  to  test  results. 
The  next  section  outlines  the  program  in  chronologi¬ 
cal  order  of  the  various  activities. 

2  .3  Technical  Approach 

Below,  in  outline  form,  is  a  brief  description  of 
the  entire  program. 

Phase  1:  Technical  Assessment 

o  Identify  microelectronic  package  problems  that 
could  be  solved  by  analytical  methods, 
o  Determine  the  extent  to  which  mechanical  and 
thermal  analyses  of  microelectronic  packages  are 
performed  by  industry. 

Phase  2:  Analysis  Methods 

o  Compare  and  assess  general  purpose  public  domain 
finite  element  computer  programs . 
o  Compare  finite  element  methods  with  closed-form 
methods . 
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Perform  sample  finite  element  analyses  on 
typi cal  problems . 

Phase  3:  Critical  Problem  Analysis 

o  Select  a  critical  problem  and  appropriate  finite 
element  computer  program  to  show  how  a  finite 
element  analysis  can  solve  a  useful  problem,  how 
it  compares  with  a  closed-form  solution,  how  it 
compares  with  test  results,  and  as  a  demonstra¬ 
tion  of  the  application  of  the  guidelines. 

Phase  4:  Correlation  to  Test  Data 
o  Instrument  and  test  an  empty  hybrid  package 
o  Prepare  guidelines  for  finite  element  analyses 
of  microel ectroni c  packages. 

3.0  TECHNICAL  ASSESSMENT  (PHASE  1) 

3 .1  Introducti on 

A  technical  assessment  was  deemed  necessary  to 
determine  which  mi croel ectroni c  package  problems  could  be  solved 
by  analytical  methods,  and  the  extent  to  which  mechanical  and 
thermal  analyses  of  microel ectronics  are  performed  by  industry. 
Included  in  this  survey  was  an  assessment  of  computer  programs, 
closed  form  techniques,  and  correlation  to  test  data. 

The  technical  assessment  consisted  of  sending  some 
43  questionnaires  to  representati ve  microelectronics  manufac¬ 
turers  and  users.  Several  government  agencies  were  also  polled. 
A  copy  of  the  questionnaire  is  presented  as  Table  1.  Other 
reference  sources  were:  a)  Report  Bibliography  from  the  Defense 
Technical  Information  Center,  b)  Shock  and  Vibration  Information 
Center,  c)  personal  references  and  contacts  at  other  companies. 
The  two  most  useful  information  sources  for  this  study  were  the 
Shock  and  Vibration  Information  Center  and  the  personal  refer¬ 
ences.  The  Questionnaire  survey,  in  theory,  should  provide  the 
most  data;  however,  this  was  not  the  case.  It  is  believed  that 
for  the  survey  to  be  effective,  each  company  must  be  visited  to 
obtain  a  meaningful  response. 
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1.  Do 


'""BLE  1  (cont'd) 

CUESTICNNAIRE  FOR  FINITE  ELEMENT  ANALYSIS  OF  MICROELECTRONIC  PACKAGES  (cont'd) 


4.  Content  on  the  impact  of  MEL-STD-883  on  your  end  product. 


5.  Are  there  any  outstanding  problems  in  the  design,  fabrication,  test 
or  utilization  practices  for  microcircuit  packages? 


6.  General  content  -  please  discuss  any  problem  that  has  not  been  addressed 
above. 


3  .2  Field  $  urvey 

The  search  for  microelectronic  package  data  (and 
finite  element  analysis  of  the  same)  was  fairly  successful.  As 
in  any  search,  it  was  not  all  inclusive;  however,  it  is  felt 
that  the  results  represent  actual  trends  in  package  problems  and 
analysis.  Below  is  a  brief  description  of  each  search  area  and 
the  pertinent  results. 

3.2.1  Questionnaire  Survey 

Forty-three  (43)  Questionnaires  (see  Table  1  for 
Questionnaire  form)  were  sent  to  various  hybrid  manufacturers, 
users,  and  package  suppliers.  Fourteen  (14)  responded  by 

returning  the  completed  forms.  This  is  an  average- to- above- 
average  response  for  this  type  of  survey.  Table  2  presents  a 

summary  of  the  survey  results.  Note  that  several  respondents 

did  not  complete  all  sections  of  the  Questionnaire. 

From  Table  2  it  is  seen  that  the  number-one  package 
related  problem  is  lack  of  hermeticity,  the  number-one  wire 
problem  is  a  broken  or  corroded  wire,  and  the  number-one 
wirebond  problem  is  broken  bonds. 

One  purpose  of  this  study  is  to  determine  which 

mi c roel ectron ic  problems  can  be  analyzed  using  FEA  or  closed- 
form  solutions,  and  now  is  an  opportune  time  to  discuss  this 
topic  as  well  as  to  describe  microelectronic  packages. 

A  definition  of  a  microelectronic  package  is  very 
difficult,  since  it  encompasses  many  types  of  parts  and 
components.  It  would,  however,  be  useful  to  give  a  list  of  some 
of  the  containers  used  to  package  microelectronic  packages, 
they  are; 

Flatpack  Package  (Figure  1) 

Metal  Plug-In  Package  (Figure  2) 

"TO"  Transistor  Outline  (Figure  3) 

D  ual  -  i  n-1  i  ne- pack  age  (OIP)  (Figure  4) 

Microwave  Flat  Package  (Figure  5) 

Multiple  Cavity  Package  (Figure  6) 

Leaded  Chip  Carrier  (Figure  7) 

Leadless  Chip  Carrier  (Figure  8) 
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Table  2.  Summary  of  Survey  Results 


1.  Use  of  mechanical  engineering  analysis 

Computer  Closed-Form 
Analysis  Analysis  Test  Data 

Structural  and  Dynamic  43  5 

Thermal  6  3  7 

2.  Devices  that  the  above  analysis  or  test  data  are 
applied  to  for  solving  mechanical  problems. 

Hybrids,  DIPS,  Flatpacks,  and  any  special  modules 

3.  The  following  matrix  chart  is  to  ascertain  the 
relationship  between  package  types  and  related 
problems  . 


Related  Problems 

Type  of 

Hybr i d 

Mi crowave 

S  t  r  i  p  1  i  n  e 

Microelectronic  Integrated 

Integrated 

Mi c  rowave 

Package  Problem 

Ci  rcui  t 

Ci  rc  ui  t 

C  i  rc  ui  t 

Other 

Lack  of  Hermeticity 

8 

3 

2 

1 

Die  bond  fail ure 

4 

1 

_ 

.. 

Broken  external  lead 

3 

1 

1 

- 

External  lead  corrosion 

5 

3 

3 

2 

External  lead  fatigue 

4 

1 

1 

Solder  Reject 

3 

1 

1 

«■ 

Broken  wire 

3 

2 

2 

1 

Corroded  wire 

3 

2 

2 

1 

Shortened  wire 

2 

2 

2 

1 

Broken  bond 

6 

2 

1 

1 

I ntermetal 1 i c  formation 

5 

2 

2 

1 

Lifted  bond 

3 

2 

1 

1 

Misplaced  bond 

4 

2 

1 

1 

Multiple  bonds 

2 

2 

2 

2 

Overbonded 

3 

1 

- 

- 

4.  Impact  of  MIL-STD-883  on  end  product 
o  Makes  for  more  reliable  product 

o  For  large  packages  (2"  x  2”)  Method  1014  seal  test  at 
30  PSIG  minimum  bomb  pressure  may  cause  breakage  of 
feed-through  glass  beads 

o  Some  use  883  as  a  general  specification  and  use 
internal  spec's  with  more  realistic  testing 
o  MIL-STD-883  Is  more  expensive  than  the  above  approach 
o  No  destructive  analysis  required 
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Table  2.  Summary  of  Survey  Results  (Continued) 


Outstanding  problems  in  the  design,  fabrication,  test,  or 
utilization  practices  for  microcircuit  packages 
o  Need  better  knowledge  on  the  sealing  of  large  custom 

packages 

o  There  exsits  a  basic  lack  of  knowledge  on  fracture 

failures  of  brittle  materials  used  in  packages  such  as 
glasses  and  ceramics 

o  Meniscus  crack  criteria  needs  to  be  established 
General  Comments  - 

o  Long  lead  times  for  finished  hybrid  products, 
o  Large  packages  usually  have  more  difficulty  in  sealing 
than  smal 1 er  ones  . 

o  Need  to  develop  a  criteria  for  inspecting  mounted 
hermetic  chip  carriers  (leaded  or  leadless), 
o  The  industry  may  be  approaching  the  position  to  have  to 
pay  for,  and  receive  rights  to  test  lot  samples  for 
seal  integrity. 


FIGURE  1  FLAXPACK  Panra/tP 


FIGURE  2  METAL  PLUG-IN  PACKAGE 


FIGURE  3  "TO"  TRANSISTOR  OUTLINE 
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Figure  4.  Outl-in-Une  Packages  Mounted 
to  a  Printed  Circuit  Board 


C  BAND  -  UHF  DOWNCONVERTER 
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Pot*  "till  Pro  Worn  Artis 
— Scnw  Sicking  out 
!— Sufcstriti  Moulting  JU  Brtaktg* 
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Potential  Problem  Area 
1— Solder  Joint  Cracking  Due 
to  Extreme  Temperature  Cycling 
or  Vibration  and  Shock 


HYBRID  SUBSTRATE 
WITH  COMPONENTS 
MOUNTED 
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These  packages  are  used  to  package  Hybrid  Integra¬ 
tion  Circuits  (HIC),  Microwave  Integrated  Circuits  (MIC),  and 
Stripline  Microwave  Circuits  (SMC),  or  any  other  mi croelectroni c 
circuits,  such  as  Large  Scale  Integration  (LSI)  or  Very  High 
Speed  Integrated  Cricuits  (VHSIC).  Figures  9,  10,  and  11 

present  other  examples  of  microelectronic  packages. 

Below  is  a  direct  Quote  from  [9]  which  further 
describes  microelectronic  packaging. 

"Microelectronic  packaging  protects  and 
supports  electronic  devices  and  circuits 
and  provides  connections  to  the  other 
parts  of  the  system.  The  protection 
function  avoids  mechanical,  electrical, 
chemical,  contamination,  and  photo- 
optical  damage,  degeneration,  and  causes 
of  malfunction.  Hybrid  microelectronic 
circuits  and  subsystem  packages  support 
the  substrates;  the  substrate  contains 
the  circuit  elements,  (semiconductor 
devices  or  IC  chips,  deposited  or  chip 
resistors  and  capacitors,  and  attached 
inductors),  as  well  as  deposited  and 
bonded  interconnection  wires.  The 
connections  to  other  parts  of  the  system 
include  electrical  leads,  heat  removal 
paths,  and  mounting  functions.  At 
present,  in  order  to  meet  the  demands  of 
VLSI,  the  emphasis  is  on  packages  with 
higher  densities  while  maintaining 
performance,  reliability,  and  low  cost." 

In  discussing  mi croel ectroni c  packages,  it  becomes 
necessary  to  discuss  microelectronic  packaging  on  the  Printed 
Circuit  Board  (PCB)  level.  This  is  required  since  this  is  one 
level  of  potential  failure.  Failures  can  occur  during  part 
screening  or  qualification  testing.  Failures  can  also  occur  at 
the  integrated  assembly  qualification  testing  (i.e.,  micro- 
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electronics  packages  mounted  to  PCB's  which  are  in-turn  secured 
to  an  electronic  housing  assembly).  Failures  may  also  occur 
during  field  service.  So  it  is  seen  that  one  must  be  able  to 
assess  potential  failures  at  any  cycle  in  the  life  of  a  micro¬ 
electronic  package.  FEA  can  be  applied  to  any  or  all  of  these 
potential  failures  where  dynamic,  static,  or  thermal  stresses 
are  imposed.  Be low  is  a  condensed  list  of  microelectronic 
packages  potential  problem  areas,  some  of  these  are  indicated  in 
Figures  1  through  11  . 

o  Lack  of  hermeticlty  -  This  can  occur  when  the 
lead  gl  ass-to-metal  seal  fails,  when  the  lid 
seal  fails,  or  when  leakage  occurs  through  any 
seal  or  connector. 

o  Die  bond  failure  -  When  the  die  (or  chip)  bond 
is  stressed  beyond  its  allowable  strength,  a 
bond  failure  will  occur.  This  can  result  from 
thermal  shock,  thermal  stress,  or  static  and 
dynamic  loading  (see  Figure  9). 

o  Broken  external  lead  -  When  the  external  lead  is 
stressed  beyond  its  endurance  or  ultimate 
strength,  a  failure  will  occur.  This  can  be 
accelerated  if  scratches  or  nicks  are  present 
from  an  improper  lead  bending  process. 

o  External  lead  corrosion  -  If  proper  lead  plating 
is  not  specified,  then  lead  corrosion  can  result 
in  lead  fail ure  . 

o  Wire  bond  and  wire  failures  -  The  wire  bond  is 
the  connection  between  the  chip  or  die  and  the 
external  lead  of  the  package.  An  example  is 
given  in  Figure  7.  These  wires  are  usually 
0.001"  diameter  gold  or  aluminum.  Wire  bond 
failures  occur  from  Improper  bonding  techniques, 
over-stress  conditions,  and  corrosion.  Wire 
failures  generally  result  from  corroded  wire  or 
an  over-stress  condition. 
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Microwave  package  failures  -  These  packages  are 
typically  very  rugged  and,  therefore,  structural 
or  lead  failures  are  rare.  Figures  5  and  6  show 
some  of  the  expected  failures.  They  are: 
substrate  tie-down  screws  backing  out  due  to 
improper  torque  or  over-stress,  substrate 
mounting  tab  breakage  due  to  over-stress  and 

feed-through  connector  breakage  due  to  thermal 

stress  build-up. 

Leadless  chip  carrier  solder  joint  failures  - 

this  is  probably  one  of  the  most  discussed, 

analyzed,  and  tested  areas  in  the  micro¬ 
electronics  area  over  the  past  four  years. 
Figure  8  presents  empty  leadless  chip  carriers 
mounted  to  a  polyimide  substrate.  Solder  joint 
failures  result  from  a  difference  in  coefficient 
of  thermal  expansion  between  the  ceramic 
leadless  chip  carrier  and  the  polyimide  or  G-10 
substrate . 

Air  wound  inductor  breakage  -  Figure  10  presents 
an  RF  hybrid  multiplier  with  air  wound 
inductors.  At  certain  dynamic  levels,  these 
inductors  will  either  physically  break  or  will 
be  out  of  the  electrical  spec i f  i  ca t  i  on  for  which 
they  were  intended. 

Connector  lead  breakage  or  connector  wear  - 
While  a  connector  is  not  a  microelectronic 
package,  it  is  a  part  of  the  overall  packaging 
concept  and  a  connector  failure  will  cause 
system  failure  the  same  as  a  package  failure 
will.  Figure  11  shows  a  connector  mounted  to  a 
PCB  populated  with  hybrids.  Two  failures 
generally  occur  in  the  connector;  they  are:  1) 
lead  breakage,  due  to  excessive  relative  motion 
between  the  connector  and  PCB,  and  2)  connector 
pin  wear  from  extended  or  high  vibration  levels 


causing  i ntermi ttents  in  the  system  operation. 
A  third  connector  failure  should  be  mentioned 
although  it  is  unique  to  a  particular  design. 
There  have  been  reports  of  the  connector  solder 
joints  cracking  or  pins  breaking  due  to  a 
mismatch  in  the  coefficient  of  thermal  expansion 
between  the  connector  body  and  the  PCB. 

3  .2  .2  Shock  and  Vibration  Information  Center 

The  Shock  and  Vibration  Bulletin  Index  was  r  viewed 
for  applicable  material  to  be  used  for  this  study.  The  first 
pass-through  resulted  in  no  obvious  references.  At  this  point 
it  was  decided  to  inspect  each  paper  from  1978  to  1980.  This 
resulted  in  eight  excellent  references  for  FEA  or  closed-form 
analysis  compared  to  test  data.  It  was  concluded  that  a  review 
of  the  reference  titles  would  not  necessarily  indicate  whether  a 
paper  discussed  FEA.  It  should  be  noted,  however,  that  no 
references  were  found  on  mi  croel  ectroni  c  FEA.  This  does  not 
mean  that  these  articles  do  not  exist;  however,  there  are 
probably  not  many  such  articles.  Many  analysts,  the  authors 
included,  perform  such  analyses  on  a  daily  basis,  but  do  not 
present  the  results  at  every  symposium.  Reference  [1]  presents 
the  dynamic  analysis  of  a  torpedo  shell  modeled  with  plate  and 
beam  elements.  NASTRAN  was  used  as  the  Finite  Element  Computer 
Program . 

Overall,  there  was  very  excellent  agreement  between 
the  FEA  solution  and  the  test,  showing  once  again  that  a 
properly  constructed  finite  element  model  will  give  accurate 
results.  This  fact  was  repeated  over  and  over  again  throughout 
this  search . 

Reference  [2]  is  another  example  of  FEA  applied  to 
large  structures  (helicopters)  with  test  data  comparisons.  The 
results  were  excellent  over  most  of  the  frequency  range,  but  at 
higher  frequencies  the  accuracy  deteriorated.  The  authors 
believe  that  this  reduction  in  accuracy  was  due  to  modeling 
assumptions  rather  than  the  models  themselves  or  the  FEA.  This 
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shows  that  one  must  understand  all  the  analytical  assumptions 
and  their  effects  on  the  final  results. 

Reference  [3]  is  a  different  application  of  FEA.  It 
compares  the  FEA  and  an  exact  solution.  The  FEA  uses  the  STAGS 
USA  (Underwater  Shock  Analysis)  Code.  The  FEA  agrees  favorably 
with  the  exact  method  when  calculating  the  transient  response  of 
an  elastic  cylindrical  shell  immersed  in  an  acoustic  media  that 
is  engulfed  by  a  plane  wave. 

Transient  analysis  of  multicomponent  structures 
impacting  rigid  barriers  is  given  by  Reference  [4].  Separate 
analyses  are  performed  on  each  component  considering  distortion 
and  failure  of  joints,  multiple  impacts,  and  e 1  as ti c-pl as t i c 
material  behavior.  Comparison  of  analytical  results  to  test 
data  demonstrates  the  applicability  and  accuracy  of  the  models. 
Next,  the  total  structure  is  analyzed  by  generating  a  new  model 
which  consists  of  the  structure  and  the  individual  component 
models.  The  overall  results  give  good  correlation  between 
analysis  and  test.  (The  program  used  was  SUPER,  a  three- 
dimensional  general  purpose  finite  element  structural  analysis 
code  . ) 

In  Reference  [5],  closed-form  equations  are  devel¬ 
oped  for  comparison  with  test  results.  Resonant  frequency 
results*  are  excellent,  giving  credence  to  the  fact  that  closed- 
form  equations,  if  properly  used,  can  provide  accurate  answers. 

The  author  of  Reference  [6]  compared  three  di  *erent 
closed-form  approximation  techni ques  to  find  an  easy,  low  cost 
method  of  evaluating  vibration  response  characteristics  for 
airborne  optical  packages  .  The  author  states,  "a  finite  element 
anal ys is  of  the  simply  supported  beam  was  conducted  in  order  to 
determine  the  accuracy  of  the  various  analytical  techniques  pre¬ 
sented  here".  The  results  show  that  one  method,  the  generalized 
coordinates  solution,  gave  very  good  estimates  of  the  first  four 
modes.  It  is  important  to  note  that  the  FEA,  and  not  a  test,  is 
used  as  the  check  here.  It  is  also  important  to  realize  that, 
at  times,  a  closed-form  solution  can  provide  acceptable  answers 
at  a  reduced  cost. 
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Reference  [7]  Is  an  excellent  application  of  FEA, 
namely  damping  predictions.  The  authors  used  MSC/NASTRAN  and 
reported  good  results. 

Reference  [8]  is  an  excellent  discussion  on  the 
STAGS  USA  code.  The  authors  say  that  this  code  is  conditionally 
acceptable  for  the  response  of  submerged  shells  to  shock 
loadings  but  improvements  are  needed.  They  have  no  hesitation 
in  recommending  this  program. 

The  following  was  learned  from  References  [1] 

through  [8]: 

o  FEA  modeling,  with  proper  assumptions,  results 

in  excellent  correlation  to  test  data. 

o  FEA  has  been  applied  to  many  problems  but  there 

are  not  many  reports  on  mi croel ectroni  c  package 

analysis . 

o  Closed-form  solutions  can  result  in  accurate, 

cost-effective  solutions. 

o  FEA  analysis  is  sometimes  used  as  a  standard  in 
evaluating  closed-form  solutions  in  lieu  of  test 
data  , 

3.2.3  Personal  References 

Personal  references  include  personal  libraries, 
personal  contacts,  and  personal  experience  in  the  FEA  and  micro¬ 
electronic  package  areas.  Since  the  number  of  personal 
references  Is  so  great,  Table  3  was  generated  to  summarize  the 
f i ndi ngs  . 

A  special  note  is  needed  to  discuss  Reference  [9], 
one  of  the  most  significant  sources  of  information  on  micro¬ 
electronic  package  problems.  This  book  was  discovered  in  a 
personal  library,  and,  after  being  reviewed,  many  interesting 
articles  were  ordered  for  this  study.  These  articles  discussed 
bonding  problems  and  analysis,  gl ass-to-metal -seal  data, 
hermeticity,  package  design,  thermal  design,  fabrication 
techniques,  repai r/rework ,  types  of  packages,  reliability,  (bond 
failure,  circuitry  failure,  encapsulation  failure,  contamination 
and  cleaning),  coatings,  substrates,  and  much  more. 


23 


REMARKS 


This  appears  to  be  a  good  ref.  for 
FEM  analysis  on  bond  joints. 

A  good  theoretical  ref.  for  ther¬ 
mal  stress  analysis. 

Closed  form  versus  Fen  bond  anal¬ 
ysis,  good  correlation. 

A  general  discussion  on  the  con¬ 
trol  of  thermal  stresses  due  to 
CTE  differences. 

Discussion  on  matching  CTE's  for 
silicon  substrates , poiyimice/ 
Kevlar. 

FEM  analysis  compared  to  test  data 
for  bond  joints,  both  good  and 
poor  correlation. 

Band  joint  closed  form  analysis 
appears  to  be  an  excellent  article. 

Thermal  design  and  analysis  of  a 
hybrid  system,  good  analysis  to 
test  correlation. 

Has  a  fortran  transient  computer 
program,  finite  difference,  good 
data. 

Finite  difference  thermal  analysis 
is  used  to  understand  thermal  shock. 
Test  results  for  semi  conductor 
pkgs.,  is  a  good  application  of 
thermal  analysis. 

Glass-to-metal-seal  evaluation, 
general  discussion. 

Glass-to-metal-seal  discussion. 


TABLE  3  -  Discussion  of  References 
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Thermal 


REMARKS 


x  Reliability  of  microcircuit  inter¬ 
connections,  FEM  can  not  be  used 
for  material  corrosion  problem 
discussed  here. 

x  Reliability  can  be  improved  by  re¬ 
ducing  the  number  of  electrical 
joints,  using  more  stable  substrate 
materials,  provide  greater  resist¬ 
ance  to  mechanical  stresses.  They 
also  reccntnend  high  temperature 

operational  testing  for  burning. 

x  Glass-to-metal-seal  information. 

x  Glass-to-metal-seal  information. 

x  Glass-to-metal-seal  information. 

x  The  adhesion  promoter  is  very  im¬ 
portant  for  moisture  protection  and 
at  elevated  temperatures  parylene 
cracked. 

x  Good  article  on  parylene  coating  for 
microcircuits. 

A  good  selection  of  FEM  and  closed 
form  analysis  with  test  data. 

x  Appears  to  be  a  good  article  look¬ 
ing  at  thermal  fatigue  of  bonded 
joints. 

Test  and  analysis  data,  good  corr¬ 
elation. 

x  Hybrid  substrate  bonding  specificat¬ 
ion. 


A  good  ref.  for  thermal  fatigue  in 
bonded  joints,  good  test  to  analysis 
correlation. 


TABLE  3  (cont'd) 


Ttenrial 


Test  Closed  ?oi 
Data  Analysis 


Good  data. 

Old  data, not  all  of  it  is  valid. 

Appears  to  be  valuable  data. 

Good  source  of  bond  joint  analysis 
and  test  for  static  4  fatigue 
loadings. 

Appears  to  be  a  good  article. 

A  good  reference  for  bond  joint 
analysis. 

Good  correlation  for  non  linear 
analysis,  bond  joint. 

Good  correlation,  presents  guide¬ 
lines  for  designing  bond  joints. 

x  Pour  out  of  thirty  were  concerned 
with  analytical  methods  for  sealing, 
good  response  from  questionnaire 
(also  made  many  visits)  the  top 
persistant  failure  modes  in  hybrid 
microcircuits  were  in  descending 
order;  wire  bonds,  defective  semi¬ 
conductor  devices,  contamination, 
moisture  related,  4  leaky  packages. 

x  Good  source  of  hybrid  articles. 

x  See  para.  3.2-3  for  detailed  write¬ 
up. 

x  Test  data  compared  to  finite  diff¬ 
erence  analysis,  good  correlation 

x  Gccd  reference  for  lead  corrosion. 

Appear  co  be  an  excellent  comparison 
between  F&l  4  test , g’-'od  correlation. 


TABLE  3  (cant'd) 


lef .  FEW 
to. _ Analysis 


Test 

Data 


Closed  Form 


Analysis 
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Thermal  design  criteria  for  sub¬ 
strates. 


49 


x 


50 


x 


51 


Thermal  analysis  of  integrated 
circuit  packages. 

x  Thermal  analysis  versus  test  data  for 
chip  carriers,  good  correlation. 

A  good  reference  article. 


52 


53 


x 


54 


55 


56 


57 


x  An  excellent  survey  of  DIP  and  flat- 
pack  mfg. ,  problems  which  can  be  re¬ 
lated  to  hybrids,  #1  problem  was 
lack  of  hermeticity  with  die  bond 
failure  #2. 

A  good  reference  on  isoparametric 
finite  elements. 

x  A  good  source  of  microcircuit  pkg. , 
problems. 

x  An  excellent  reference  on  surface 
treatment  for  bonding  surfaces. 

x  A  very  good  source  of  state-of-the- 
art  packaging  problems. 

An  excellent  collection  of  FEM  app¬ 
lications  &  evaluations,  this 
reference  is  discussed  in  detail  in 
the  text. 


TABLE  3  (cant'd) 
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Reference  [44]  is  deserving  of  some  attention  since 
it  has  some  very  pertinent  data  for  mi c roel ec tron i c  packages. 
Shown  below  is  a  list  of  information  gained: 

o  Commercial  houses  are  claiming  lower  cost  and 
higher  reliability  than  military  houses  (pages  6 
and  7  )  . 

o  Die  rework  and  testing  account  for  75%  of  hybrid 
cost  ( page  21 )  . 

o  S pec i f i ca ti on  relaxation  can  reduce  hybrid  cost 
10%  (may  be  an  application  for  computer 
analysis)  . 

o  Recommenda*  *  ons  were  made  for  deletion  of  the 
following  tests  (page  29): 

Thermal  cycling 
Mechanical  shock 
SEM  Inspection 
Lead  bond  testing 
Centri f  uge 
Die  shear  test 

o  Only  67%  of  those  who  responded  to  the  survey 

used  thermal  analysis.  No  mention  was  made  of 
stress  and  dynamics  which  means  they  ignore  it 
or  use  engineering  judgement  (page  37). 

Reference  [57]  is  a  collection  of  many  FEA  articles 
and  is  also  deserving  of  some  special  discussion.  Below  are 
discussions  on  pertinent  papers: 

(1)  Paper  5,  Volume  1  "Practical  Aspects  of  the 

Finite  Element  Method,"  Hoggenmacher ,  G.G.,  and  Lahey,  R.S.  The 
important  conclusions  or  considerations  are  listed  below: 

o  The  accuracy  of  some  membrane  elements  may 

deteriorate  at  higher  aspect  rati. os.  An 

arbitrary  system  of  elements  covering  a  given 

region  may  not  result  in  useful  data;  sound 
judgment  is  needed. 
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o  Model  integrity  can  be  verified  by 

1.  Line-by-line  check  of  the  input  data 

2.  Computer  plots  to  detect  ill-shaped 
structures 

3.  A  "common  sense"  review  to  detect 
unreasonable  results 

4.  Ecuilibrium  checks  on  all  nodes 

5.  Constraint  compatibility  check 

o  The  key  point  is  to  understand  the  limitations 
of  the  elements  you  are  using  in  order  to 
achieve  the  best  possible  model.  Then  apply  all 
possible  checking  techniques. 

(2)  Paper  [9],  "Utilization  of  Isoparametri c  Shell 
Elements  in  Solution  of  Practical  Problems,"  Citipitiogl u,  E., 
Nicholas,  V.T.,  and  Ecer,  A.  This  is  an  excellent  article  on 
the  practical  usage  of  i  soparametric  shell  elements:  both 

middle  surface  shell  (MSS)  and  solid  type  shell  (STS). 
Instructive  examples  are  presented  which  show  the  pitfalls  in 
using  these  element  types.  One  example  presented  is  the 
comparison  between  analysis  and  experiment  for  a  steel  wheelwell 
from  the  front  end  of  an  automobile.  The  wheelwell  was  stamped 
from  0.033  inch  thick  steel  into  a  complex  geometric 
configuration.  The  thickness  variation  caused  by  the  stamping 
process  varied  from  0.030  to  0.035  inches.  There  were 
thirty-nine  holes  of  different  sizes  with  many  types  of  ribs. 
The  results  for  the  first  natural  frequency  are  summarized 
below: 

Number  of  Elements  in  Model:  78  142  204  289  520 

%  Error:  +145%  +94%  +51%  21%  7% 

This  example  demonstrates  the  analysis  by  an 
iterative  modeling  method.  Basically,  a  coarse  model  is  chosen 
initially,  and  then  it  is  systematically  improved. 

Many  other  interesting  articles  are  presented  in 
Reference  [57],  some  of  the  more  pertinent  ones  are  listed 
below: 
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Paper  11: 

Paper  12: 

Paper  16 : 
Paper  19: 


"Higher  Order  Versus  Lower  Order  Elements;  Economics 
and  Accuracy,"  MacNeal,  R.H. 

"Element  Evaluation  -  A  Set  of  Assessment  Points  and 
Standard  Tests,"  Robinson,  J. 

"Thin  Shell  Isoparametri c  Elements,"  Berkovic,  M. 
"Overview  and  Evaluation  of  Some  Versatile  General 
Purpose  Finite  Element  Computer  Programs," 
Frederikson,  8.,  Mackerle,  J. 


3  .3  Trips 

Two  trips  were  made  for  the  survey  portion  of  this 
study.  One  trip  was  a  conference  with  a  large  representati on  of 
microelectronic  package  manufacturers  and  users.  The  other  trip 
was  to  a  microelectronic  package  manufacturer.  Many  key  points 
were  learned  on  these  trips;  the  results  are  summarized  below 
based  on  these  two  trips  and  the  experience  of  the  many 
companies  polled. 

o  Most  companies  do  not  use  FEA  to  assess  the 

structural  and  thermal  integrity  of 
microelectronic  packages.  Most  large  companies 
use  the  finite  element  method  for  analysis  of 

structures  . 

o  FEA  for  m i c roe  1 e c t r on i c  packages  has  great 

potential  . 

o  Very  few  recommend  the  "full  blown"  screening  in 
MIl-STD-883.  Most  utilize  an  abbreviated 

version  that  has  worked  for  them. 

o  Each  company  has  their  own  preference  for  finite 
element  computer  programs,  but  the  most  widely 
used  are  STARDYNE,  ANSYS,  and  NASTRAN. 

3  .4  Analysis  and  Test  Cost  Comparison  for  Qualification 

Veri f i cati on 

The  purpose  here  was  two-fold;  1)  to  determine  which 
typical  mi croel ectroni c  military  environmental  tests  can  be 
replaced  by  finite  element  analysis  and  2)  the  relative  cost 
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comparisons  between  closed-form  analysis,  finite  element 
analysis  and  testing. 

3.4.1  Finite  Element  Analysis  Applied  to  Mi c roel ec tron i c 

Package  Qualification  Verification 

The  purpose  here  was  to  determine  which,  if  any, 
military-specified  test  environments  for  electronic  packages 
might  lend  themselves  to  a  structural  or  thermal  analysis  prior 
to  a  test  being  performed  on  an  actual  piece  of  hardware. 
Obviously,  schedules  and  budgets  could  benefit  from  the  use  of 
analysis  to  detect  design  shortcomings  prior  to  expensive  and 
time  consuming  test  sequences.  Therefore,  the  test  methods  of 
four  different  military  specified  tests  were  reviewed  to 
determine  the  possibility  of  supplementing  the  test  secuences 
with  either  a  finite  element  analysis  (structural  or  thermal)  or 
a  finite  difference  analysis  (thermal  only). 

The  four  military  test  sequences  which  were  reviewed 
included  MIL-STD-883  (Test  Methods  and  Procedures  for 
Microel ectronics ) ,  MIL-STD-810  (Environmental  Test  Methods), 
MIL-STD-202  (Test  Methods  for  Electronic  and  Electrical 
Component  Parts),  and  MIL-M-38510  (Digital  Microcircuits).  For 
each  of  these  four  test  sequences,  the  test  methods  which  could 
be  suitable  for  a  finite  element  and/or  finite  difference 
analysis  were  noted. 

3. 4. 1.1  Applicable  Test  Methods 

Tables  4  through  7  indicate  the  test  methods 
suitable  for  analysis  from  each  of  the  four  military 
specifications  considered.  Also  inlcuded  is  a  brief  description 
of  the  type  of  analysis  which  might  be  used  to  simulate  the 
various  test  methods. 

MIL-STD-883 

Table  4  summarizes  those  test  methods  of 

MIL-STD-883B  which  were  considered  suitable  for  analysis: 
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Table  4.  MIL-STD-883B  Test  Methods 

Method  No. 

1010  .2 
1011  .2 
1012 

1015  .2 

1016 
2001  .2 
2002  .2 

2004  .2 

2005  .1 

2007  .1 

2011  .2 

2019  .1 

As  shown  in  Table  4,  twelve  of  the  test  methods  in 
MIL-STD-883  were  considered  suitable  for  structural  and  thermal 
analysis.  Methods  1010.2  and  1011.2  could  use  a  finite  element 
heat  transfer  program  to  define  component  temperatures ,  which 
could,  in  turn,  be  used  to  define  the  loading  for  a  thermal 
stress  analysis  of  the  same  system.  Methods  1012,  1015.2,  and 
1016  could  use  either  a  finite  element  heat  transfer  program,  or 
a  finite  difference  program,  to  solve  for  specific  component 
temperatures  based  on  the  boundary  conditions  specified  in  the 
test  methods.  In  each  of  the  above-mentioned  cases,  analysis 
could  be  used  to  predict  a  failure  or  to  spot  a  marginal  design 
(excessive  junction  temperatures,  excessively  large  thermal 
stresses,  etc.)  prior  to  the  test  sequence. 

From  a  structural  point  of  view,  test  methods 
2001.2,  2004.2,  2011.2,  and  2019.1  could  be  analyzed  using  the 
static  analysis  capabilities  of  a  finite  element  program. 


Desc  ri pti on 


Temperature  Cycling 
Thermal  Shock 
Thermal  Charactcri sti cs 
Burn-In  Test 
Life/Reliability  Tests 
Constant  Acceleration 
Mechanical  Shock 
Lead  Integrity 

Vibration  Fatigue 

Vibration,  Variable 
Fre  quency 
Bond  Strength 

Die  Shear  Strength 


ilicable  Analyses 


Thermal  ,  thermal  stress 

Thermal,  thermal  stress 

Thermal 

Thermal 

Thermal 

Static  stress 

Dynamic  stress  (shock) 

Static  stress,  static 

nonl i near  stress 

Dynamic  stress 

(  s  i  nusoi  dal ) 

Dynamic  stress 
( si nusoi dal ) 

Static  stress,  static 
nonlinear  stress 
Static  stress 


Methods  2001.2  and  2019.1  would  be  relatively  straight-forward 
statics-type  problems;  however,  methods  2004.2  and  2011.2  would 
require  the  program  to  have  nonlinear  analysis  capabilities. 
These  nonlinear  capabilities  would  be  required  in  order  to 
s ati sf actori ly  analyze  the  lead  bending  of  method  2004.2  (large 
deflection  theory  with  plastic  deformation),  and  to  analyze  the 
various  bonds  of  method  2011.2  (a  common  bond,  solder,  tends  to 
plastically  relieve  itself  at  stresses  exceeding  about  500  psi). 

Similarly,  test  methods  2002.2,  2005.1,  and  2007.1 
could  be  analyzed  using  a  structural  finite  element  program  with 
dynamic  analysis  capabilities.  For  both  the  shock  and  the 
sinusoidal  vibration  environments,  the  finite  element  program 
would  be  used  to  calculate  the  system’s  natural  frequencies,  and 
to  then  calculate  the  dynamic  response  stresses  based  on  user- 
provided  transmissibil i ty  assumptions. 

Using  the  stresses  calculated  by  the  finite  element 
analyses  of  the  various  mechanical  tests,  marginal  or  inadequate 
designs  could  be  spotted  (excessively  large  stresses  or 
unsati sf actory  fatigue  coefficients)  prior  to  expensive  testing 
se  ouences . 

MIL-STD-810 

Table  5  summarizes  those  test  methods  of 
MIL-STD-810C  which  were  considered  suitable  for  analysis. 

*As  shown  in  Table  5,  thirteen  of  the  test  methods  in 
MIL-STD-810  were  considered  suitable  for  structural  and  thermal 
analysis.  Methods  500.1,  504.1  ,  517.2,  and  518.1  would  require 
a  combination  of  "static"  structural  analysis  ( di f f erenti al 
press ure/venti ng  analysis),  heat  transfer  analysis  (varying 
convective  and  conductive  characteristics)  ,  and  thermal  stress 
analysis.  Methods  50  1  .  1  ,  502.1  ,  503.1  ,  and  50  5  .  1  would 
similarly  require  a  heat  transfer  analysis  and  a  thermal  stress 
analysis  to  determine  the  effects  of  the  specified  thermal 
1 oads  . 
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Table  5.  MIL-STD-810C  Test  Methods 


Method  No. 

Oescri pti on 

Applicable  Analyses 

500  .1 

Low  Pressure  (Altitude) 

Static  stress , thermal 

501  .1 

High  Temperature 

The rmal ,  thermal 

stress 

502 .1 

Low  Temperature 

Thermal  ,  thermal 

stress 

503  .1 

Temperature  Shock 

The  rmal ,  thermal 

stress 

504  .1 

Temperature-Al titude 

Static  stress, 

thermal ,  thermal 

stress 

505  .1 

Solar  Radiation 

Thermal ,  thermal 

stress 

513.2 

Acceleration 

Static  stress 

514  .2 

Vibration 

Dynamic  stress 

(sine  A  random) 

515  .2 

Acoustical  Noise 

Dynamic  stress 

516  .2 

Shock 

Dynamic  stress 
( shock ) 

517  .2 

Space  Simulation 

Static  stress, 
thermal ,  thermal 

stress 

518  .1 

Temperature-Humidity- 

A1  ti tude 

Static  stress, 

thermal ,  thermal 

stress 

519  .2 

Aircraft  Gunfire 

Vibration 

Vibration  (sine 
superimposed  over 

random) 

From  a  structural  viewpoint,  method  513.2  would 
require  a  static  analysis,  while  methods  514.2,  515.2,  516.2, 
and  519.2  would  all  require  some  type  of  dynamic  stress  analysis 
(either  sinusoidal  or  random  vibration,  acoustical  noise,  or 
shock).  For  method  519.2,  it  would  be  necessary  to  use  super- 
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position  to  determine  the  effects  of  the  combined  sinusoidal  and 
random  vibration  environment. 

MIL-ST0-202 

Table  6  summarizes  those  test  methods  of 
MIL-STD-202E  which  were  considered  suitable  for  analysis. 

Twelve  of  the  test  methods  in  MIL-STD-202  were 
considered  suitable  for  structural  and  thermal  analysis.  Method 
105C  would  require  a  static  stress  analysis  to  account  for  the 
behavior  of  any  pressurized  assemblies,  and  a  heat  transfer 
analysis  to  account  for  the  varying  convective  and  conductive 
characteristics.  Methods  107D,  108A  and  210A  would  all  require 
a  combined  heat  transfer  and  thermal  stress  analysis  to 
determine  the  effects  of  the  specified  thermal  loads. 


Method  No. 

Table  6  .  MIL-STD-202E  Test 

Descri pti on 

Methods 

A ppl i c abl e  Analyses 

105C 

Barometric  Pressure 

( reduced) 

Static  stress .thermal 

1070 

Thermal  Shock 

Thermal ,  thermal 

stress 

108A 

Life  (at  elevated 

Thermal ,  thermal 

ambient  temp) 

stress 

20 1A 

V i brati on 

Dynamic  stress 

(si nusoi dal ) 

203E 

Random  Drop 

Dynamic  stress 
(repeated  shock) 

204C 

Vibration,  High 

Dynamic  stress 

Fre  quency 

( si n  usoi dal ) 

207A 

High  Impact  Shock 

Dynamic  stress 
( shock ) 

2 1 0A 

Resistance  to  Soldering 

Thermal ,  thermal 

Heat 

stress 

21 1A 

Terminal  Strength 

Static  stress,  static 

nonl i near  stress 

212A 

Accel erati on 

Static  stress 

213B 

Shock  (specified  pulse) 

Dynamic  stress 

( shock ) 

214 

Random  Vibration 

Dynamic  stress 

( random) 
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From  a  structural  point  of  view,  methods  201A,  2038, 
204C ,  207A,  2138,  and  214  would  all  require  some  type  of  dynamic 
stress  analysis  --  either  sine  vibration,  random  vibration,  or  a 
specified  shock  or  series  of  shocks.  Finally,  methods  2 1 1 A  and 
212A  would  require  a  static  stress  analysis  to  determine  the 
resulting  stresses  .  However,  it  would  be  necessary  to  use  a 
non’inear  program  for  method  2 1 1 A :  the  "Terminal  Strength" 
method  involves  the  plastic  deformation  of  terminal  leads. 

MIL-M-38510 

Table  7  summarizes  the  test  procedures  of 
MIL-M-38510  which  were  considered  suitable  for  analysis: 

Table  7.  MIL-M-38510  Test  Procedures 


Descri pti on 
Stabilization  Bake 


Temperature  Cycling  and/or  Shock 
Constant  Acceleration 


Bond  Strength 


Lead  Integrity 

Steady  State  Life  Tests 
Die  Shear  Strength  Test 
Barometric  Pressure 
Intermittent  Life 


Applicable  Analyses 
Thermal,  thermal  stress 
Thermal,  thermal  stress 
Static  stress 
Static  stress,  non¬ 
linear  static  stress, 
thermal  stress 
Static  stress,  non¬ 
linear  s tati c  stress  , 
The  rmal 
Static  stress 
Static  stress,  thermal 
Thermal 


Nine  of  the  test  procedures  in  MIL-M-38510  were 
considered  suitable  for  structural  and  thermal  analysis.  The 
"Stabilization  Bake"  and  the  "Temperature  Cycling  and/or  Shock" 
tests  would  require  a  combined  heat  transfer  and  thermal  stress 
analysis  to  determine  the  effects  of  the  specified  thermal 
loads.  In  addition,  the  "Steady  State  Life  Tests"  (accelerated 
and  normal  tests),  the  "Intermittent  Life",  and  the  "Barometric 
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Pressure"  methods  would  all  require  a  heat  transfer  analysis  of 
the  system. 

From  a  structural  viewpoint,  the  "Constant 
Acceleration",  the  "Bond  Strength",  the  "Lead  Integrity",  the 
"Die  Shear  Strength  Test",  and  the  “Barometric  Pressure"  methods 
could  all  benefit  from  a  static  stress  analysis.  Also,  it  would 
be  necessary  to  use  a  nonlinear  program  to  analyze  the  stresses 
in  the  "Bond  Strength"  and  "Lead  Integrity"  methods  (solder 
bonds  plastically  relieve  themselves,  and  leads  are  plastically 
deformed  in  these  tests). 

3  .4  .2  11  Front-End"  Analysis  Cost  Savings 

The  purpose  of  this  section  is  to  compare  the 
relative  costs  of  closed-form  analysis  and  finite  element 
analysis  supplementing  qualification  testing  of  microel ectroni  c 
packages.  Obviously  not  all  qualification  tests  can  be 
replaced;  however,  based  on  the  previous  discussion,  let  us 

assume  that  the  following  tests  can  be  replaced  by  analysis; 
shock  (0.6  to  11.0  msec,  up  to  100  G's),  acceleration  (up  to 

30,000  G's),  variable  frequency  vibration  (20  to  70  G's), 

f  tigue  vibration  (20  G's  to  50  G's),  and  random  vibration  (up 
to  50  Grms).  With  this  in  mind,  the  following  assumptions  are 
given  for  an  example  in  cost  savings: 

1.  Hybrid  cost  -  $50Q/each 

2.  A  military  program  has  a  requirement  for  1,000 
hybrids  per  MIL-STD-883. 

3.  Compliance  to  MIL-STD-883,  on  a  recent  Harris 

program  with  11  hybrid  types,  required  23%  more 
non  de  1  i  v  e  r  ab  1  e  hybrids  to  be  used  for  5008 
testing.  Therefore,  23%  will  be  assumed  for 
this  example  as  the  number  of  extra  hybrids  for 
MIL-STD-883  compliance. 

4.  230  hybrids  are  required  for  the  full  compliance 

to  MIL-STD-883.  Part  of  this  number  of  test 
samples  Is  required  for  dynamic  testing,  let  it 

be  assumed  that  100  hybrids  can  be  deleted  if 


w 


analysis  is  substituted  for  dynamic  testing. 
This  is  only  an  assumption  and  will  vary  with 
each  individual  program.  The  re fore,  a  cost 
reduction  of  $100  x  $500  or  $50,000  would  be 
saved  if  analysis  was  used  to  supplement  the 
tests.  Below  are  typical  current  costs  for 
testing  and  analysis  based  on  engineering 
experience  at  many  different  companies. 

5.  Dynamic  Testing  cost  is  $3,200 

6.  Finite  element  analysis  cost  is  $7,000 

7.  Closed-form  analysis  cost  is  $3,000 

Table  8  presents  a  summary  of  the  cost  comparisons 
between  analysis  and  test.  It  is  seen  that  a  substantial  amount 
can  be  realized  based  on  a  company's  overhead  structure.  The 
savings  for  this  hypothetical  case  is  $46,200  for  finite  element 
analyses  and  $50,200  for  closed-form  analysis.  Now  this  is,  of 
course,  assuming  that  all  of  the  dynamic  tests  are  deleted  which 
is  not  recommended.  What  is  recommended  is  to  perform  an 
analysis  using  finite  element  analysis,  and  develop  a  relative 
comparison  of  resulting  stress  levels  for  each  dynamic  test. 
Then  the  worst  case  tests  can  be  selected  and  the  equipment 
subjected  to  these  tests.  It  is  estimated  that  this  procedure 
would  result  in  a  real  cost  savings  of  roughly  $25,000  to 
$30,000,  not  an  insignificant  amount.  Additional  cost  savings 
can  be  realized  by  using  pre  and  postprocessor  software 
programs,  and  by  utilizing  Computer  Aided  Design/Computer  Aided 
Manufacturing  (CAD/CAM)  techniques. 


Table  8.  Cost  Comparison  for  Analysis  and  Test  of  a  Hybrid 

Component 


Dynamic  Testing 
Compliant  to 
MIL-STD-883 


Finite  Element 
Analysis  Substituted 
for  Oyanmic  Testing 


Closed-Form  Analysis 
Substituted  for 
Dynamic  Testing 


$53,200 


S7 ,000 


$3,000 
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3 .4 .3  $  ummary 

The  previously  discussed  tables  list  those  test 
methods  of  MIL-STD-883,  MIl-STD-810,  MIL-STD-202,  and 
MIL-M-38510  which  we  consider  to  be  suitable  for  a  finite 
element  structural  or  thermal  analysis,  or  for  a  finite 
difference  thermal  analysis.  It  is  believed  that  the  use  of 
analysis  in  conjunction  with  testing  will  benefit  schedules  and 
budgets;  design  shortcomings  may  be  identified  and  solved  prior 
to  expensive  testing  seouences.  It  has  been  shown,  by  use  of  an 
example,  that  cost  savings  from  525,000  to  $50,200  can  be 
realized  by  substituting  analysis  for  dynamic  testing  per 
MIL-STD-883  on  hybrids. 

3  .5  General  Summary  for  Techicai  Assessment  (Phase  1) 

The  more  germane  findings  for  the  technical 
assessment  are  presented  below: 

o  Although  FEA  is  primarily  a  tool  for  designers 
of  large  complex  structures,  a  great  potential 
exists  if  applied  to  mi c roel ec tron i c  packages, 
o  By  using  FEA  for  early  warnings  of  potential 

failures  the  mechanical  engineer  can  play  a 
significant  role  in  assuring  or  improving  the 
reliability  of  microelectronic  systems, 
o  The  survey  showed  that  the  areas  giving 

microelectronic  packages  the  most  problems  were 
hermeticity,  broken  or  corroded  internal  wires, 
and  broken  wi rebonds. 

o  FEA  can  provide  accurate  answers  if  proper 

assumptions  are  made  and  if  one  understands  the 
finite  element  theory. 

o  Commercial  hybrid  houses  claim  lower  cost  and 

higher  reliability  than  similar  products 
produced  at  military  hybrid  houses, 
o  Most  companies  do  not  use  FEA  to  assess  the 

structural  and  thermal  integrity  of  micro¬ 
electronic  packages. 
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Tne  "full  blown"  screening  in  MIL-STD-883  is  not 
recommended  by  most  of  the  companies  polled. 
Instead,  an  abbreviated  version  that  works  for  them 
is  utilized. 

o  Closed-form  solutions  can  result  in  accurate, 
cost-effective  solutions. 

o  Considerable  cost  savings  can  be  realized  by 
substituting  FEA  for  MIL-STD-883  screening  tests. 
An  example,  showed  cost  savings  from  $25,000  to 
$50,000  for  a  lot  of  1000  hybrids  with  11  hybrid 
types  . 

4 -0  ANALYSIS  METHOOS  (PHASE  2) 

4 .1  Introducti on 

This  phase  of  the  study  consisted  of  taking  the 
information  gained  in  Phase  1  in  the  form  of  programs,  closed-form 
solutions,  test  data,  and  problems,  and  deciding  which  to  use.  The 
computer  programs  studied  were  those  most  pertinent  to  the 
mi c roel ectroni c  problems.  Sample  problems  were  then  developed  for 
evaluating  these  computer  programs. 

The  topic  of  general  purpose  programs  has  been  one  of 
considerable  interest  since  the  power  of  the  finite  element  first 
became  recognized.  The  capabilities  of  today's  programs  are 
numerous  and  varied  --  many  programs  have  identical  capabilities  in 
several  areas  of  analysis.  Several  programs  provide  excellent 
ability  to  solve  problems  in  a  limited  area  of  structural 
mechanics,  whereas  others  provide  the  ability  to  solve  problems  in 
many  areas  of  structural  mechanics.  This  makes  the  evaluation  of 
general  purpose  programs  a  very  difficult  task.  There  are  so  many 
parameters  involved  that  it  is  almost  impossible  to  do  a  parameter 
study  with  all  parameters  included.  Furthermore,  many  of  the 
parameters  involved  are  of  a  "subjective"  nature  and  are  difficult 


to  Quantify  [69].  Therefore,  the  selection  of  a  general  purpose 
finite  element  program  reouires  the  user  to  identify  his  analysis 
needs  --  it  also  usually  involves  some  degree  of  personal  prefer¬ 
ence.  It  is  important  that  the  user  not  restrict  himself  to 
selecting  one  program  --  this  is  especially  true  if  the  user  has 
wide-ranging  analysis  needs. 

4 .2  Criteria  for  Program  Selection 

Listed  below  are  the  requirements  for  program  selection 
[58].  An  ideal  program  with  analysis  capabilities  in  all  areas 
would  possess  all  the  requirements  below.  However,  no  program  has 
"complete”  analytical  abilities.  A  good  program  would,  however, 
possess  most  of  these  requirements: 

1.  It  should  have  proven  capability  for  both  linear 
static  and  linear  dynamic  response,  including 
computation  of  natural  frequencies.  Desirable 
capabilities  include  (a)  nonlinear  static  response, 
(b)  nonlinear  dynamic  response,  (c)  stability 
(buckling),  both  static  and  dynamic,  (d)  thermal 
loading,  and  (e)  heat  transfer  analysis. 

2.  A  good  library  of  planar  and  three-dimens ional 
elements  should  be  available.  It  should  include: 

a.  One-dimensional  straight  elements  for  axial 
1  o  a  d  i  n  g  , 

b.  One-dimensional  straight  and  curved  elements, 
capable  of  resisting  forces  and  moments  in 
three  di mens  i  or.s  , 

c.  Plate  and  shell  elements  of  various  shapes, 
both  for  membrane  and/or  bending  loading, 

d.  Three-dimensional  elements  of  various  shapes 
and  numbers  of  nodes, 

e.  Axisymmetric  elements  and  axisymmetric  thin  and 
thick  shell  elements,  and 
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f.  Geometrically  nonlinear  elements,  such  as  gaps 
and  tension-only  or  compression-only  elements. 

3.  The  capability  for  handling  linear  and  nonlinear 
materials  . 

4.  It  should  provide  accurate  answers. 

5.  It  should  solve  problems  quickly,  with  a  minimum 
cost . 

6.  It  should  be  readily  available  and  have  well- 
documented  User's,  Examples,  and  Theoretical 
Man ua 1 s  . 

7.  Competent  technical  assistance  must  be  available, 
either  from  the  host  computer  center,  from  the 
program  developers,  or  from  the  user's  group. 

The  above-mentioned  list  of  program  attributes  are 
necessary  for  program(s)  performing  the  entire  gamut  of  finite 
element  analyses.  For  the  specific  case  of  finite  element  analysis 
of  electronics  equi  pmen t  and  mi croe 1 ec tron i cs  ,  in  particular,  the 
following  program  prerequisites  are  offered.  The  structural 
loadings  should  include  statically  applied  forces,  thermal 
loadings,  sinusoidal  vibration,  random  vibration,  shock  loading, 
and  other  generalized  transient  dynamic  loadings.  Therefore,  a 
program  must  have  capabilities  for  solving  linfear  statics  and 
dynamics  problems  with  a  wide  variety  of  loading  conditions. 

When  temperature  distributions  are  necessary  for  a 
thermal  stress  analysis,  steady  state  and  transient  heat  transfer 
analyses  become  desirable  program  capabilities.  Finally,  many 
electronics  equipment  and  materials  ( i  .  e  .  ,  solder)  exhibit 
nonlinear  behavior.  This  makes  the  ability  to  solve  nonlinear 
statics  and  dynamics  problems  a  desirable  program  attribute.  The 
area  of  nonlinear  material  behavior  (especially  that  of  solder)  is 
a  promising  area  for  future  work  and  needs  more  investigation. 
Typically,  microelectronics  problems  involve  models  of  the 
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sma 1 1  - to-med i urn  size  range.  This  precludes  the  necessity  for 
solving  large  models  with  huge  numbers  of  elements. 

4  .3  Preliminary  List  of  Programs  Reviewed 

Reviews  of  the  majority  of  well-known  programs 
available  in  the  United  States  and  Europe  were  completed.  Below  is 
the  initial  list  of  21  programs  that  were  reviewed  during  this 
study.  Obviously,  the  best  way  to  evaluate  a  program  is  to  use  it 
constantly  over  a  long  period  of  time.  This  would  involve  a 
substantial  effort,  and  was  considered  beyond  the  scope  of  this 
study.  Based  on  their  capabilities,  12  programs  from  this  initial 
list  were  then  eliminated  from  further  consideration  and  an  interim 
list  of  nine  programs  was  made.  This  interim  list  included  most  of 
the  widely  known  general  purpose  programs  .  The  programs  in  the 
interim  list  were  reviewed  in  detail  and  five  of  these  programs 
were  then  eliminated,  forming  a  final  matrix  of  four  programs.  The 
programs  in  the  final  matrix  were  further  reviewed,  including 
benchmark  computer  runs  of  these  programs. 

The  programs  reviewed  are  available  on  a  wide  variety 
of  computer  network  systems.  This  general  sense  of  wide  avail¬ 
ability  makes  it  relatively  easy  for  the  analyst  to  use  the  program 
of  his  choice.  Program  capabilities  are  constantly  being  upgraded 
and  increased.  This  makes  it  virtually  impossible  to  easily 
compare  program  capabilities.  Therefore,  the  program  comparisons 
were  based  on  available  information  that  was  current  at  the  time  of 
this  writing. 

The  initial  list  of  programs  is  shown  below: 

1.  ABAQUS 

2.  ADINA 

3.  ANSYS 

4.  ASKA 

5.  COSA 

6.  EASE2 
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7.  FINITE 

8.  MARC-CDC 

9.  MINIELAS 

10.  NASTRAN 

11 .  NISA 

12.  NONSAP 

13.  PISCES 

14.  SAP  IV 

15 .  SES AM-69 

16.  SPAR 

17  .  STAGS 

18  .  STAROYNE 

19  .  STRUDL  1 1 

20.  STRUOL  DYNAL 

21 .  TEXGAP  3-D 


4 .4  Development  of  Interim  Program  List 

The  following  is  a  list  of  the  12  programs  eliminated 
from  further  consideration  after  reviewing  the  original  21 
programs.  After  each  program  is  a  brief  synopsis  of  program 
capabilities,  along  with  reasons  for  that  program's  elimination 
from  further  consideration: 

1.  ASKA  -  This  is  a  fairly  large  program  that  can 
perform  both  linear  and  nonlinear  analyses  [69]. 
It  was  developed  in  Germany  [58]  and  has 
capabilities  in  dynamics  analysis.  It  has  an 

extensive  element  library,  good  mesh  generation, 
and  is  especially  effective  with  substructuring. 
However,  it  is  a  very  large  system  requiring 
considerable  user  experience  [58].  It  is  not 

readily  available  in  the  United  States.  A  recent 
study  [78]  has  further  evaluated  this  program. 
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ASKA  does  not  appear  on  any  major  computer  network 
in  the  U.S.  It  appears  not  to  have  random 
vibration  capabilities.  This  would  make  this 
program  undesirable  for  mi croe 1 ec t ron i c  analysis. 

2.  C QSA  -  This  program  appears  to  have  been  developed 
in  Germany.  It  has  a  good  element  library  and 
appears  suitable  for  dynamic  analysis.  Very  little 
information  was  available  on  this  program  [58]. 
Performing  only  dynamic  analysis  makes  this 
program's  use  very  limited  for  microelectronic 
analysis  . 

3.  EASE2  -  This  is  a  program  limited  to  linear  statics 
and  dynamics  analysis  [59].  It  seems  to  be 
well-suited  to  large  problems.  Little  information 
was  found  on  this  program.  It  does  not  possess 
heat  transfer  capabilities.  However,  it  appears  to 
have  limited  dynamic  analysis  features.  A  recent 
study  [79]  revealed  that  EASE2  does  not  have  random 
vibration  capab i 1 i ti es  .  It  is,  therefore,  not  a 
desirable  program  for  mi c roel ec tron i cs  analysis. 

4.  MINIELAS  -  This  is  part  of  the  EUAS  system  [58]  and 
is  limited  to  dynamic  analysis.  It  is  a  special- 
purpose  program  designed  for  only  random  vibration 
analysis.  It,  therefore,  has  very  limited  capabil¬ 
ities  and  is  not  well  suited  for  microelectronic 
analys i s  . 

5.  NONSAP  -  This  program  is  limited  to  static  and 

dynamic  analysis  of  shell  structures  [58],  It  is 

the  forefunner  of  the  ADINA  program,  to  be 
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discussed  later.  Since  microelectronic  packages 
are  not  shell  structures,  this  program  is  not 
s ui table  . 

6.  PISCES  -  This  is  a  program  that  solves  linear  and 

nonlinear  problems.  It  is  especially  strong  where 
nonlinear  materials  are  being  analyzed  [58]. 
However,  it  employs  the  finite  difference  method  of 
analysis,  rather  than  the  finite  element  method. 
The  finite  difference  method  employs  nodes,  as  does 
the  finite  element  method.  However,  the  finite 

difference  method  requires  all  the  information 
concerning  node  connectivity  (i.e.,  stiffness)  to 
be  fully  defined.  The  finite  element  method  uses 
model  geometry,  element  properties,  and  material 
properties  to  compute  these  stiffnesses.  The 
result  is  that  structural  analysis  via  finite 
elements  is  generally  cheaper  than  that  done  with 
finite  differences  [79].  Structural  analysis  with 
finite  differences  has  very  limited  applications 
and  would  not  be  cost-effective  when  analyzing 
mi croel ectronics  . 

7.  SESAM-69  -  It  was  developed  in  Norway  [58]  and  has 
applications  primarily  designed  for  super  elements. 
It  has  limited  nonlinear  and  dynamics  capabilities 
[69],  This  program  has  recently  been  evaluated 
[78],  It  appears  to  be  a  s pec i al ty- type  program 
with  strong  fracture  mechanics  capabilities.  It 
does  not  appear  to  have  random  vibration 
capabilities,  nor  can  prescribed  accelerations  be 
applied.  This  would  make  this  program  undesirable 
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for  microelectronic  analysis.  A  more  current 
version  of  this  program  is  SESAM-80. 


8.  SPAR  -  This  is  a  series  of  small  interactive 
programs  for  linear  static  and  dynamic  analysis 
[60].  It  appears  to  have  limited  capabilities  and 
does  not  compare  well  with  the  larger  general 
purpose  programs.  Since  this  program  concentrates 
on  beam-like  structures,  it  is  unsuitable  for 
mi croel ectron i c  analysis. 

9.  STAGS  -  This  program  is  limited  to  linear  and 
nonlinear  shell  analysis  [60].  It  is  similar  to 
ANSYS  [61]  and  is  a  good  program  for  this  limited 
application.  It  is  also  a  finite  difference 
program.  The  finite  difference  comments  relating 
to  the  PISCES  program  also  apply  to  STAGS.  Since 
STAGS  concentrates  on  shell  structures,  it  is  not 
suitable  for  microelectronic  analysis. 

10.  STRUDL  II  -  This  is  a  program  developed  for  the 
civil  engineering  community  [62].  It  has  limited 
element  properties  and  loadings  with  poor  mesh 
generation  capability  [58].  It  was  a  very  popular 
program  in  the  building  industry.  However,  further 
developments  on  the  program  were  stopped,  due  to 
lack  of  support.  Me  Donne  1 1 -Do ugl  as  has  developed 
and  marketed  its  own  version,  and  its  use  appears 
to  be  declining  [58].  Since  this  program  concen¬ 
trates  on  civil  engineering  structures  (especially 
frameworks  [79]),  it  is  not  suitable  for 
mi croel ec tron i c  analysis. 
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11.  STRUDL  DYNAL  -  This  is  the  proprietary  version  of 
STRUDl  II  with  emphasis  on  dynamic  analysis  [58]. 
All  the  previous  comments  on  STRUDL  II  apply  here, 
except  that  its  capabilities  have  been  improved  by 
Me  Donnel 1 -Dougl as  .  However,  it  suffers  from 
several  deficiencies,  including  no  transient 
analysis  capability  using  direct  integration.  This 
makes  STRUDL  DYNAL  unsuitable. 

12.  TEXGAP  3-D  -  This  program  was  developed  by  the  Air 
Force  and  is  limited  to  linear  and  nonlinear 
statics  analysis  [61].  It  has  an  element  library 
oriented  towards  continuum  analysis.  It  is  a  good 
program  for  performing  thermal  stress  analysis. 
This  program  appears  to  be  excellent  when  solving 
statics-only  problems.  Since  this  program  cannot 
solve  dynamics  problems,  its  use  would  be  very 
limited  when  performing  microelectronics  analysis. 


The  above  programs  were  then  reviewed  in  detail  and 
five  were  eliminated  from  further  cons i derati on ,  forming  a  final 
matrix  of  four  programs.  These  eliminated  programs  are  summarized 
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below,  along  with  reasons  for  their  elimination.  The  capabilities 
of  the  interim  list  of  nine  programs  are  considered  "very  good." 
This  made  it  difficult  to  eliminate  five  of  these  programs  from 

further  consideration.  The  analyst  should  consider  using  one  of 
the  five  eliminated  programs  when  his  specific  analysis  needs 
correspond  with  the  specific  strengths  of  one  of  these  programs. 
For  certain,  specific,  uncommon  problems,  the  five  eliminated 

programs  possess  capabilities  that  are  considered  "excellent"  when 
performing  microelectronic  analysis. 

1.  AD  I NA  -  This  is  a  program  that  performs  linear  and 
nonlinear  static,  dynamic,  and  heat  transfer 
analysis  [71].  It  is  particularly  suited  to 
nonlinear  problems  involving  large  displacements, 
large  strains,  and  nonlinear  materials  [  6  9  j  . 
However,  it  has  a  limited  element  library;  it  is 
difficult  to  obtain  technical  support;  and  is  not 

especially  user- ori ented  [68].  A  recent  study 
[72,7oj  of  ADINA  revealed  the  following:  (1)  it 

has  well-written  users  and  theoretical  manuals,  (2) 
it  is  easy  for  the  user  to  add  newly  developed 
elements  or  materials,  (3)  it  is  efficient  in 
handling  simpler,  linear  problems,  (4)  it  is 

relatively  machi ne- independent ,  (5)  it  can  be  used 
for  large  problems,  (6)  it  is  well-suited  as  a 
research  tool,  (7)  it  has  poor  data  generation  and 
error  checking  features,  (8)  it  has  trouble  with 
complex  structures  due  to  a  limited  element 
library,  (9)  its  method  for  handling  large 
deformation  problems  is  incomplete;  (10)  it  cannot 
handle  the  situation  where  load-deformation 
response  exhibits  softening  behavior  and  then 

suddenly  becomes  stiffened,  and  (11)  it  is 
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difficult  to  determine  the  optimum  time  increment 
when  solving  nonlinear  dynamics  problems. 

When  considering  microelectronics  nonlinear 
problems  (solder  behavior,  bottoming  out,  nonlinear 
springs,  etc.),  the  AD  I NA  program  had  deficiencies 
as  noted  by  comments  (7)  through  (11)  above.  Based 
on  the  findings  of  Dr.  Chang  [68,72,73]  it  was  felt 
that  A8AQUS  was  a  more  suitable  program  for 
nonlinear  microel ectroni  c  analysis  than  ADINA. 

2.  FINITE  -  This  is  a  program  similar  to  STRUDL  and 

is,  therefore,  oriented  towards  civil  engineering 
problems.  It  solves  both  linear  and  nonlinear 

statics  problems.  Its  dynamics  capabilities  are 
still  being  developed.  It  has  a  good  element 
library  and  seems  best  equipped  at  solving 
nonlinear  problems.  However,  elements  for 
geometric  nonlinearities  are  not  yet  available. 
Its  strengths  are  in  material  nonlinear  problems 
and  in  substructuring  [70]  and  appears  to  be  a  good 
program  when  solving  these  types  of  problems.  When 
geometric  nonlinear  elements  and  a  dynamic  analysis 
capability  are  added,  it  may  become  a  good 
multi-purpose  program.  Until  then,  it  is  not 
considered  a  useful  program  for  microelectronic 
analysis  . 

3.  MARC-CDC  -  This  program  offers  the  most  advanced 
technology  for  nonlinear  static  analysis  [58,69]. 
Although  it  does  have  linear  static  and  dynamic 
capabilities  [64],  it  should  be  used  primarily  as  a 
nonlinear  analysis  program.  It  has  a  good  element 
library  for  one  and  two-dimensional  problems.  It 


50 


has  a  limited  library,  however,  of  three-dimen¬ 
sional  elements  [58,66,67].  Its  dynamic 

capabilities  are  limited  [58,64,69].  Some  users 
reported  the  program  was  not  easy  to  use  and  poorly 
documented  [58,64,65,78],  however,  improvements 
have  been  made  in  recent  versions.  Agreement  of 
results  as  calculated  by  MARC-CDC  with  experimental 
results  indicates  that  a  high  degree  of  confidence 
can  be  placed  in  a  MARC  analysis  [64],  When  the 
dynamic  capabilities  are  brought  on  a  par  with  the 
static  capabilities,  the  program  may  become  a 
powerful  tool  for  nonlinear  dynamics  analysis. 
Until  then,  MARC  should  only  be  used  for  solving 
those  nonlinear  microelectronic  problems  that 
cannot  be  handled  by  ABAQUS,  ADINA,  or  ANSYS . 

4.  NASTRAN  -  This  is  probably  the  largest  general 
purpose  program  available  [63,64].  Three  versions 
[78]  of  NASTRAN  are  available:  (1)  COSMIC,  a 

vers  i  oh  distributed  by  NASA,  (2)  a  version 
distributed  by  S pe r ry - U n i v ac ,  and  (3)  MSC,  a 
version  distributed  by  McNeal-Schwendler 
Corporation.  It  has  been  the  subject  of  more 
discussion  than  any  other  program.  Part  of  the 
reason  for  NASTRAN's  popularity  could  be  because 
use  of  this  program  has  been  reouired  on  many  NASA 
and  nuclear  power  contracts.  It  provides  a  wide 
base  of  analysis  capabilities  with  emphasis  on 
aerospace  problems.  This  program,  however, 
requires  extensive  user  experience  [58,64].  It  is 
basically  a  linear  program,  capable  of  static, 
dynamic,  buckling,  heat  transfer,  hydro-el  as ti c , 
and  aeroelastic  analysis  [61,63].  Its  documenta- 
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tion  is  extensive,  but  its  size  is  likely  to  scare 
the  potential  user  [58,64,65].  It  is  expensive  for 
solving  smaller  problems  [58],  and  is,  therefore, 
better  suited  for  larger  problems.  NASTRAN  is  not 
considered  a  user-oriented  program  [78].  Its  use 
should  be  limited  to  large  problems,  or  those 
involving  special  capabilities,  such  as  aeroelas- 
ticity.  It  does  not  do  a  good  job  of  solving 
sinusoidal  and  random  vibration  response  problems 
[61]. 

Concerning  microelectronic  analysis,  the  following 
reasons  are  given  explaining  why  NASTRAN  is  not  a 
suitable  package.  It  is  difficult  to  use  and 
requires  a  lot  of  experience  for  effective  use. 
This  makes  NASTRAN  hard  to  use  for  the 
inexperienced  analyst  or  the  experienced  analyst 
not  familiar  with  NASTRAN.  Microelectronic 
problems  are  typically  of  the  smal 1  - to-med i urn  size 
range.  NASTRAN  is  not  cost  effective  in  solving 
smaller  problems.  Its  wide-ranging  capabilities 
make  it  difficult  for  the  user  to  locate  the 
portion  of  the  program  suited  to  his  needs  (bigger 
is  not  necessarily  better  for  finite  element 
software).  Microelectronics  packages  nearly  always 
experience  sinusoidal  or  random  vibration,  and 
NASTRAN  is  not  particularly  effective  when  solving 
these  problems  . 

5.  SAP  IV  -  A  more  current  version  of  this  program  is 
SAP  VI.  However,  applicable  information  was 
obtained  only  for  the  SAP  IV  version.  This  program 
performs  linear  static  and  dynamic  analysis.  Users 
report  that  it  is  easy  to  use  [58,64,68].  It  is 
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the  onl y  program  that  is  simple  enough  that  users 
have  tried  to  modify  it  to  suit  their  needs 
[58,64].  It  is  therefore  useful  for  research-type 
problems  requiring  considerable  interface  with  user 
routines.  However,  it  has  limited  dynamic  capa¬ 
bilities  [58].  These  limited  capabilities  include 
no  Guyan  reduction,  only  one  method  of  eigenvalue 
extraction  (subspace  iteration),  only  one  value  of 
damping  can  be  used  for  transient  analysis,  and 
centrifugal  loadings  are  not  available.  SAP  IV 
does  not  appear  to  have  heat  transfer  analysis 
capabilities.  It  is  most  useful  when  applied  to 
small  or  medium -sized  linear  problems.  When 
compared  with  other  programs,  STARDYNE  was  ranked 
first  while  SAP  IV  was  ranked  third  in  order  of 
preference  when  consideration  to  solving  small-to- 
medium  range  dynamics  problems  was  given  [58].  Due 
to  its  limited  dynamics  capabilities  and  its 
comparison  with  other  programs  solving  small-to- 
medium  range  dynamics  problems,  SAP  IV  was  found  to 
be  a  less  desirable  program  than  STARDYNE  for 
microelectronic  analysis. 

The  five  programs  that  have  just  been  eliminated  from 
the  interim  program  list  are  all  very  good  general  purpose  programs 
within  the  limitations  already  discussed.  Their  use  should  only 
include  problems  that  address  these  program’s  strengths. 

4 .6  Fi nal  P  rogram  List 

Below  is  the  1st  of  the  four  programs  that  remain, 
along  with  a  brief  description  of  each  program's  p r i m a ry  functions: 

1.  A8AQUS  -  nonlinear  statics  and  dynamics  analysis 

2.  ANSYS  -  general  purpose  (linear  and  nonlinear 
statics  and  dynamics  and  heat  transfer  analysis) 
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3.  N I S  A  -  general  purpose  (linear  and  nonlinear 

statics  and  dynamics  and  heat  transfer  analysis) 

4.  STARDYNE  -  linear  statics  and  dynamics  analysis 

The  following  is  a  detailed  summary  of  each  program's 
capabilities,  strengths,  and  weaknesses. 

1.  ABAQUS  -  This  program  [74]  performs  linear  and 

nonlinear  static,  dynamic,  and  heat  transfer 
analysis.  It  is  best  suited  for  nonlinear  problems 
encountered  in  static  and  dynamic  analysis.  It  can 
solve  the  three  types  of  nonlinear  problems:  (1) 

material  nonl ineari ties  ,  (2)  geometric  nonlineari¬ 
ties,  and  (3)  la-ge  displacements  and  rotations. 

It  has  a  good  element  library  including  axisym- 

metric  elements.  It  has  a  full  range  of  material 
models  including  viscoelastic,  elastic-plastic, 
creep,  and  volumetric  swelling.  It  has  the 
capability  of  following  a  static  analysis  with  a 
dynamic  analysis  (and  vice  versa),  such  as  in  the 
dynamic  response  of  preloaded  structures.  It  can 
solve  linear  or  nonlinear  heat  transfer  problems 
and  can  provide  temperatures  for  loading  a  thermal 
stress  analysis  model. 

Dr.  Chang  [68]  had  the  following  comments  about 
ABAQUS:  (1)  it  is  the  best  program  for  nonlinear 

problems,  (2)  it  has  a  hotline  which  provides 
excellent  technical  support,  (3)  it  is  easier  to 
use  than  ANSYS,  (4)  it  has  a  better  element  library 
than  ADINA,  and  (5)  it  is  use  r- or  i  ented  and  is 
easier  to  use  than  ADINA.  Dr.  Chang  is  currently 
doing  a  study  on  ABAQUS.  His  interim  report  was 
not  yet  available  at  the  time  of  this  writing. 
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The  user  manual  appears  to  be  well-written  with 
good  example  problems.  Sligh*  difficulty  was  noted 
in  interpreting  program  commands.  Free  format 
input  makes  data  generation  simple.  It  appears 
that  ABAQUS  is  a  user-ori ented  program.  Its  use 
should  be  limited  to  nonlinear  mi c roel ec tron i cs 
behavior  such  as  solder  behavior,  bottoming  out,  or 
nonlinear  springs. 

2.  ANS  Y  S  -  This  program  [75]  is  a  large  general 
purpose  package  that  provides  linear  and  nonlinear 
static  and  dynamic  and  heat  transfer  analysis.  It 
is  used  heavily  in  the  nuclear  industry.  It 
performs  linear  and  nonlinear  elastic  analysis  of 
structures  subjected  to  static  loads  as  well  as 
transient  and  harmonic  dynamic  excitations.  It 
does  not  perform  random  vibration  response 
analysis.  The  program  considers  the  nonlinear 

effects  of  plasticity,  creep,  swelling,  and  large 
deformations.  Transient  and  steady  state  heat 
transfer  analyses  consider  conduction,  convection, 
and  radiation  effects.  Coupled  thermal -el ectric 
and  wave  motion  analysis  capabilities  are 
available.  ANSYS  also  predicts  steady  state  and 
transient  fluid  flow  in  one  -  dimens i on al  networks. 
Temperatures  obtained  from  thermal  analyses  can  be 
input  as  loadings  for  static  stress  analyses. 

It  is  considered  a  versatile  program  [58]  with 
a  wide  variety  of  engineering  applications 
possible.  Its  strength  is  its  comprehensive 
element  library.  It  appears  to  have  the  most 
complete  library  of  heat  transfer  elements 
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available.  It  has  excellent  node  and  element 
generation  capabilities. 

It  has  a  well-written  user's  manual  and  an 
excellent,  thorough  examples  manual.  It  is 
relatively  user- or i en te d  ,  and  excellent  technical 
support  is  available.  ANSYS  will  be  best  utilized 
for  microelectronics  problems  involving  heat 
transfer  solutions  as  inputs  for  thermal  stress 
c a  1 c ul a t i on s .  We  have  used  the  axisymmetric 
elements  for  microelectronics  thermal  stress 
analysis  [81]  and  found  the  program  to  be 
excel  lent. 

3.  NISA  -  This  program  [76]  is  a  large  general  purpose 
program  that  provides  linear  and  nonlinear  static 
and  dynamic  and  heat  transfer  analysis.  It  has  an 
extensive  element  library  including  isoparametric 
linear,  parbolic,  cubic,  linear  parabolic,  linear 
parabolic  cubic,  etc.  for  plane  stress  and  plain 
strain  problems.  Additional  elements  include 
axisymmetric  (with  symmetric  or  unsymmetric 
loading),  general  shells,  laminated  composite  or 
sandwich  shells,  thick  shells,  solids,  beams, 
spars,  springs,  mass  elements  and  rigid  elements. 
A  unique  program  feature  is  its  treatment  of 
composite  and  sandwich  structures.  It  has  good 
node  and  element  generation  capabilities  [58]. 

NISA  can  solve  nonlinear  problems  involving 
geometric  or  material  nonlinearity  or  both. 
Dynamic  analysis  capabilities  include  eigen- 
val  ue/ei  genvector  extraction,  transient  analysis, 
shock  spectrum  analysis,  harmonic  analysis,  and 
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random  vibration  analysis.  Heat  transfer  capa¬ 
bilities  include  steady  state  and  transient 
problems  involving  conduction,  convection,  and 
radiation.  Temperatures  obtained  from  thermal 
analyses  can  be  input  as  loadings  to  static  stress 
analyses.  A  recent  study  [79]  found  NISA  to  be  the 
best  program  for  analyzing  tactical  shelters. 

It  appears  to  have  a  well-written  user's  manual 
and  seems  to  be  a  user-oriented  program.  An 
extensive  examples  and  verification  manual  is 
available  with  several  good  example  problems 
illustrated.  Good  technical  support  is  available. 
NISA  appears  well-suited  for  solving  the  myriad  of 
m  i  c r oe 1 ec t r on i c s  problems.  Its  only  apparent 
limitations  are  in  solving  highly  nonlinear 
problems  . 

4.  STARO YNE  -  This  program  [77]  performs  linear  static 
and  dynamic  analysis.  It  has  limited  nonlinear 
static  and  dynamic  capabilities.  Nonlinear  problem 
applications  are  limited  to  geometric  nonlineari¬ 
ties  such  as  gaps  and  bottoming  out  of  adjacent 
structures.  It  has  a  strong  ability,  however  to 
^olve  dynamics  problems.  The  program  has  a  simple, 
yet  good,  element  library  with  enough  elem*nts  to 
describe  most  structures  under  static  or  dynamic 
loading.  The  same  model  can  be  used  for  static  and 
dynamic  analysis,  saving  time,  eliminating  errors, 
and  saving  computer  costs. 

Dynamic  analysis  capabilities  include 
eigenvalue/eigenvector  extraction,  transient 
analysis,  shock  analysis,  harmonic  analysis,  random 
vibration  analysis,  and  shock  spectrum  analysis. 
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It  is  considered  to  have  the  best  ability  for 
solving  sinusoidal  vibration  response  and  random 
vibraton  response  problems  [61].  It  has  a  very 
efficient  method  for  extracting  natural  frequencies 
and  mode  shapes  [69].  It  was  rated  the  best 
general  purpose  program  for  solving  small  to 
medium-sized  dynamics  problems  [58].  STARDYNE  does 
not  have  heat  transfer  capabilities.  However,  it 
is  capable  of  performing  thermal  stress  analysis. 

The  program's  strengths  are  its  simplicity,  its 
efficient  solution  of  dynamics  problems,  and  the 
fact  that  it  is  extremely  user-oriented.  The 
user's  manual  is  well-written  and  provides  a 
simplified  tool  for  learning  the  program.  There  i s 
an  excellent  learner's  manual,  examples  manual,  and 
theoretical  manual  available.  STARDYNE  is  one  of 
few  programs  that  can  be  learned  by  simply  reading 
the  user ‘ s  manual  . 

Since  STARDYNE  is  particularly  adept  at  solving 
small-to-medium  range  statics  and  dynamics 
problems,  it  is  probably  the  best  program  for  per¬ 
forming  the  majority  of  microelectronics  finite 
element  analyses.  The  writers  have  used  STARDYNE 
for  solving  general  electronics  packaging  and 
microelectronics  packaging  problems  [81],  and  have 
found  the  program  to  be  excellent. 


4  .7  Comparison  of  Final  Program  Capabilities 

To  compare  program  capabilities,  tables  of  parameters 
were  made  in  the  following  areas: 

1.  Types  of  Analysis 

2.  Material  Properties 


3.  Element  Library 


4.  Types  of  Loading 

5.  Data  Input  and  Output 

6  .  Ooc  umentati on 

These  program  capabilities  are  summarized  in  Tables  9  through  14 
which  follow. 

4  .8  General  Summary  of  Analysis  Methods  (Phase  2) 

1.  A  review  of  21  general  purpose  finite  element 
programs  resulted  in  an  interim  list  of  nine 
programs  .  Five  programs  -  ADINA,  FINITE,  MARC-CDC, 
NASTRAN,  and  SAP  IV  -  were  then  eliminated  from 
further  consideration.  These  five  programs  were 
rated  "very  good";  however,  their  use  should  only 
include  those  unique  microelectronic  problems  that 
address  specific  program  strengths. 

2.  The  four  recommended  programs  are  ABAQUS ,  ANSYS, 
NISA,  and  STAROYNE.  These  recommended  programs 
will  best  be  utilized  when  performing  the  following 
types  of  microelectronic  analyses: 

a.  ABAQUS  -  nonlinear  statics  and  dynamics. 
Nonlinear  effects  may  include  material 
n on  1  i ne a r i t i e s  such  as  a  solder  creep, 
geometric  nonlinearities,  such  as  a  substrate 
bottoming  out  on  a  module,  and  large 
displacement  nonlinearities,  such  as  a  hybrid 
lead  bending  under  load. 

b.  ANSYS  -  linear  and  nonlinear  statics  and 
dynamics  and  heat  transfer.  ANSYS  is  best 
suited  for  heat  transfer  analysis  where 
generated  temperatures  can  be  used  as  Inputs 
for  a  thermal  stress  analysis  such  as  required 
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Table  9.  Types  of  Analysis 


PROGRAM 

CAPABILITY 

ABAQUS 

ANSYS 

NISA 

STARDYNE 

Linear  Statics 

X 

X 

X 

X 

Thermal  Stress 

X 

X 

X 

X 

Nonl i near  Stati cs 

X 

X 

X 

X 

Geometri  c 

X 

X 

X 

X 

Materi al 

X 

X 

X 

Large  Deflection 

X 

X 

X 

Linear  Dynamics 

X 

X 

X 

X 

Modal  Extraction 

X 

X 

X 

X 

T  rans i ent 

X 

X 

X 

Harmoni c 

X 

X 

X 

Random  Yibration 

X 

X 

Shock  Spectra 

X 

X 

Nonl i near  Dynamics 

X 

X 

X 

X 

Transient 

X 

X 

X 

X 

Geometri c 

X 

X 

X 

X 

Materi al 

X 

X 

X 

Large  Deflection 

X 

X 

X 

Buckl i ng 

X 

X 

Linear  Heat  Transfer 

X 

X 

X 

Nonlinear  Heat  Transfer 

X 

X 

X 

FI  uid  Flow 

X 
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Table  10.  Material  Properties 


CAPABILITY 

ABAQUS 

PROGRAM 

ANSYS  NISA 

STARDYNE 

Linear  Elastic 

X 

X 

X 

X 

Isotropi c 

X 

X 

X 

X 

Ani sotropi c 

X 

X 

X 

X 

Orthotropi c 

X 

X 

X 

Temperature-Dependent 

X 

X 

P 1 asti c 

X 

X 

X 

V i scoel asti c 

X 

X 

Creep 

X 

X 

Swell i ng 

X 

X 

Sandwl ch 

X 

X 

X 

Composi te 

X 
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Table  11.  Element  Library 


PROGRAM 

CAPABILITY 

ABAQUS 

ANSYS 

NISA 

STARDYNE 

Bar  Element 

X 

X 

X 

X 

Beam  Element 

X 

X 

X 

X 

Membrane  Plate  Element 

X 

X 

X 

X 

Bending  Plate  Element 

X 

X 

X 

X 

Thin  Shell  Element 

X 

X 

X 

X 

Thick  Shell  Element 

X 

X 

Isoparametric  8  Node  Solid 

X 

X 

X 

X 

Isoparametri c  20  Node  Solid 

X 

X 

X 

Axisymmetric  Shell  Element 

X 

X 

Axlsymmetric  Solid  Element 

X 

X 

X 

Pipe  Element 

X 

X 

Gap  Element 

X 

X 

X 

X 

Friction  Element 

X 

X 

Spring,  Mass,  Damper  Elements 

X 

X 

X 
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Table  12. 

CAPABILITY 

Point  Loads 

Line  Loads 

Surface  Loads 
Volume  Loads 
Thermal  Loading 
Centrifugal  Loads 
Axisymmetric  Loads 
Prescribed  Displacements 
Elastic  Foundation 
Time  Dependent 
Deformation  Dependent 
Contact 
Fri cti on 


Types  of  Loading 

PROGRAM 

ABAQUS  AMSYS  NISA 
XXX 

X  X 

XXX 

X  X 

XXX 

XXX 

XXX 

XXX 

X 

XXX 

X 

XXX 

XXX 


STARDYNE 

X 

X 

X 

X 

X 

X 

X 

X 
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Table  13.  Data  Input  and  Output 


CAPABILITY 


PROGRAM 

ABAQUS  ANSYS  NISA 


Input 

Node  Generation 
Element  Generation 
Load  Generation 
Interactive  Graphics 
Restart  Capability 
Free  Format  Input 
Plot  Routines 
Data  Input  Check 


X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 


Cutput 

Numerical  Results  X 
Graphical  Results  X 
Plotting 

Oeformed, Undeformed  X 
Temperature,  Stress  Contours  X 
Oynamic  Response  X 


X  X 

X  X 


X  X 

X  X 

X  X 


STAROYNE 


X 

X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 


Table  14.  Documentation 


PROGRAM 


CAPABILITY 

ABAQUS 

ANSYS 

NISA 

STARDYNE 

User’s  Manual 

X 

X 

X 

X 

Programmer's  Manual 

X 

Theoretical  Manual 

X 

X 

Sample  Problems  Manual 

X 

X 

X 

X 

Verification  Manual 

X 

X 

Technical  Support  Available 

X 

X 

X 

X 
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when  analyzing  bonded  MIC  or  stripline 
substrates  . 

c.  NISA  -  linear  and  nonlinear  statics  and 

dynamics  and  heat  transfer.  NISA  is  best 
suited  for  analyzing  composite  or  sandwich 
structures  such  as  a  circuit  board  with  layered 
copper  ci  rcuit  runs  . 

d.  STARDYNE  -  linear  statics  and  dynamics. 

STARDYNE  is  best  suited  for  analyzing  dynamics 
problems  such  as  circuit  board  response  and 
hybrid  lead  dynamic  stresses  due  to  random  or 
sinusoidal  vibration. 

3.  STARDYNE  is  the  best  program  for  analyzing  the 

majority  of  microelectronic  problems. 

4.  The  area  of  nonlinear  analysis  -  particularly 

solder  material  effects  -  is  a  promising  area  for 
future  work  and  should  be  investigated. 

5  .0  SAMPLE  PROBLEM  ANALYSIS  (PHASE  2) 

5  .1  Introducti on 

This  section  of  the  report  will  discuss  each  of  the 
four  sample  problems.  The  finite  element  models  used  for  each 
problem  will  be  explained,  and  the  results  from  the  different 
finite  element  programs  will  be  compared  with  each  other,  and 
also  with  the  theoretically  correct  results. 

The  sample  problems  were  chosen  using  the  following 
criteria:  (1)  they  must  represent  typical  analysis  problems 

associated  with  microelectronic  packaging,  and  (2)  they  should 
be  simple  enough  that  a  theoretically  "correct”  answer  is 
available  for  comparison.  Consistent  with  these  criteria,  four 
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problems  were  chosen  from  the  general  categories  of  linear 
statics,  linear  dynamics,  linear  steady-state  heat  transfer,  and 
nonlinear  statics. 

Theoretically  correct  answers  were  only  available 
for  the  problems  of  the  first  three  categories  mentioned  above. 
The  analysis  approach  was  to  determine  program  accuracy  by 
comparing  computed  and  theoretical  answers,  and  to  note  any 
program  idiosyncrasies.  Noting  the  simplicity  of  these 
problems,  the  approach  was  not  to  determine  "blanket" 
conclusions  about  a  specific  program's  ability  to  solve  a 
specific  class  of  problems.  Instead,  it  was  desirable  to  note 
any  situation  where  a  program's  solution  was  radically  different 
from  the  other' s  . 

The  NISA  program  is  capable  of  solving  all  the  above 
types  of  problems.  However,  it  was  not  available  on  either  the 
Control  Data  Cybernet  or  the  Harris  computer  networks. 
Therefore,  NISA  was  not  evaluated  using  sample  problems.  Of  the 
remaining  three  programs,  all  can  solve  each  of  the  above 
problems  (with  the  exception  of  STARDYNE,  which  cannot  solve 
heat  transfer  problems).  Below  is  a  description  of  the  various 
sample  problems  . 

5 .2  Sample  Problem  Descriptions 

Each  of  the  sample  problems  is  described  in  the 
following  paragraphs;  detailed  discussions  of  these  problems 
will  follow  in  Sections  5.3  through  5.6. 

5.2.1  Linear  Statics  Problem 

Hybrid  designs  are  frecuently  assessed  by  applying  a 
uniform  external  pressure  to  the  package.  Therefore,  it  follows 
that  the  design  engineer  may  need  to  determine  the  stresses  and 
deflections  In  the  hybrid  lid  due  to  the  pressure  loading. 
Theoretically,  this  is  the  classic  problem  of  a  plate  bending 
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under  static  loading.  As  a  first  (and  possibly  conservative) 
approximation,  the  engineer  may  assume  the  lid  has  four  simply 
s upported  edges  . 

To  compare  the  finite  element  solution  to  the 
theoretically  correct  answer,  only  the  deflections  will  oe 
considered.  The  finite  element  model  was  purposely  made  with  a 
relatively  coarse  grid,  and  the  stresses  obtained  with  this 
model  will  reflect  this  coarseness.  This  problem  is  physically 
described  in  Figure  12. 

5.2.2  Linear  Dynamics  Problem 

Many  mi c roel ec tron i c  packages  are  mounted  to  printed 
circuit  boards.  Since  most  electronics  eouipment  is  evaluated 
by  dynamic  testing,  it  is  desirable  to  determine  the  dynamic 
characteri s ti cs  of  a  particular  package.  If  the  package  is 
mounted  to  a  circuit  board,  the  first  step  in  this  process  is  to 
determine  the  board's  natural  freauencies  and  mode  shapes. 

For  this  problem,  the  first  four  natural  frequencies 
of  the  circuit  board  will  be  calculated  and  compared  to  the 
theoretically  correct  solutions  (note  that  for  most  situations, 
it  is  the  first  few  1  ow-  f  re  fluency  modes  that  represent  the 
greatest  damage  potential  to  the  system).  In  addition  to  the 
natural  freauencies,  the  correspondi ng  jnode  shapes  will  also  be 
described . 

This  problem  is  shown  pictorially  in  Figure  13;  it 
represents  a  0.50  lb.,  G-10  epoxy /fi berg! ass  printed  circuit 
board . 

5 .2 .3  Heat  Transfer  Problem 

Pin  fins  are  commonly  used  in  electronics  applica¬ 
tions  to  remove  heat  from  temperature-sensitive  areas.  It  Is 
commonly  desired  to  determine  the  effectiveness  of  a  particular 
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Material :  KQVAR  , 

Young's  Modulus  =  20x10°  psi. 
Poisson’s  Ratio  =0.3 


Loading:  Uniform  exterior  pressure  of 

30  psi. 


Support:  All  four  edges  simply 
supported 


Material :  G-10  EPOXY  FIBERGLASS  , 

Young's  Modulus  =  2.5  x  10°  psi 
Poisson's  Ratio  =  0.12 
Total  Weight  =  0.50  lb. 


Support:  All  four  edges  simply  supported 


Fiaure  13.  Linear  Dynamics  Problem 
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pin  fin  design;  i.e.,  it  is  desired  to  know  the  heat  flux 

allowed  by  a  particular  fin  design. 

I  he o re t i c a  1 1 y  ,  this  problem  is  that  of  a  one¬ 
dimensional  bar  with  both  conductive  and  convective  heat 
transfer  capabi 1 i ti es  .  We  will  calculate  both  the  heat  flux  and 
the  temperature  distribution  of  the  fin,  and  the  solution  will 
be  compared  to  the  theoretically  correct  results.  The  particu¬ 
lar  fin  used  for  this  problem  is  described  in  Figure  14. 

5  .2  .4  Nonlinear  Statics  Problem 

This  problem  is  concerned  with  a  typical  design 
where  a  substrate  is  mounted  in  a  module.  The  substrate  is 

attacned  to  the  module  with  four  screws  (one  per  substrate 

corner),  and  it  is  mounted  to  small  bosses  in  the  module  such 

that  the  substrate  is  supported  off  the  module  floor. 

This  problem  is  considered  to  be  nonlinear  because 
it  exhibits  different  responses  depending  upon  the  direction  of 
loading.  For  instance,  in  the  event  of  an  upward  load  (produced 
by  a  downward  acceleration  of  the  module),  the  substrate  is 

restrained  by  the  four  mounting  screws  only.  However,  in  the 
case  of  a  downward  load,  it  is  possible  that  the  center  of  the 

substrate  may  deflect  enough  to  contact  the  module  floor.  Thus, 

should  the  s ubstrate- to-mod ul e  gap  be  closed,  the  boundary 
conditions  must  consider  the  effect  of  the  module  floor  in 

addition  to  the  four  mounting  screws. 

From  a  static  standpoint,  the  most  likely 
environment  to  produce  the  above-discussed  responses  is  a 
constant  acceleration.  The  20,000g  acceleration  chosen  for  this 
problem  is  consistent  with  that  found  in  MIL-STD-883B,  Method 
2G01.2,  Test  Condition  0.  The  substrate  is  physically  described 
in  Figure  15  . 
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MATERIAL  6061 -T6  Aluminum 
Conductivity  =  96  Btu/hr  ft°F 


AIR  FLOW 

Flow  Speed  =  200  fpm 

Flow  Temperature  =  131°F  ? 

Heat  Transfer  Coefficient  =  12.7  Btu/hr-ft' •  F  (For  both 
the  circumferential  surface  area  and  the  end  surface  area) 

IMPOSED  BOUNDARY  CONDITION 

i 

Constant  base  temperature  =  176°F 


Figure  14.  Heat  Transfer  Problem 
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0.020" 


SUBSTRATE 


Support;  Fixed  support  at  all 
four  corners 


Displacement  Restraints: 

1)  Upward  load  -  none 

2)  Downward  load  -  vertical 
translational  restraint  if 
gap  is  closed 


Figure  15.  Nonlinear  Statics  Problem 
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5  .3  Sample  Problem  #1  -  Linear  Statics 

The  linear  statics  sample  problem  involved  the 
determination  of  the  deflection  for  a  hybrid  lid  due  to  a  30  psi 
pressure  loading.  This  problem  was  solved  using  the  following 
finite  element  programs:  (1)  STARDYNE,  (2)  ANSYS,  and  (3) 

ABAQUS  .  Due  to  the  symmetry  of  the  chosen  hybrid  lid,  a 

"cuarter"  model  was  used  without  any  loss  of  accuracy;  this 
model  is  shown  in  Figure  16. 

As  shown  in  Figure  16,  the  cuarter  model  for  the 

hybrid  lid  consists  of  sixteen  nodes  and  nine  euadri  1  ateral 
plate  elements.  The  boundary  conditions  for  the  four  edges  of 
the  model  were  grouped  into  two  classifications:  Those 

necessary  to  represent  the  simply-supported  edges,  and  those 
necessary  to  provide  continuity  along  the  two  lines  of  symmetry. 
Edge  1-2-3-4  and  edge  1-5-9-13  were  allowed  to  rotate  about  the 
X-axis  and  the  Y-axis  respectively;  all  other  degrees-of -freedom 
(DOF)  were  assumed  to  be  fully  restrained  for  these  nodes. 
Also,  due  to  the  hybrid  lid  symmetry,  edge  4-8-12-16  was  not 
allowed  to  translate  in  the  X-direction,  nor  was  it  allowed  to 
rotate  about  the  Y-axis  (the  slope  of  edge  4-8-12-16  in  the 
X-direction  is  constrained  to  be  zero  by  symmetry) .  Similarly, 
edge  13-14-15-16  was  restrained  against  Y-direction  translations 
and  X-axis  rotations.  These  edge  restraints  resulted  in  node  16 
being  free  only  to  translate  in  the  Z-di  recti  on  since  it  is 
located  at  the  intersection  of  the  two  lines  of  symmetry. 

In  addition  to  the  boundary  conditions  listed  above, 
one  additional  restraint  was  required  to  prevent  any  nodal 
rotations  about  the  axis  perpendicular  to  the  plane  of  the  lid. 
STARDYNE  utilizes  a  plate  element  which  has  5  DOF  (degrees  of 
freedom)  per  node;  there  is  no  stiffness  associated  with 
rotations  about  an  axis  pe rpend i c ul  a r  to  the  plate.  Therefore, 
when  using  STARDYNE,  it  is  necessary  to  restrain  this  rotational 
DOF  to  prevent  singularity  problems  from  arising  during  the 
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/bri d  Lid  Outline  (0.92"  x  0.92") 


Axes  of  Symmetry 


"Quarter"  Finite  Element  Model 
(0.46"  x  0.46") 

-Element  #'s  are  circled  - 


Figure  16.  Hybrid  Lid  Finite  Element  Model 


matrix  solution.  Similarly,  the  plate  elements  of  ABAQUS  have 
only  5  DOF  per  node,  and  the  rotational  DOF  about  the  normal 
axis  must  be  restrained  (ABAQUS  will  automatically  restrain  this 
DOF  if  necessary).  On  the  other  hand,  ANSYS  provides  an  option 
for  including  a  small  rotational  stiffness  for  rotations  normal 
to  the  plane  of  the  plate.  However,  as  discussed  in  the  ANSYS 
manual,  this  rotational  stiffness  is  usually  the  least  important 
of  all  stiffness  components  for  the  element  and  may  be  safely 
neglected.  Therefore,  to  keep  the  finite  element  model  of  the 
lid  consistent  for  each  of  the  three  programs,  rotations  about 
an  axis  perpendicular  to  the  plane  of  the  lid  were  restrained 
for  all  sol uti ons . 

The  remaining  physical  and  material  properties 
reouired  for  this  analysis  are  shown  below  in  Table  15. 


Table  15  . 

Physical  Dimensions 

Plate  Type 

Material 

Young  '  s  Mod  ul  us 

Poisson's  Ratio 

Loading 


Physical  and  Material  Properties 

0  .92"  x  0  .92"  x  0  .015" 
Bending  and  Membrane 
Kov  ar 

20  x  106  psi 
0  .3 

30  psi  Pressure 


As  shown  in  Table  15,  the  plates  used  for  this 
finite  element  model  had  both  bending  stiffness  and  membrane 
stiffness.  In  each  case,  the  ouadri 1 ateral  plate  element  which 
was  used  had  four  nodes  per  element.  The  element  used  for  the 
STARDYNE  analysis  was  the  "QUAD”  element,  the  element  used  for 
the  ANSYS  analysis  was  the  “STIF43"  element,  and  the  element 
used  for  the  ABAQUS  analysis  was  the  "S4R"  element.  Although 
both  ANSYS  and  ABAQUS  offer  8  -  node- pe r-el ement  quadrilateral 
plate  elements,  the  more  simple  4-node- per-el ement  plate 


elements  were  used  in  order  to  provide  a  valid  comparison  with 
STARDYNE's  4-node  ouadri 1 ateral  element. 

A  listing  of  the  required  input  for  each  of  the 
three  finite  element  programs  is  shown  in  Figures  17  through  19. 
Both  STARDYNE  and  ANSYS  require  a  fixed  format  input,  while 
ABAQUS  has  a  free  format  option.  For  a  complete  understanding 
of  these  listings,  refer  to  the  appropriate  user  manual  [82], 
[83],  [84]  . 

The  nodes  and  elements  for  this  model  were  generated 
more  rapidly  by  AK'SYS  than  by  STARDYNE  or  ABAQUS.  The  "second 
level  generation"  capability  of  ANSYS  was  used  for  generating 
the  model's  elements  and  nodes  (ANSYS  also  provides  "third  level 
generation"  capability  for  three-dimensional  gridworks).  For 
example,  nodes  1  through  4  were  defined  by  specifying  the 
locations  of  node  1  and  node  4,  and  then  telling  the  program  to 
generate  a  total  of  four  eoually  spaced  nodes  along  the  defined 
line.  Next,  the  remaining  twelve  nodes  were  defined  by  telling 
the  program  to  repeat  the  nodal  set  of  nodes  1-4  an  additional 
three  times,  incrementing  the  Y-coordinate  by  0.1533"  each  time. 
Similarly,  elements  1-3  were  generated  by  copying  element  1 
twice,  and  then  repeating  the  element  set  1-3  twice. 

STARDYNE  has  less  node  and  element  generation  capa¬ 
bility  than  ANSYS,  although  for  this  problem  the  difference  was 
not  great.  STARDYNE  allowed  second  level  generation  of  the 
nodes  using  its  MODEGRD  option,  but  this  program  is  limited  to 
first  level  generation  of  elements  (three  separate  commands  were 
reculred  to  generate  elements  1-3,  4-6,  and  7-9). 

ABAQUS  has  the  least  amount  of  generation  capability 
among  the  three  programs  of  interest;  ABAQUS  provides  first 
level  generation  of  nodes  and  elements  only.  However,  this 
should  not  be  considered  a  significant  limitation  of  the 
program.  For  many  structures  analyzed  with  finite  elements,  the 
geometry  is  too  complex  to  lend  itself  to  a  large  amount  of  node 


or  element  generation.  It  is  only  for  large,  repetitive, 
"gri dwork-type"  models  that  generation  can  save  the  analyst  a 
significant  amount  of  time  and  effort. 

5  .3  .1  Results  for  Sample  Problem  #1 

Table  16  shows  the  computed  transverse  deflections 
for  each  of  the  nine  nodes  with  an  unrestrained  z-transl ati onal 
DOF.  The  deflection  at  node  16  is  the  maximum  deflection  at  the- 
center  .f  the  hybrid  lid. 


Table  16. 

Deflections  of 

a  Hy  b  r  i  d  Lid 

(inches) 

Node 

STARDYNE 

ANSYS 

ABAQUS 

Theoreti cal 

6 

.003779 

.003960 

.003772 

- 

7 

.006337 

.006620 

.006383 

- 

8 

.007224 

.007540 

.007289 

- 

10 

.006337 

.006620 

.006383 

- 

11 

.010679 

.011149 

.010857 

- 

12 

.012194 

.012726 

.012426 

- 

14 

.007224 

.007540 

.007289 

- 

15 

.012194 

.012726 

.012426 

- 

16 

.013931 

.014538 

.014239 

.014137 

Dev* 

-1  .4577, 

+  2  .837% 

+  0  .722% 

*Percent 

Jevi ati on 

of  the  computed 

deflection  for  node  16  from 

the  theoretical  maximum  deflection  at  the  center  of  the  lid 

As  shown  in  the  above  table,  the  computed  values  of 
the  maximum  lid  deflection  are  all  within  3%  of  the  theoretical 
solution.  The  theoretical  solution  was  that  found  in  the  5th 
edition  of  Formulas  for  Stress  and  Strain  by  R.  J.  Roark  [85]. 
The  solution  is  detailed  on  page  83. 
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Figure  17.  STARDYNE  Input  Listing  for  Problem  #1 
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Figure  18,  ANSYS  Input  Listing  for  Problem  #1 
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ABACUS  INPUT  ECHO 


PAGE  1 


5 


10  IS  20  25  30 


35  AO  55  50 


55  60  65  70  75 


5 


10 


15 


20 


25 


30 


35 


50 


55 


50 


♦  HE A3  IMG 


RAOC 
♦  NODE 

SAMPLE  PROBLEM  1  .. 

1 

0.0  0.0 

5 

.56  0.0 

13 

0.0  .56 

16 

.56  .56 

♦  NGEN 

1.13.5 

5.16.5 

1.5 
5,8 
9,12 

13,16 

♦ELEMENT,  TYPE-S5R 

1.1. 2. 6. 5 

♦ELGEN,  ELSET-ALL 

1.3. 5. 3 

1.3 

5.3 

7.3 

♦SHELL  SECTION,  ELSET-ALL 
.015 

♦MATERIAL,  ELSET-ALL 
♦ELASTIC,  TYPE-ISO 

20. E6  0.3 

♦BOUNDARY 

1  ,1,6 

2.2.3 

2.5.6 

3.2.3 

3.5.6 

5.1.3 

5.5.5 

5.1 

5.3.5 

5.6 

3.1 

7.3.5 

3.6 

13.1.5 

13.6 

8.1 

3.5.5 

12,1 

12.5.6 
15,2 

15.5 

15.6 
15,2 


HYBRID  LIO  UNDER  UNIFORM  PRESSURE 

0.0 

0.0 

0.0 

0.0 


( ABAOUSI 


Figure  19.  ABAQUS  Input  Listing  for  Problem  #1 
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2 


where: 


y  =  -  a  q  b4 

Jmax  ^t3 

«  =  0  .0444  for  a/b  =1.0 

0  =  30  psi  pressure 

b  =  0.92" 

E  =  20  x  106  psi 

t  =  0.015" 

The  results  shown  in  Table  16  indicate  that  ABAQUS 
most  closely  approached  the  value  of  the  theoreti cal  1 y- cor rec t 
maximum  deflection,  followed  closely  by  STARDYNE  and  ANSYS. 

A  second  comparison  which  could  be  made  between  each 
of  the  three  numerical  solutions  concerned  the  deflection  shapes 
of  the  deformed  lid.  The  deflection  shapes  presented  in  Table 
17  are  based  on  the  maximum  deflection  for  each  solution  being 
normalized  to  a  value  of  1.0. 

As  shown  in  Table  17,  the  deflection  shapes  are  all 
very  similar.  STARDYNE  a.id  ANSYS  have  the  least  deviation 
between  their  deflection  shapes;  the  largest  deviation  between 
their  normalized  values  is  0.516%  at  node  6.  The  deviations 
between  ABAC'JS  and  both  STARDYNE  and  ANSYS  are  only  slightly 
1 arger  . 

Table  17.  Deflection  Shape  of  a  Hybrid  Lid 


Node 

STARDYNE 

ANSYS 

ABAQ'  S 

6 

0  .2713 

0 .2724 

0  .2649 

7 

0  .4549 

0 .4554 

0  .4483 

8 

0  .5186 

0 .5186 

0  .5119 

10 

0  .4549 

0 .4554 

0  .4483 

11 

0  .7666 

0 .7669 

0  .7625 

12 

0  .8753 

0 .8754 

0  .8727 
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Table  17.  Deflection  shape  of  a  Hybrid  Lid  -  Continued 


14 

0 .5186 

0 .5186 

0 .5119 

15 

0 .8753 

0 .8754 

0 .8727 

16 

1 .0000 

1 .0000 

1 .0000 

Finally,  a  comparison  was  made  between  the 
calculated  values  for  the  maximum  stress  in  the  plate.  Although 
the  model  was  relatively  coarse,  it  was  felt  that  a  comparison 
of  the  stress  values  between  the  three  programs  would  be  valid, 
regardless  of  how  closely  these  values  approached  the 

theoretically  correct  stress.  The  maximum  bending  stress  which 
occurred  at  the  center  of  the  lid  (plate  #9)  is  summarized  below 
for  each  finite  element  program. 

Table  18.  Maximum  Hybrid  Lid  Stress 
Solution  Stress  (psi)  Deviation  (%) 

STARDYNE  29,980  -7.57% 

ANSYS  30,595  -5.67% 

ABAQUS  30,780  -5.10% 

Theoretical  32,434 

The  theoretical  solution  shown  in  Table  18  was  once 
again  taken  from  Roark's  Formulas  for  Stress  and  Strain  [85]. 
The  ecuation  is: 

0  c  b^ 

max  =  t2 

where  0  -  0  .2874  for  a/b  =  1.0 

c  =  30  psi  pressure 

b  =  0.92" 

t  =  0.015" 
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As  shown  in  Table  13,  all  of  the  computed  maximum 
stresses  were  within  8"  of  the  theoretically  correct  value.  It 
should  be  noted  that  the  maximum  stress  occurred  at  the  center 
of  plate  *9  in  the  outer  fibers  of  the  plate  surface.  However, 
the  theoretical  solution  shews  the  maximum  stress  at  the  center 
of  the  hybrid  lid  (node  16).  A  finer  element  grid  would  have 
allowed  the  finite  element  programs  to  calculate  stresses 
similar  to  the  theoretical  answer;  however,  for  most  structural 
analyses,  the  numbers  shown  in  Table  18  would  be  completely 
ade  cuate  . 

This  completes  the  discussion  of  the  first  sample 
problem.  The  next  section  of  this  report  will  deal  with  Sample 
Problem  #2:  the  calculation  of  natural  frequencies  and  mode 
shapes  for  a  printed  circuit  board. 

5  .4  Sample  Problem  #2  -  Linear  Dynamics 

This  problem  involved  the  determination  of  the  first 
four  natural  frequencies  (and  the  corresponding  mode  shapes)  for 
a  printed  circuit  board.  The  problem  was  solved  using  the 
following  finite  el ement  programs:  (1)  STAR  DYNE,  (2)  ANSYS,  and 
(3)  ABAQUS  .  The  finite  element  model  of  the  circuit  board  is 
shown  in  Fig ure  20  . 

It  should  be  noted  that,  although  the  finite  element 
model  was  symmetric  about  a  pair  of  axes  parallel  to  the  X-Y 
axes  passing  through  node  13,  a  "cuarter"  model  was  not  used  for 
this  solution.  The  reason  for  not  using  a  cuarter  model  was 
that  we  desired  to  calculate  the  first  four  natural  frequencies 
of  the  system;  the  quarter  model  would  only  allow  us  to 
accurately  determine  the  natural  frequencies  for  modes  with 
shapes  symmetric  about  the  central  axes  (in  this  case,  only  the 
first  natural  frequency  would  have  been  correct). 

The  boundary  conditions  for  this  model  were  similar 
to  those  chosen  for  the  first  sample  problem;  all  four  edges 
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Figure  20.  Printed  Circuit  Board  Finite  Element 
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were  assumed  to  be  simply  supported,  or  hinged,  along  their 
length.  Also,  all  nodal  rotations  perpendicular  to  the  PC  board 
were  restrained  for  the  reasons  discussed  in  the  description  of 
the  first  sample  problem.  Other  physical  and  material 
parameters  are  summarized  in  Table  19. 


Table  19.  Physical  and  Material  Properties 


Physical  Dimensions 

Plate  Type 

Material 

Yo  ung  '  s  Modul  us 

Poisson's  Ratio 

Weight 


5  .0"  x  5  .0"  x  0  .060" 
Bending  and  Membrane 
G-10 

2.5  x  106  psi 
0.12 
0.5  lb 


As  shown  in  Table  19,  the  plates  used  for  this 
analysis  had  both  bending  stiffness  and  membrane  stiffness.  For 
each  solution,  4-node  cuadri 1 ateral  plates  were  used  to  model 
the  PC  board  ("QUAD"  elements  for  STARDYNE,  "STIF43"  elements 
for  ANSYS,  and  " S 4 R "  elements  for  ABAQUS). 

Also  shown  in  Table  19  is  the  weight  of  the  PC 
board.  This  weight  represents  the  sum  of  the  G-10  weight  and 
the  weights  of  the  individual  components  mounted  to  the  board. 

For  this  analysis,  it  was  assumed  that  the 

components  were  fairly  evenly  distributed  about  the  board,  thus 

allowing  us  to  include  their  weight  Into  the  calculation  for  the 
overall  PC  board  density.  It  was  also  assumed  that  the 

components  did  not  affect  the  overall  stiffness  of  the  board. 
For  PC  boards  with  large,  heavy  components,  neither  of  these 
assumptions  would  be  realistic,  and  a  more  detailed  finite 

element  model  would  be  recuired.  Also,  for  multilayer  circuit 
boards  with  several  layers  of  copper,  it  would  be  necessary  to 
calculate  an  "equivalent"  Young's  Modulus  for  the  copper/G-10 
composi te  . 
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A  listing  of  the  required  input  for  each  of  the 
three  finite  element  programs  is  shown  in  Figures  21-23.  Note 
that  both  STARDYNE  and  ANSYS  require  a  fixed  format  input,  while 
ABAQUS  allows  a  free  format  input.  Also,  note  that  STARDYNE 
requests  a  weight  density,  while  ANSYS  and  ABAQUS  both  request  a 
mass  density. 


5.4.1  Results  for  Sample  Problem  #2 

The  following  Table  20  shows  the  first  four  natural 
frequencies  as  computed  by  each  program.  Also,  the 
theoretically  correct  natural  frequencies  are  included  for 
comparison  purposes.  These  values  were  calculated  using  the 
equations  in  the  reference  Formulas  for  Natural  Frequency  and 
Mode  Shape  [86];  these  equations  are  summarized  below: 


2 

*  ij 
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For  mode  l,  i  =  1,  j  =  1 

For  mode  2,  i  =  2,  j  =  1  Defined  by  Blevins, 

For  mode  3,  i  =  1,  j  =  2  Table  11-4  [86] 

For  mode  4,  i  *  2,  j  *  2 


whe »  e : 


a 

b 

E 

h 


length  of  plate 
width  of  plate 
Yo ung ’ s  Modul us 
plate  thickness 
mass /unit  area 
Poisson's  Ratio 


a 

5  .0" 

a 

5  .0" 

a 

2.5  x  106 

psi 

a 

0  .060" 

a 

5.17081  x 

1 

O 
*— « 

a 

0 .12 

lbf-sec2/in3 
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Figure  21.  STARDYNE  Input  Listing  for  Problem  #2 
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RAOC  SAMPLE  PROBLEM  2  ...  PC  BOARD  VIBRATION  IANSYSI 
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Figure  22.  ANSYS  Input  Listing  for  Problem  #2 
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12  ROT 2 
IT  ROT Z 


a  a  a  o  u  S 


INPUT 


ECHO 


P_A  J,_  E _ 1 


5  10  15  20  2b  30  3b  AO  A5  bO  5b  60  65 


♦HEADING 

R ADC  SAMPLE  PROBLEM  2  ...  PC  BOARD  VIBRATION 
♦NDDE,  NSET-CORN 
1*0. 0.0 
5. 5, 0,0 
21,00,0  _ 

25,5,5,0 

♦NGEN,  NSET-EDGE 

1.5.1  _ _ _ 

21.25.1 

♦NGEN,  NSET-SIDE 

1.21.5  _ 

5.25.5 

♦NGEN,  NSET-MIO 

6,10,  1 

11.15.1 

16.20.1 

♦ELEMENT,  TTPE-SAR 
1,1,2 ,7,6 
♦ELGEN,  EL  SET-PC8 
1  «  A  ,  1 , 1 
1  ,  A  ,  5  ,  A* 

2,  A, 5, A 

3,  A, 5,  A 
A  ,  A  ,  5  ,  A 

♦SHELL  SECTION,  ELSEWC6 
.060 

♦MATERIAL,  ELSET-PCS 

♦DENSITY 

8.618  0  IE  —A 

♦ELASTIC,  TYPE-ISO 

2.5=6, .12 

♦BOUNDARY 

EDGE, 2, 3 

E  D GE ,5,6 

SIDE.l 

SIDE, 3, A 

SI DE  ,  6 

♦STEP,  L  I  NE  Art -Nt  * 

N4TURAL  F  RE  OUENC I E  S  AND  MODE  SnAPeS  OF  A  PC  BOARD 

♦FRE3UENCY 

5,2000 

♦PRINT 

•END  STEP 


5  10  lb  20  2b  30  3b  AO  A5  50  55  60  65 


Figure  23.  ABAQUS  Input  Listing  for  Problem  #2 
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As  indicated  in  Table  20,  STARDYNE  appears  to  most 
accurately  predict  the  first  four  natural  frequencies.  ANSYS  is 
the  second  most  accurate  of  the  three,  with  ABAQUS  being  the 
least  accurate.  The  percent  deviation  of  each  program  from  the 
theoretical  solution  is  summarized  in  Table  21. 


Table  20. 

Natural  Frequencies  of  PC 

Board  (Hz) 

Mode 

STARDYNE* 

ANSYS 

ABAQUS 

Theoretical 

1 

120  .0 

115  .1 

124  .2 

118  .1 

2 

300  .4 

285 .5 

377  .2 

295  .2 

3 

300  .4 

285  .5 

377  .2 

295  .2 

4 

484 .7 

438  .8 

572  .6 

472  .4 

*Lanczos  Method;  other  options  are  Househol der-QR  Method  or 
Inverse  Iteration  Method 


Table  21  . 

Deviation  from 

Theoretical 

Fre  quencies 

Mode 

STARDYNE 

ANSYS 

ABAQUS 

1 

+  1  .61% 

-2  .54% 

+  5  .17% 

2 

+  1  .76% 

-3  .29% 

+  27  .78% 

3 

+  1  .76% 

-3  .29% 

+  27  .78% 

4 

+  2  .60% 

-7  .11% 

+  21  .21% 

Average 

+  1  .93% 

-4  .06% 

+  20  .49% 

At  this  point,  it  would  be  timely  to  make  two  points 
concerning  idiosyncrasies  of  both  ANSYS  and  ABAQUS.  First,  when 
ANSYS  is  requested  to  determine  natural  frequencies,  it  will 
compute  as  many  natural  frequencies  as  there  are  DOF  in  the 
structure.  For  instance,  the  model  of  the  PC  board  had  69 
active  DOF;  69  frequencies  were  thus  calculated,  regardless  of 
the  fact  that  we  were  interested  in  only  the  first  four 
frequencies.  The  only  way  to  reduce  the  number  of  calculated 
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frecuencies  is  to  compose  a  list  of  "master  degrees  of  freedom"; 
these  are  the  DCF  considered  to  be  the  most  significant 
contributors  to  the  structural  behavior.  These  master  DOF  must 
be  carefully  selected  in  order  to  ensure  an  accurate  modal 
analysis.  For  our  analysis  of  the  PC  board,  all  69  DOF  were 
allowed  to  contribute  to  the  modal  analysis,  thus  remaining 
consistent  with  the  STARDYNE  and  ABAQUS  solutions. 

Secondly,  ABAQUS  requires  that  more  natural 
frecuencies  be  recuested  than  the  actual  number  which  are  of 
interest  to  the  analyst.  Unfortunately,  this  is  not  indicated 
in  the  manual.  The  first  time  that  the  ABAQUS  model  was  run, 
four  natural  frecuencies  were  recuested.  The  first  three 
frecuencies  were  the  same  as  those  indicated  in  Table  20; 
however,  the  fourth  frequency  was  computed  to  be  1293.9  Hz. 
Since  no  obvious  errors  could  be  found  in  the  model,  ABAQUS1 
customer  service  was  consulted.  They  volunteered  to  build  and 
run  their  own  PC  board  model  with  ABAQUS  to  determine  the  source 
of  the  problem.  After  several  iterations  it  was  determined  that 
there  were  no  model  errors,  but  that,  instead,  the  last 
eigenvalue  calculated  by  ABAQUS  was  not  an  accurate  solution. 
In  order  to  determine  four  natural  frequencies  accurately,  it 
was  necessary  to  request  a  solution  of  five  or  more  frequencies. 
As  a  general  rule  of  thumb,  it  was  suggested  that  the  analyst 
should  always  request  two  or  three  more  natural  frequencies  than 
those  in  which  be  has  an  interest  when  using  the  modal  analysis 
option  of  ABAQUS. 

As  a  footnote  to  the  above  discussion  of  the  modal 
analysis  option  of  ABAQUS,  it  is  of  interest  to  note  that 
STARDYNE  automatically  calculates  more  natural  frequencies  than 
recuested  (how  many  more  is  a  function  of  the  original  number  of 
frecuencies  requested).  These  additional  frecuencies  are  then 
presented  to  the  analyst  as  "approximate"  eigenvalues,  without 
their  corresponding  mode  shapes.  For  example,  in  this  PC  board 


analysis,  four  frequencies  were  requested.  STARDYNE  then 
calculated  four  natural  freouencies  and  mode  shapes,  and 
"estimated"  two  additional  natural  frequencies.  In  this  way, 
STARDYNE  ensures  the  accuracy  of  the  requested  eigenvalues  and 
ei genvectors  . 

Returning  to  the  discussion  of  the  PC  board  results, 
the  next  comparison  which  could  be  made  between  the  programs 
concerns  the  computed  mode  shapes.  Table  22  summarizes  the 
normalized  mode  shapes  for  each  of  the  three  programs.  Note 
that  the  tabulated  values  represent  normalized  deflections  in 
the  transverse  direction.  Also  notice  that  only  modes  1  and  4 
were  calculated  to  have  the  same  normalized  shape  by  all  three 
programs;  differences  do  exist  among  the  calculated  shapes  for 
modes  2  and  3  . 

As  shown  in  Figures  24  through  27,  the  mode  shapes 
calculated  by  each  of  the  three  programs  are  basically  the  same. 
The  magnitudes  shown  in  Table  22  indicate  that  differences  do 
exist  between  the  programs  for  modes  2  and  3;  however,  the 
sketches  of  the  mode  shapes  reveal  that  the  basic  shape 
descriptions  are  the  same.  It  is  interesting  to  note  that  mode 
2  and  mode  3  are  the  same  mode,  except  for  being  rotated  90° 
around  the  PC  board  from  each  other. 

In  summary  of  Sample  Problem  #2,  STARDYNE  was  the 
most  accurate  program  for  determining  the  first  four  natural 
frequencies.  It  also  had  the  most  direct,  us e r- o r i en te d 
procedures  for  requesting  these  frequencies  --  no  i dosyncras i es 
were  apparent  as  was  the  case  for  both  ANSYS  and  ABAQUS  . 

5  .5  Sample  Problem  #3  -  Heat  Transfer 

The  heat  transfer  sample  problem  involved  the 
calculation  of  the  temperature  distribution  and  heat  flux  of  a 
0.6"  long  by  0.1"  diameter  aluminum  pin  fin.  The  base  of  the 
fin  was  assumed  to  be  held  at  a  constant  temperature  of  176°F, 
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Table  22.  PC  Board  Normalized  Mode  Shapes 


Mode 

Node 

STARDYNE 

ANSYS 

ABAQUS 

1 

7 

0  .500 

0  .500 

0  .500 

8 

0  .707 

0.707 

0  .707 

9 

0  .500 

0  .500 

0  .500 

12 

0.707 

0.707 

0  .707 

13 

1  .000 

1  .000 

1  .000 

14 

0  .707 

0  .707 

0  .707 

17 

0  .500 

0  .500 

0  .500 

18 

0  .707 

0.707 

0.707 

19 

0  .500 

0  .500 

0  .500 

2 

7 

-0  .414 

-0  .706 

0  .496 

8 

0  .414 

0  .002 

1  .000 

9 

1  .000 

0 .706 

0.918 

12 

-1  .000 

-1  .000 

-0  .298 

13 

0  .000 

0  .000 

0  .000 

14 

1  .000 

1  .000 

0  .298 

17 

-1  .000 

-0  .706 

-0  .918 

18 

-0  .414 

-0  .002 

-1  .000 

19 

0  .414 

0  .706 

0  .496 

3 

7 

1  .000 

-0  .706 

0  .327 

8 

1  .000 

-1  .000 

0  .938 

9 

0  .414 

-0 .706 

1  .000 

12 

0  .414 

-0  .002 

-0  .476 

13 

0  .000 

0  .000 

0  .000 

14 

-0  .414 

0  .002 

0  .476 

17 

-0  .414 

0  .706 

-1  .000 

18 

-1  .000 

1  .000 

-0  .938 

19 

-1  .000 

0  .706 

-0  .327 

4 

7 

-1  .000 

-1  .000 

-1  .000 

8 

0  .000 

0  .000 

0  .000 

9 

1  .000 

1  .000 

1  .000 

12 

0  .000 

0  .000 

0  .000 

13 

0  .000 

0  .000 

0  .000 

14 

0  .000 

0  .000 

0  .000 

17 

1  .000 

1  .000 

1  .000 

18 

0  .000 

0  .000 

0  .000 

19 

-1  .000 

-1  .000 

-1  .000 
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and  the  heat  transfer  coefficient  between  the  fin  and  the 
surrounding  air  flow  was  assumed  to  have  a  constant  value  of 
12.7  Btu/hr . f •  F .  Also,  it  was  assumed  that  the  temperature 
distribution  was  one-dimensional  along  the  length  of  the  fin; 
the  temperature  across  any  given  cross-section  was  assumed  to  be 
constant . 

It  was  originally  planned  to  solve  this  problem 
using  both  ANSYS  and  ABAQUS  (STARDYNE  does  not  have  any  heat 
transfer  elements).  Unfortunately,  ABAQUS  did  not  lend  itself 
well  to  solving  this  one-dimensional  problem;  the  reasons  will 
be  discussed  later  in  this  section  of  the  report.  Therefore, 
only  a  comparison  of  the  theoretical  closed-form  solution  and 
the  ANSYS  numerical  solution  was  made.  The  ANSYS  finite  element 
model  of  the  pin  fin  is  shown  below  in  Figure  28. 


Node  1 

Base  Temperature= 
1 76°F 


Air  Temperature  =  131°F 


Figure  28.  Pin  Fin  Finite  Element  Model 
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As  shown  in  Figure  28,  the  pin  fin  finite  element 
del  consisted  of  14  nodes  and  13  elements.  Modes  1-7  were 
used  to  define  the  pin  fin  structure,  while  nodes  8-14  were  used 
to  represent  the  ambient  airflow  "sink”.  It  should  be  noted 
that  nodes  8-14  could  very  easily  be  replaced  by  a  single  node 
to  which  the  convective  link  elements  7-13  would  be  connected. 
At  the  time  of  the  model  creation,  it  was  felt  that  having 
distinct  ambient  nodes  would  simplify  the  generation  of  the 

elements  7-13.  However,  this  slight  advantage  in  element 

generation  was  partially  offset  by  the  need  to  create  the 
additional  nodes  9-14. 

Elements  1-6  are  one-dimensional  conduction  elements 
having  the  conductivity  of  aluminum.  Element  7  is  a  convective 
link  element  with  a  heat  transfer  coefficient  of  12.7  Btu/hr  ft^ 
°F,  and  a  heat  transfer  area  eoual  to  the  surface  area  of  a 

0.05"  long  by  0.1"  diameter  cylinder.  Similarly,  elements  8-12 
are  convective  link  elements,  but  their  heat  transfer  areas  are 
each  ecual  to  the  surface  area  of  a  0.1"  long  by  0.1"  diameter 
cylinder.  Finally,  the  last  convective  link  element,  number  13, 
has  a  heat  transfer  area  ecual  to  the  sum  of  the  surface  area  of 

a  0.05"  long  by  0.1"  diameter  cylinder  and  the  surface  area  at 

the  end  of  the  pin  fin. 

As  mentioned  earlier,  ABAQUS  was  not  used  to  model 
the  pin  fin.  The  primary  reason  for  the  decision  not  to  use 
ABAQUS  was  that  ABAQUS  did  not  provide  a  convective  link  element 
which  could  be  used  in  conjunction  with  its  one-dimensional 
conductive  link  element  "DC1D2".  Although  ABAQUS  does  provide 
for  convection  heat  transfer  away  from  its  one-d imens i ona 1 
conduction  element,  this  convection  occurs  at  each  face  of  the 
element  (a  face  being  defined  as  the  cross-sectional  surface 
area  at  a  node)  .  As  may  be  seen  by  referring  to  the  previous 
paragraph,  the  cross-sectional  area  of  the  fin  is  not  the 
convective  heat  transfer  area  required  for  this  model. 


101 


A  method  of  solution  for  this  problem  using  ABAQUS 
was  available  using  the  program's  4-node,  axisymmetric  heat 
transfer  element  "DCAX4"  .  However,  this  would  have  resulted  in 
a  two-dimensional  solution  to  the  problem,  which  would  not  have 
provided  a  totally  valid  comparison  with  the  one-dimensional 
solution  of  ANSYS.  Also,  we  felt  that  a  problem  of  this 
simplicity  did  not  warrant  a  second  solution  by  ANSYS  using 
two-dimensional  heat  transfer  elements  solely  for  the  purpose  of 
comparison  with  the  two-dimensional  ABAQUS  solution.  Therefore, 
only  ANSYS  was  used  to  solve  this  thermal  problem  (a  listing  of 
the  recuired  ANSYS  input  is  shown  in  Figure  29). 

As  indicated  in  the  input  data  listing,  the  ANSYS 
"STIF32"  element  was  used  to  represent  the  pin  fin,  and  the 

ANSYS  "STIF34"  element  was  used  to  provide  the  convection  link. 
One  should  not  be  confused  by  the  fact  that  ANSYS  labels  its 

one-dimensional  heat  conducting  bar  (STIF32)  as  a  two-dimension¬ 
al  element.  By  using  the  term  "two-dimensional",  ANSYS  is 
referring  to  the  fact  that  this  element  is  used  only  in  planar 
or  axisymmetric  analyses  where  each  node  is  defined  by  two 
coordinates.  ANSYS  also  has  a  one-dimensional  heat  conducting 
bar  that  is  labeled  as  "three-dimensional":  both  of  the 

elements’  nodes  are  defined  by  a  three-dimens ional  coordinate 
system . 

5.5.1  Results  for  Sample  Problem  #3 

As  a  basis  of  comparison  for  the  ANSYS  solution,  the 
temperature  distribution  was  also  determined  based  on  the 

theoretical  solution.  The  closed-form  solution  was  taken  from 
the  third  edition  of  Principles  of  Heat  Transfer  by  Frank  Kreith 
[87],  and  is  presented  below: 
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LND 

F IN  I  SM 


T(y)  -  T  .  _  cosh  m(L~y)  +  (h/mk)  sinh  m(L-y) 

Tbase  '  Tair  cosh  +  sinh  (mb) 


here:  T(y)  =  temperature  at  distance  "y"  from  the  base  (°F) 

Tair  =  air  temperature  ( 1 3 1 0 F > 

Tbase  =  base  temperature  (176°F) 

m  _  /hP 

m  -  \/TJ 

h  =  heat  transfer  coefficient  (12.7  Btu/hr  fx  °F) 

P  =  pin  fin  perimeter  (0.02618  ft) 

K  =  Conductivity  (96  Btu/hr  ft  °F) 

-  fi  2 

A  =  pin  fin  cross-sectional  area  (54.54  x  10  ft  ) 
L  =  pin  fin  length  (0.05  ft) 

In  addition,  the  closed  form  solution  for  the  heat  transfer  rate 
through  the  fin  is  again  taken  from  Kreith: 


c  =v/PhAK ~  T 


sinh  (mL)  +  (h/mK)  cosh  (mL) 


Table  23  summarizes  the  temperature  distribution  along  the  fin 
at  0.1"  intervals  (and  also  indicates  the  heat  flow  rate)  as 
determined  by  ANSI'S  and  by  the  closed-form  solution. 


Table  23.  Temperature  Distribution  and  Heat  Flow 


Distance  from  Base 

ANSYS 

Theoretical 

0.0" 

176  .00° 

F 

176.00°  F 

0  .1" 

174  .93° 

F 

174.92°  F 

0  .2" 

174  .04° 

F 

174.04°  F 

0  .3” 

173  .35° 

F 

173.35°  F 

0  .4" 

o 

in 

CO 

CM 

F 

172  .85°  F 

0  .5" 

172  .53° 

F 

172  .53°  F 

0  .6" 

1  72  .39° 

F 

1  72  .39°  F 

Heat  Flow  Rate: 

0  .7377  ! 

Btu/Hr 

0  .7374  Btu/Hr 
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As  shown  in  Table  23,  the  temperature  distribution 
computed  by  ANSYS  is  virtually  identical  to  the  temperature 
distribution  defined  by  the  closed-form  solution.  Also,  ANSYS' 
computed  heat  flow  rate  is  within  0.04*  of  the  theoretically 
correct  value.  Obviously,  the  one-dimensional  heat  transfer 
element  of  ANSYS  is  capable  of  sufficient  accuracy  for  this  type 
of  analysis. 

This  concl udes  the  discussion  of  the  third  sample 
problem.  The  next  section  of  this  report  will  deal  with  the 

last  sample  problem:  a  nonlinear  statics  problem. 

5  .6  Sample  Problem  #4  -  Nonlinear  Statics 

As  discussed  in  the  sample  problem  " Introducti on" , 
the  nonlinear  statics  problem  involved  the  determination  of  the 
deflections  for  an  electronic  substrate.  This  substrate  was 

assumed  to  be  mounted  in  a  module  such  that  a  0.012"  gap  existed 
between  the  bottom  surface  of  the  substrate  and  the  top  surface 
of  the  module  floor.  The  substrate  was  assumed  to  be  mounted  by 
screws  at  each  of  its  four  corners,  and  was  subjected  to 
constant  accelerations  of  (+)  and  (-)  20,000  G's.  The 

nonlinearity  of  this  problem  falls  into  a  class  of  problems 
called  "one-way  structures".  The  substrate  is  tree  to  deflect 
away  from  the  module,  but,  when  the  acceleration  is  reversed, 
the  substrate  may  only  deflect  freely  until  it  encounters  the 
module.  At  this  point,  further  distortion  of  the  substrate 
occurs  until  the  restraining  forces  equal  the  applied  load. 

All  of  the  following  three  programs  are  capable  of 

solving  this  type  of  nonlinear  problem:  (1)  STARDYNE,  (2) 

ANSYS,  and  (3)  ABAQUS  .  However,  each  program  requires  the 
analyst  to  go  through  different  solution  procedures.  The  finite 
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element  model  shown  in  Figure  30  was  used  by  each  of  the 
programs  to  solve  for  the  substrate  deflections. 

As  shown  in  Figure  30,  the  substrate  finite  element 
model  consisted  of  25  nodes  and  16  elements.  It  should  be  noted 
that  in  addition  to  the  nodes  and  elements  shown  in  Figure  30, 
each  of  the  three  programs  also  required  several  other  "gap" 
nodes  and/or  "gap"  elements  peculiar  to  each  program;  these  will 
be  described  to  a  limited  degree  later  in  this  section  of  the 
report  (a  complete  description  of  the  necessary  gap  elements 
and/or  nodes  for  each  program  is  given  in  the  corresponding  user 
manuals)  . 

The  following  Table  24  summarizes  the  physical  and 
material  properties  of  the  substrate  required  for  this  analysis. 


Table  24. 

Physical  Dimensions 

Plate  Type 

Material 

Young  1  s  Modul  us 

Poisson’s  Ratio 

Wei ght 

Restraints 

Loading 


Physical  and  Material  Properties 

1  .0"  x  1  .0"  x  0  .020" 

Bending  and  Membrane 
771  Alumina 
43  x  106  psi 
0  .3 

0.0028  lb 

Each  corner  fully  restrained 
20,000  G' s 


As  shown  in  Table  24,  the  plates  used  for  this 
analysis  had  both  bending  stiffness  and  membrane  stiffness.  In 
each  case,  the  quadri 1 ateral  plate  element  which  was  used  had 
four  nodes  per  element.  The  element  used  for  the  STARDYNE 
analysis  was  the  "QUAD"  element,  the  element  used  for  the  ANSYS 
analysis  was  the  "STIF43"  element,  and  the  element  used  for  the 
ABAQUS  analysis  was  the  "S4R"  element.  Again,  although  both 
ANSYS  and  ABAQUS  offer  e i g h t- no de - pe r- e 1 emen t  quadrilateral 
plates  elements,  the  more  simple  fo ur- node- per-el ement  plate 


i 

I 
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elements  were  used  in  order  to  provide  a  valid  comparison  with 
STARYDNEs  four-node  cuadilateral  element. 

A  listing  of  the  required  input  for  each  of  the 

three  programs  is  shown  in  Figures  31  -  35.  It  should  be  noted 
that  STARDYNE  reouired  a  total  of  three  different,  seouential 
runs  to  solve  this  problem;  hence,  the  three  separate  input 
listings.  Also,  the  ANSYS  input  listing  is  only  one  of  many 

different  ANSYS  runs  which  were  required  to  determine  the 
model's  necessary  "gap"  element  stiffness  for  the  downward 
loading  case.  Finally,  a  representative  ABAQIIS  input  is  shown. 

5.6.1  Solution  Technicues  for  Sample  Problem  #4 

The  first  portion  of  this  section  of  the  report  will 
be  divided  into  three  main  topics,  each  will  be  devoted  to  a 

discussion  of  the  solution  procedure  required  by  the  different 
programs.  It  should  be  mentioned  that  a  closed-form  solution 

did  not  exist  for  this  problem;  therefore,  the  numerical 
solutions  could  not  be  compared  with  a  theoreti cal  ly  correct 
sol  uti on  . 

5. 6. 1.1  STARDYNE  Sol ution 

STARDYNE's  "NUBOP"  program  was  used  to  solve  for  the 
deflections  of  the  substrate  when  exposed  to  the  acceleration 
loadings.  NUBOP  is  used  to  obtain  static  solutions  for  models 
which  have  nonlinear  connections,  boundaries,  or  elements. 
Typical  uses  of  NUBOP  include  " bo t tom i ng- o ut "  or  separation  of 
adjacent  structural  members,  or  the  analysis  of  tension-only  or 
compr es s i on- on  1 y  bar  members. 

In  general  terms,  the  following  procedure  was 
reouired  to  run  the  NUBOP  program.  First,  the  model  of  the 
substrate  was  examined  for  potential  "bottom- out"  points 
(referred  to  as  a  "BOP").  It  was  assumed  that  any  of  the 
following  nodes  might  possibly  bottom-out  against  the  module 
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Figure  31.  First  STARDYNE  Input 


111 


Listing  for  Problem 
(continued) 
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Figure  35.  ABAQUS  Input  Listing  for  Problem  f4 
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surface:  nodes  7-9,  12-14,  and  17-19.  Therefore,  each  of  these 
nodes  was  defined  as  a  "BOP". 

Next,  static  load  cases  were  run  for  the  substrate 
model  with  upward  unit  loads  applied  separately  to  each  of  the 
ECP  nodes  (it  was  known  that  the  module  surface  would  resist  the 
downward  deflection  of  the  substrate  with  an  upward  force) . 
These  unit  static  load  cases,  referred  to  as  BOP  cases,  were 
then  used  by  the  NUBOP  program  to  create  a  flexibility  matrix 
for  the  substrate.  In  addition,  a  pair  of  load  cases  were  run 
in  which  the  two  accel erati ons  of  interest  ( ^2  0,000  G‘s)  were 
applied  to  che  structure.  This  resulted  in  the  creation  of  a 
pair  of  vectors  representing  the  relative  displacements  at  each 
BOP  node  due  to  the  two  accelerations. 

The  third  step  in  the  NUBOP  procedure  was  to  set  the 
displacement  bounds  for  each  BOP  node  in  the  displacement 
vectors.  For  our  substrate  problem,  each  BOP  node  was 
restricted  to  a  maximum  downward  travel  of  0.012". 

Finally,  to  determine  the  deflection  shapes  of  the 
substrate  when  exposed  to  the  ^+20,000  G  accelerations,  a  static 
load  case  was  run  which  combined  the  applied  acceleration  loads 
with  the  necessary  BOP  loads.  It  should  be  notea  that  when  the 
substrate  deflected  under  the  -20,000  G  acceleration  (which  is 
eeui valent  to  an  upward  load  on  the  substrate),  the  displacement 
bound  was  not  crossed,  and  no  BOP  loads  were  allowed  to 
contribute  to  the  final  load  vector. 

5. 6. 1.2  ANSYS  Sol  uti  on 

ANSYS  did  not  recuire  the  series  of  procedural  runs 
necessary  for  the  STARDYNE  solution,  but  ANSYS  did  require  that 
nine  additional  nodes  and  elements  be  added  to  the  finite 
element  model.  The  additional  elements  (ANSYS  "STIF52" 
elements)  were  nonlinear  interface  elements  capable  of 
supporting  only  compression  in  the  direction  normal  to  the 
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surfaces  and  shear  in  the  tangential  directions.  This  element 
may  be  given  a  gap  specification,  and  a  specified  stiffness  acts 
in  the  normal  and  tangential  directions  when  the  gap  is  closed. 

The  ST3F52  elements  required  for  this  analysis  were 
each  defined  by  a  unique  pair  of  nodes.  The  upper  node  for  each 
of  these  elements  was  already  in  existence  (nodes  7-9,  12-14, 
and  17-19);  however,  nine  additional  nodes  were  reeuired  to 
represent  the  module  surface  (these  nodes  were  fully 
restrai ned)  . 

Once  the  gap  elements  had  been  defined,  it  was 
necessary  to  specify  a  gap  interface  stiffness.  The  manual 
suggested  that  for  most  problems  the  local  surface  deformation 
is  not  of  importance,  and  that  the  stiffness  of  the  interface 
may  be  estimated  as  an  o-'der  of  magnitude  or  two  greater  than 
the  adjacent  element  stiffness.  However,  it  also  warned  against 
the  use  of  unreasonably  high  stiffness  values  due  to  large 
increases  in  the  iterative  solution  time.  Ideally,  we  would 
have  liked  the  interface  stiffness  to  be  infinite  so  that  the 
solution  could  be  compared  directly  with  STARDYNE's  solution  (in 
which  the  BOP's  were  infinitely  stiff  when  they  encountered  the 
transverse  deflection  limit  of  -0.012").  Therefore,  after 
running  the  program  several  times,  an  interface  stiffness  of  1  x 

Q 

10  lb/in  was  finally  decided  upon  as  a  value  approaching  the 
maximum  stiffness  that  would  still  allow  static  convergence  of 
the  sol  uti on  . 

5. 6. 1.3  ABAQUS  Sol  uti  on 

The  ABAQUS  solution  procedure  is  probably  the 
simplest  of  the  three  programs'  methods.  Nine  additional  nodes 
were  defined  0.012"  below  the  nine  potential  bottom-out  nodes  on 
the  substrate.  Next,  while  specifying  the  boundary  conditions, 
the  gap  condition  was  defined  between  the  nine  pairs  of 
potential  interface  nodes.  At  this  point,  all  that  remained  was 
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to  indicate  to  A6AQUS  that  a  nonlinear  statics  run  was  to  be 
initiated  (several  iteration  and  convergence  parameters  were 
required  by  the  program) . 

Unf ortunately ,  the  ABAQUS  runs  were  consistently 
aborted  from  the  system.  Finally,  the  ABAQUS  Customer  Service 
analysts  were  consulted  with  our  problem.  After  several 
telephone  conversations,  it  was  concluded  that  the  S4R 
quadrilateral  plate  elements  which  had  been  used  to  model  the 
substrate  could  not  be  used  for  this  model  without  resulting  in 
nodal  singularities.  The  S4R  element  could  not  handle  the 
boundary  condition  restraints  imposed  upon  the  substrate  (all 
four  corners  fully  restrained)  ever,  for  a  linear  static  load 
case.  If  we  wished  to  run  this  model,  it  would  have  been 
necessary  to  change  the  elements  to  "S8R"  elements  (quadri¬ 
lateral  plate  elements  with  mid-side  nodes).  However,  we 
decided  that  these  higher  order  elements  would  not  have  provided 
a  valid  comparison  with  the  f our-node-per-el ement  plate  elements 
used  by  both  STARDYNE  and  ANSYS.  Therefore,  no  solution  was 
found  using  ABAQUS  for  the  nonlinear  statics  problem. 

5  .6  .2  Results  from  Sample  Problem  #4 

The  following  Table  25  summarizes  the  deflections 
due  to  both  the  -20,000  G  acceleration  and  the  +20,000  G 
acceleration,  as  determined  by  STARDYNE  and  ANSYS. 

As  shown  in  Table  25,  the  STARDYNE  and  ANSYS 
solutions  differ  by  a  considerable  amount  (it  should  be  noted 
that  only  one-quarter  of  the  model  is  represented  in  Table  25, 
due  to  the  symmetry  of  the  substrate).  Unfortunately,  even  the 
linear  part  of  the  solution  ("OPEN  GAP"  column)  differs  between 
the  programs.  A  clue  to  which  of  the  "linear"  solutions  is  the 
more  accurate  might  lie  in  the  fact  that  the  deflections  for  the 
symmetric  node  pairs  (6  8  2,  11  &  3,  and  12  8  8)  of  the  ANSYS 
model  are  not  identical,  as  they  are  for  the  STARDYNE  model. 
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Table  25.  Transverse  Deflections  of  Nonlinear  Substrate  Model 

( inches} 

STARDYNE  ANSYS 


Node 

Open  Gap 

Close  Gap 

★ 

Open  Gap 

Close  Gap* 

2 

0  .00476 

-0 .00403 

0  .00626 

-0  .00595 

3 

0  .00852 

-0 .00718 

0  .01055 

-0  .01001 

6 

0  .00476 

-0 .00403 

0 .00625 

-0.00526 

7 

0  .01025 

-0 .00825 

0  .01313 

-0  .01167 

8 

0  .01320 

-0 .01048 

0 .01621 

-0.01451 

11 

0  .00852 

-0 .00718 

0  .01054 

-0  .00875 

12 

0  .01320 

-0.01048 

0 .01626 

-0.01413** 

13 

0  .01566 

-0  .01200 

0  .01885 

-0  .01650** 

*  Refers 

to  the  case 

where  the 

acceleration  attempts 

to  close 

the  initial  gap  (+20,000  G) 

**  Exceeds  initial  gap  clearance 

This  could  indicate 

that  the 

ANSYS 

sol uti on 

had 

not  fully 

converged 

for  this 

parti c  ul ar 

r  un  . 

Further 

iterations  with 

various 

convergence 

criteria 

might 

possibly 

correct  this 

discrepancy.  However,  based  on  the  solutions  presented  in  Table 
25,  it  is  believed  that  the  STARDYNE  solution  is  probably  the 
more  reliable  of  the  two. 

When  comparing  the  "closed  gap"  solutions,  it  was 
immediately  obvious  that,  in  the  ANSYS  solution,  the  maximum 
deflection  was  not  limited  to  -0.012  inches.  This  was  due  to 
the  fact  that  the  interface  element  did  not  have  infinite 
stiffness;  when  the  substrate  encountered  the  module  surface 
nodes,  it  actually  deformed  the  module  surface  to  some  degree. 
However,  we  had  hoped  that  the  large  stiffness  that  had  been 
given  to  the  interface  elements  would  have  prevented  this  from 
happening.  It  might  be  possible  to  increase  the  interface 
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stiffness  (and  then  increase  the  maximum  number  of  iterations  to 
allow  the  solution  to  converge)  in  order  to  more  closely  model 
the  real  world  situation. 

The  STARDYNE  "closed  gap"  solution  appears  to  be  a 
very  believeable  solution:  the  maximum  deflection  is  -0.012 

inches,  and  the  symmetric  node  pairs  have  identical  deflections. 
In  addition  t.o  being  the  more  believable  solution,  it  must  be 
noted  that  the  STARDYNE  solution  was  very  straightforward.  The 
STARDYNE  solution  worked  the  first  time  that  it  was  run,  and 
there  was  no  uncertainty  about  interface  stiffness  or 
convergence  criteria.  STARDYNE  appears  to  have  the  most 
user-oriented  nonlinear  "gap"  routine  of  the  three  programs  that 
were  considered.  This  concludes  the  discussion  of  the  nonlinear 
statics  sample  problem.  The  next  section  of  this  report  will 
present  a  summary  of  the  program  capabilities  and  idiosyncrasies 
which  became  apparent  during  the  sample  problem  analyses. 

5  .7  General  Summary  for  Sample  Problems  Analysis ( Phase2  ) 

This  summary  will  explain  the  major  conclusions  from 
each  of  the  sample  problem  analyses,  and  it  will  also  list  the 
more  general  observations  of  program  idiosyncrasies  and/or 
useful  features  which  were  discovered  during  the  sample  problem 
analysis  . 

5 .7 .1  Linear  Statics  Problem 

(1)  ABAQUS  most  closely  approximated  the  theoretical 
maximum  deflection  at  the  center  of  the  lid. 
However,  it  should  be  noted  that  all  solutions 
(STARDYNE,  ANSYS,  and  ABAQUS)  were  within  31  of 
the  theoretical  answer. 

(2)  ABAQUS  also  came  nearest  to  approximately  the 
theoretical  maximum  stress  in  the  center  of  the 
simply  supported  lid,  followed  closely  by  ANSYS 
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and  STARDYNE.  All  stresses  were  within  8%  of 
the  theoretical  answer;  this  error  is  partially 
a  function  of  the  coarse  nodal  gridwcrk  used  in 
our  model  . 

(3)  The  deflection  shapes  calculated  by  STARDYNE  and 
A  N  S  Y  S  compare  well  with  each  other.  The 
deflection  shape  calculated  by  ABAQUS  o'viates 
slightly  f^om  the  other  two. 

Linear  Dynamics  Problem 

(1)  STARDYNE  most  accurately  predicted  the  first 
four  natural  freouencies  of  the  circuit  board 
(average  deviation  of  1.9*),  followed  by  ANSYS 
(average  deviation  of  4.1%),  and  lagged  con¬ 
siderably  by  ABAQUS  (average  deviation  of 
20  .5%)  . 

(2)  STARDYNE's  Lanczos  routine  was  the  most  user- 
oriented  modal  analysis  routine  of  the  three 
programs;  it  accurately  calculated  the  desired 
number  of  natural  frecuencies  and  mode  shapes, 
and  estimated  several  natural  frecuencies  beyond 
those  re  cues  ted  . 

(3)  ANSYS  calculated  a  natural  frecuency  for  every 
"Master-DCF"  in  the  structure.  However,  it 
provided  mode  shapes  only  for  the  first  four 
recuested  natural  frecuencies.  It  should  be 
noted  that  we  did  not  reduce  the  number  of  DOF 
in  this  model  to  a  few  select  DOF,  as  suggested 
by  the  ANS^s  manual.  Instead,  we  allowed  all 
the  structural  DOF  to  contribute  in  the  modal 
extraction  process  ;n  order  to  remain  consistent 
with  the  STARDYNE  and  ABAQUS  solutions. 
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(4)  ABAQUS  Initially  appeared  to  have  a  very 
straightforward  modal  analysis  routine. 
However,  ABAQUS  recuired  the  user  to  request 
more  natural  frequencies  and  mode  shapes  then 
the  four  of  interest  to  insure  the  accuracy  of 
the  higher  modes.  Unfortunately,  this  is  not 
discussed  in  the  manual;  large  errors  were 
evident  until  this  idiosyncrasy  was  discovered. 
Yet,  even  the  best  modal  solution  by  ABAQUS  did 
not  approach  the  accuracy  of  STARDYNE  or  ANSYS  . 

(5)  It  should  be  mentioned  that  STARDYNE  provides 
three  other  modal  extraction  methods  in  addition 
to  Lanczos  .  These  are  "Householder-QR," 
"Inverse  Iteration,"  and  "Householder-QR  with 
Guyan  Mass  Condensation."  The  "H-QR  with  Guyan 
Mass  Condensation"  method  is  somewhat  similar  to 
the  method  used  by  ANSYS  in  that  it  artificially 
reduces  the  number  of  DOF  in  the  structural 
model  to  a  few,  select  "master-DOF . "  For  a 
complete  understanding  of  these  extraction 
methods,  refer  to  the  STARDYNE  User  Information 
Manual  .[82] 

5.7.3  Heat  Transfer  Problem 

(1)  ANSYS  provided  an  accurate  solution  for  the 
temperature  distribution  and  the  heat  flux  using 
relatively  simply  one- di mens i ona 1  heat  conduc¬ 
ting  elements.  The  model  creation  was  very 
strai ghtf orward . 

(2)  ABAQUS  did  not  provide  an  obvious  method  of 
including  the  convective  heat  transfer  effects 
from  the  sides  of  the  pin  fin  --  at  least  not 
using  its  one-dimensional  heat  conducting  bar 
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element.  The  problem  could  have  been  solved 
with  its  two-dimensional  element  but  this  would 
not  have  provided  a  valid  comparison  with  the 
ANS  YS  solution. 

(3)  STARDYNE  has  no  heat  transfer  elements  in  its 
element  library,  and,  therefore,  could  not  be 
used  to  solve  this  problem. 

Non! i near  Statics  Problem 

(1)  The  STARDYNE  solution  method  was  relatively 
us e  r- c  •  i  e n  t e d  .  It  did  re cui re  three  separate 
proc  jural  computer  runs,  and  the  storage  and 
retrieval  of  data  on  files.  However,  the 
program  required  no  knowledge  of  convergence 
criteria,  the  program  ran  properly  the  first 
time  it  was  attempted,  and  the  solution  appears 
to  be  more  accurate  than  the  ANSYS  solution 
(considering  the  symmetry  of  the  calculated 
deflections,  and  the  maxim  urn  downward  deflection 
of  0 .012" ) . 

(2)  The  AK'.SYS  solution  was  not  difficult;  however, 
it  did  expect  the  first-time  user  to  have  some 
knowledge  of  the  proper  relation  between 
gap-stiffness  and  convergence  criteria.  The 
program  recuired  several  iterative  runs  before 
it  appeared  to  converge;  further  runs  could 
result  in  an  even  more  accurate  solution.  The 
solution  presented  in  this  report  does  not 
display  the  expected  deflection  symmetry,  nor 
does  it  limit  the  maxim  urn  downward  travel  of  the 
module  to  0.012".  Also,  a  comparison  of  the 
upward  deflection  shape  calculated  by  STARDYNE 
and  ANSYS  shows  a  considerable  deviation  for  the 
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linear  portion  of  this  problem.  It  is  believed 
that  the  ANSYS  solution  would  reouire  further 
manipulation  before  it  approaches  the  results  of 
the  STARDYNE  solution. 

(3)  The  ABAQUS  solution  outlined  in  the  user  manual 
appeared  to  be  very  user-oriented.  However,  it 
became  evident  that  the  S4R  element  used  to 
model  the  substrate  was  not  capable  of  solving 
even  the  linear  portion  of  this  problem  without 
aborting  from  the  CDC  system,  due  to  nodal 
singularities.  A  discussion  with  ABAQUS1 
customer  service  indicates  tha.  the  S3R 

Quadr i 1 ateri al  plate  element  (8  nodes  per  plate) 
should  be  able  to  solve  this  problem;  it  was  not 
used  because  it  would  not  have  provided  a  velia 
comparison  with  the  4-node  elements  of  STARDYNE 
and  ANSYS . 

General  Observations 

(1)  STARDYNE  offers  only  a  4  node  quadrilateral 
element;  bo*h  ANSYS  and  ABAQUS  offer  8  node 
quadri 1 ateral  elements. 

(2)  STARDYNE  and  ANSYS  require  a  fixed  format  input; 
ABAQUS  offers  a  free  format  option. 

(3)  ANSYS  offers  third  level  element  and  node 
generation;  STARDYNE  offers  first  level  element 
generation  and  second  level  node  generation; 
ABAQUS  provides  first  level  element  and  node 
generati on  only. 

(4)  For  modal  analyses,  STARDYNE  requires  a  weight 
density  while  ANSYS  and  ABAQUS  reouire  a  mass 
dens i ty  . 


(5)  For  plate  elements,  ANSYS  allows  the  user  to 
specify  a  small  rotational  stiffness  about  an 
axis  perpendicular  to  the  plate.  The  plate 
elements  of  STARDYNE  and  ABAQUS  do  not  have  any 
rotational  stiffness  in  this  direction;  the 
corresponding  DOF  must  be  restrained. 

6.0  CRITICAL  PROBLEM  ANALYSIS  AND  CORRELATION  TO  TEST  DATA 

(PHASE  3) 

6  . 1  I ntroducti on 

Throughout  this  study,  we  were  searching  for  a  critical 
problem  that  would  be  of  current  interest  to  the  microelectronic 
industry,  and  yet  one  with  existing  closed-form  and  test  results. 
From  the  survey  results  of  Section  3.0,  it  was  seen  that  lid  seal 
of  large  hybrid  packages  seemed  to  be  the  problem  which  concerned 
most  respondents.  Fortunately,  reference  [88]  contains  the 
closed-form  solution  and  test  results  for  this  type  of  problem. 
Libove  [88]  provides  formulas  for  the  maximum  tensile  stress  in  the 
lid-to-wall  seal,  the  maximum  lid  deflection,  and  the  lid 
collapsing  pressure  for  a  rectangular  flat-pack  under  external 
pressure  . 

Two  separate  problems  were  considered.  First,  finite 
element  analysis  was  applied  to  the  flat-pack  lid  problem  described 
in  "Example  1"  of  Libove's  report  [88],  Basically,  this  example 
involved  a  wide-seal,  uniform-wall,  constant-thickness-lid 
flat-pack  which  was  exposed  to  an  external  pressure  of  30  psi;  it 
was  necessary  to  solve  for  both  the  maximum  seal  stress  and  the 
maximum  deflection  of  the  lid. 

Second,  an  available  I s otron i cs/Dev i ce  Closures  package 
was  analyzed  using  the  finite  element  method  to  ^olve  for  the 
stresses  and  deflections  due  to  several  loading  conditions:  (1) 

external  pressure,  (2)  discrete  point  loads,  and  (3)  thermal  loads. 
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The  hybrid  which  was  considered  was  a  uniform-wall  package  with  a 
welded,  stepped-thickness  lid.  Our  numerical  results  for  the 
external  pressure  load  cases  were  compared  to  the  values  predicted 
by  Libove,  while  our  results  for  the  discrete  point  load  case  and 
the  thermal  load  case  were  compared  to  experimental ly-deri ved 
results  using  this  particular  package. 

6  .2  Critical  Problem  =1 

Critical  problem  -1  concerned  a  Kovar  flat-pack  which 
was  exposed  to  an  external  pressure  of  30  psi.  In  order  to  easily 
make  a  comparison  between  our  analysis  and  the  work  of  Libove,  the 
flat-pack  that  was  chosen  was  the  same  wide-seal,  uni f orm-wa 1 1  , 
constant-thickness-lid  Kovar  package  described  by  Libove  in  his 
"Example  1"  of  RADC-TR-79  -  1  38  .  This  package  is  further  described 
in  Figure  36. 

LID: 

E=  20x1 06  psi 
0.3 


E  =  20  x  106  psi 
v  =  0.3 

(wall s  are  0 . 04C  ’  thick; 

Eigure  36.  Critical  Problem  #1  Flatpack  Configuration 
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6.2.1  “First-Cut"  Finite  Element  Model 

To  begin  our  analysis  of  this  problem,  we  decided  tc 
generate  a  preliminary  “first-cut"  finite  element  model  which  could 
be  used  to  calculate  the  flat-pack  lid  deflections.  We  wished  to 
demonstrate  that  a  simple,  inexpensive  model  could  be  used  tc 
determine  the  maximum  lid  deflection,  and  that  this  numerical 
solution  would  compare  well  with  Libove's  analytical  results.  In 
all  of  the  following  analyses,  it  should  be  noted  that  the  STARDYNE 
finite  element  program  was  used  to  determine  the  recuired  solu¬ 
tions  . 

The  initial  finite  element  model  is  shown  in  Figures  37 
and  38.  It  should  be  noted  that  all  finite  element  model  views 
shown  in  this  critical  analysis  section  were  generated  with  the 
Unistruc  preprocessor.  The  model  was  a  “ouarter  model"  which  took 
advantage  of  both  the  structural  symmetry  and  the  load  symmetry  to 
simplify  the  analysis.  As  shown,  this  model  consisted  of  34  nodes 
and  24  cuadri  1  ateral  plate  elements;  it  recuired  only  about  one 
hour  to  generate  this  model  and  obtain  the  results. 

It  should  be  mentioned  again  that  the  purpose  of  this 
model  was  to  obtain  the  maximum  lid  deflection;  we  were  not  looking 
for  seal  stresses  at  this  time.  Therefore,  the  model  included 
only  the  package  walls  and  the  lid;  no  elements  were  included  to 
represent  the  lid  seal  between  the  walls  and  the  lid.  The  lid  was 
connected  to  the  package  walls  through  the  use  of  common  nodes. 

The  model  was  fully  restrained  at  the  nodes  along  the 
bottom  edges  of  the  walls.  This  was  a  realistic  assumption  since 
the  walls  are  effectively  cantilevered  from  the  bottom  of  the 
flat-pack  {assuming  a  flat-pack  bottom  thickness  of  at  least  0.040" 
thick,  the  rotational  stiffness  of  the  wall/bottom  j unction  is 
considerably  greater  than  the  rotational  stiffness  of  the  wall/lid 
j  uncti on )  . 

The  next  section  of  this  discussion  will  explain  our 
results,  and  will  compare  them  with  the  solution  presented  by 
L  i  b  o  v  e  . 
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Figure  37.  "First-Cut"  Finite  Element  Model 
(Node  Numbers  Shown) 
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6  .2. 1.1  " Fi rs  t-C  ut "  Results 

The  maximum  lid  deflection  and  lid  stress  calculated  by 
our  finite  element  analysis  are  compared  to  the  deflection  and 
stress  predicted  by  Libove  [88]  in  Table  26.  In  addition,  the 
maximum  lid  deflection  and  the  maximum  lid  stress  determined  by  our 
finite  element  analysis  are  compared  to  the  corresponding  values 
calculated  with  the  formulas  provided  by  Roark  for  a  clamped 
rectangular  plate  under  external  pressure  [85].  The  deflection  and 
stress  calculations  using  the  methods  of  Libove  and  Roark  are  shown 
bel ow: 

(1)  Libove  (refer  to  "Example  1"  of  RADC-TR-79-138 , 
page  30) 

( a )  Deflection 

=  12  (1-V2)  |  (f)3  (a)  (n4)  ( n 5 ) 

maximum  lid  deflection 
Poisson's  ratio  *  0.3 
pressure  =  30  psi 
Young's  Modulus  *  20  x  106  psi 
lid  width  =  0.92" 
lid  thickness  =  0.015" 
function  of  K  and  b/a  =  0.00125 
ratio  of  wall's  to  lid's  flexural  stiffness  = 
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lid  length  =  0.92" 

large  deflection  theory  correction  factor=0.981 
0  .00426" 

(b)  Stress 

2 

6  n,  P  a  ,  ,  .  ,L  , 

-  1 _  (page  21  of  Libove) 

t2 


<5 

max 

where:  6 

max 

V 

P 

E 


a 

t 


6 

max 
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where : 


where : 


where: 


max  =  maximum  lid  stress  at  middle  of  long  edge 
n  ^  =  function  of  K  and  b/a  =  0.051 

P  =  pressure  =  30  psi 

a  =  lid  width  =0.92" 

t  =  lid  thickness  =  0.015" 


c 

max 


34,533  psi 


(2)  Roark  (refer  to  Formulas  for  Stress  and  Strain, 
page  392) 

(  a )  Def 1 ecti on 


•^max 


ac  b 
ET 


4 

7 


^max 

c 

b 

E 

t 

a 

■^max 


maximum  plate  deflection 
pressure  =  30  psi 
pi  ate  width  =  0.92" 

Young's  modulus  =  20  x  106  psi  , 
plate  thickness  =  0.015" 
function  of  a/b  =  0.0138 
0  .00439" 


(b)  Stress 


B1  c  b< 


max 


max 
9  = 

b 
t 

B1  = 
max 


maximum  plate  stress  at  middle  of  long  edge 

pressure  =  30  psi 

plate  width  =  0.92" 

plate  thickness  =  0.015" 

function  of  a/b  =  0.3078 

34,736  psi 
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Table  26.  "First-Cut"  Results  for  Critical  Problem  *1 
Max  Lid  Deflection  Max  Lid  Stress 

Li bove  STARDYNE  Roark  Libove  STARDYNE*  Roark 

0.00426"  0.00442"  0.00439"  34,533  17,019  psi  34,736  psi 

*Von  Mises  combined  stress 

As  shown  in  Table  26,  the  maximum  lid  deflections  as 
calculated  by  Libove,  STARDYNE,  and  Roark  correlate  well;  STARDYNE 
calculates  a  deflection  which  is  only  3.7%  greater  than  the  value 
published  by  Libove.  Interestingly  enough,  the  solution  by  Roark 
also  compares  well  with  the  other  two  solutions  --  remember  that 
the  Roark  solution  was  based  on  an  analysis  of  a  totally  clamped 
rectangular  plate.  This  tends  to  support  the  Libove  proposal  that 
the  walls  of  this  flat-pack  are  effectively  clamping  the  edges  of 
the  lid. 

Also  shown  in  Table  26  is  a  comparison  of  the  maximum 
lid  stresses  as  calculated  by  Libove,  STARDYNE,  and  Roark.  As 
indicated,  Roark  and  Libove  agree  within  0.6%;  this  small 
difference  could  be  due  to  the  fact  that  Roark's  solution  assumes 
completely  clamped  edges,  while  Libove's  solution  considers  the 
edges  as  something  just  slightly  less  than  fully  clamped. 

However,  both  Libove  and  Roark  disagree  with  the 
STARDYNE  solution.  The  primary  reason  for  this  discrepancy  is  that 
the  finite  element  program  calculates  the  outer-surface  plate 
stresses  at  the  centroid  of  each  element.  Therefore,  although  it 
is  known  that  the  maximum  stresses  in  the  lid  will  occur  at  the 
midpoints  of  the  lid's  edges,  STARDYNE  considered  the  maximum 
stress  to  occur  at  the  centroid  of  elements  12  and  21  (the  elements 
nearest  the  edge  midpoints).  It  would  be  necessary  to  generate  a 
finer  element  grid  near  the  edge  midpoints  (nodes  14  and  30)  to 
approach  the  actual  stress  value  as  predicted  by  Libove  and  Roark. 
(It  should  be  remembered  that  the  primary  purpose  of  this  model  was 
to  calculate  deflections.) 
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One  advantage  of  the  STARDYNE  model  is  that  it  can 
consider  the  effect  of  the  pressure  applied  to  the  package  walls  as 
well  as  the  package  lid.  However,  our  analysis  indicated  that 
there  was  no  significant  difference  in  the  calculated  lid  deflec¬ 
tions  when  the  pressure  on  the  walls  was  included  in  this  model. 

The  next  section  of  this  discussion  will  describe  a 
more  detailed  finite  element  model  of  the  flat-pack  which  includes 
elements  that  represent  the  lid  seal. 

6  .2  .2  Finite  Element  Model  with  Lid  Seal 

To  continue  our  analysis  of  Critical  Problem  #1,  we 
felt  that  it  was  necessary  that  we  consider  a  finite  element  model 
of  the  flat-pack  which  included  the  lid  seal.  This  would  allow  us 
to  determine  the  maximum  seal  stress  in  addition  to  the  maximum  lid 
deflection  and  the  maximum  lid  stress.  Also,  we  were  interested  in 
determining  the  effect  of  different  seal  materials  on  the  deflec¬ 
tion  and  stress  values. 

This  finite  element  model  is  shown  in  Figure  39. 
Again,  the  model  is  a  "cuarter  model"  which  took  advantage  of  both 
the  structural  and  the  loading  symmetry.  As  indicated,  the  model 
consisted  of  112  nodes,  43  cube  elements,  and  43  Quadrilateral 
plate  elements.  The  cube  elements  were  used  to  model  the  package 
walls,  the  lid,  and  the  lid  seal.  The  43  thin  membrane  plate 
elements  (0.00001"  thick)  were  applied  to  the  top  and  bottom 
surfaces  of  the  lid's  cube  elements,  and  to  the  outside  surfaces  of 
the  wall  elements.  These  plate  elements  were  necessary  for  two 
reasons:  (1)  to  determine  the  maximum  lid  bending  stresses  on  the 
top  and  bottom  surfaces  of  the  lid  (the  stresses  calculated  for  the 
cube  elements  are  calculated  at  the  cube  centroids,  not  on  the 
exterior  faces),  and  (2)  to  allow  a  pressure  to  be  applied  to  the 
top  surface  of  the  lid  elements  and  to  the  exterior  surface  of  the 
wall  elements  (the  STARDYNE  program  does  not  allow  a  pressure  to  be 
applied  to  a  cube  element). 


134 


Figure  39.  Finite  Element  Model  with  Lid  Seal 


r. 
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As  previously  mentioned,  cube  elements  were  used  to 
model  the  lid  seal;  it  was  not  possible  to  accurately  model  the 
structural  geometry  with  plate  elements.  The  seal  was  assumed  to 
be  C.005"  thick,  and  as  wide  as  the  top  of  the  wall.  Three 
different  seal  materials  were  considered:  Kovar,  solder,  and 
glass.  The  necessary  properties  of  these  materials  are  listed  in 
Table  27  . 


Table  27.  Seal  Material  Properties 
Seal  Material  Young's  Modulus(E)  Poisson's  Ratio(V) 


Kovar 

20  x 

10°  psi 

0  .3 

Solder 

(62%  Sn ,  36%  Pb , 

6  .4 

,  n  6 

x  10  psi 

0.4 

2%  Ag) 

Glass 

(Corning  "9010") 

9  .8 

x  106  psi 

0  .22 

This  flat-pack 

model 

,  as  in 

the  previous 

"first-cut" 

model  , 

was  fully  restrained 

at 

the  nodes 

along 

the 

?dges  of  the 

bottom 

of  the  wal 1 s .  However, 

in  this 

case , 

the 

cubical  seal 

elements  were  used  to  connect  the  flat-pack  walls  to  the  lid 
el ements  . 

The  next  section  of  this  discussion  will  explain  the 
results  obtained  with  this  model,  and  will  compare  them  with  the 
solution  presented  by  Libove. 

6. 2. 2.1  Results  of  Model  with  Lid  Seal 

The  maximum  lid  deflections  and  stresses  obtained  with 
each  of  the  three  seal  materials  are  compared  to  the  value 
predicted  by  Libove  in  Table  28. 

As  shown  in  Table  28,  the  maximum  lid  deflections  as 
calculated  by  Libove  and  STARDYNE  correlate  fairly  well.  Note  that 
Libove's  solution  makes  no  allowance  for  variations  in  seal 
material,  and  he  predicts  a  single  lid  deflection  regardless  of 
seal  material.  However,  as  shown  by  the  STARDYNE  solution, 
Libove's  neglect  of  seal  material  effects  might  be  appropriate: 
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the  lid  deflection  variation  due  to  the  three  different  seal 
materials  analyzed  with  this  model  may  be  considered  insignificant. 

Also,  in  regard  to  the  maximum  lid  deflections,  the 
deflection  calculated  by  STARDYNE  for  each  of  the  three  seal 
materials  is  less  than  the  deflection  calculated  by  S1ARDYNE  for 
the  previous  "first-cut"  analysis.  This  is  not  what  was  expected, 
especially  for  the  glass  and  the  solder  cases.  The  addition  of  the 
seal  to  the  model  should  have  increased  the  maximum  lid  deflec¬ 
tions.  An  explanation  for  this  apparent  discrepancy  is  that 
different  elements  were  used  in  the  generation  of  the  two  models. 
Any  variation  in  the  stiffness  formulations  for  the  cube  elements 
and  the  cuadr i  1  ateral  plate  elements  could  result  in  different 
calculated  deflections  for  the  same  flat-pack  lid  (rememoer  that 
the  cuadr i 1 ateral  plate  elements  are  5  D.C.F.  per  node,  and  that 
the  cube  elements  are  3  D.C.F.  per  node). 

When  considering  the  maximum  lid  stresses  in  Table  ?8, 
it  appears  that  there  is  a  substantial  disagreement  between 
Libove's  value  and  the  STARDYNE  values.  However,  as  in  the 
"first-cut"  analysis,  the  stress  calculated  by  Libove  is  located  at 
the  midpoints  of  the  lid  edges;  the  stresses  computed  by  STARDYNE 
are  located  at  the  centroids  of  the  two  plate  elements  closest  to 
the  mid-side  nodes.  A  finer  element  grid  would  be  required  for 
STARDYNE  to  predict  a  stress  closer  to  the  mid-side  node  locations. 
For  this  reason,  the  lid  stress  comparison  of  Table  28  is  not 
completely  valid.  Yet,  it  is  interesting  to  note  that  the  effect 
of  varying  seal  materials  is  almost  insignificant  as  far  as  the  lid 
stress  is  concerned. 
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Tat 1 e  28  . 

Results  of  Model 

with  Lid  Sea1 

Seal  Materiel 

Parameter 

Kov  a  r 

Solder 

Glass 

Maximum 

Lid  Deflections: 

Libove 

0  .00426  " 

0 .00426" 

0  .00426 

ii 

STARDYNE 

0  .00  391  " 

0 .00299  " 

0  .00297 

" 

-•  di  ff  . 

-  8  . 2 r 

-6  .3". 

-6  .8% 

Maxi  mum 

Lid  Stresses: 

Libove 

34,533  psi 

34,533  psi 

34,533 

psi 

STARDYNE* 

13,943  psi 

13,790  psi 

13,829 

DS  1 

T  d  i  f  f  . 

-59  .61- 

-60  .1% 

-60  .0% 

M  a  x  i  m  urn 

Seal  Stresses: 

Libove  Vertical  Stress 

5,443  psi 

5,443  psi 

5,443 

psi 

STARDYNE : 

(a)  Von  Mi ses  Stress 

4,888  psi 

2,155  psi 

2,820 

psi 

(b)  Lateral  Stress 

-5,737  psi 

-2,662  psi 

-3,244 

psi 

(c)  Longitudinal  Stress 

-1,903  psi 

1,222  psi 

-828 

psi 

( d )  Verti cal  Stress 

-335  psi 

-  3  3  2  psi 

-333 

psi 

*  Von  Mi ses  Stress 

Finally,  Table  28  compares  the  maximum  seal  stresses  as 
predicted  by  Libove  and  by  STARDYNE.  For  the  STARDYNE  solution, 
four  different  stress  quantities  are  tab ul a  ted:  (1)  the  Von  Mi ses 
combined  stress  which  should  be  used  to  predict  seal  failure,  (2) 
the  lateral  normal  stress  across  the  width  of  the  seal,  (3)  the 
longitudinal  normal  stress  down  the  length  of  the  seal,  and  (4)  the 
vertical  normal  stress  across  the  thickness  of  the  seal.  Although 
the  STARDYNE  vertical  stresses  vary  insignificantly  with  seal 
material,  they  do  not  compare  well  with  Libove’s  calculated 
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vertical  tensile  stress  of  5,443  psi.  However,  once  again,  we  are 
not  able  to  compare  stresses  at  a  common  point.  The  STARDYNE 
stresses  are  all  calculated  at  the  centroid  of  the  seal  elements 
nearest  the  midpoints  of  the  lid's  edges,  while  Libove's  maximum 
stress  is  calculated  along  the  most  outward  edge  of  the  seal  at  tnt 
midpoints  of  each  of  the  lid's  edges.  It  is  very  possible  that  the 
seal  stress  distribution  across  the  width  of  the  seal  is  such  that 
the  vertical  stress  at  the  outside  edge  of  the  seal  is  about  5,000 
psi  (corresponding  to  Libove's  value)  and  the  vertical  stress  in 
the  center  of  the  seal  is  about  -300  psi  (corresponding  to 
STARDYNE's  value).  Figure  40  compares  Libove's  assumed  stress 
distribution  with  STARDYNE's  central  vertical  seal  stress;  note 
that  it  would  recuire  a  more  complicated  seal  model,  one  with 
several  elements  placed  across  the  seal  width,  to  determine  the 
actual  stress  distribution. 

Continuing  with  the  discussion  of  the  STARDYNE  seal 
stresses,  it  can  be  seen  that  the  lateral  and  longitudinal  stresses 
at  the  center  of  the  seal  are  larger  than  the  vertical  stresses. 
These  stresses  cannot  be  ignored,  and  therefore,  they  were  combined 
with  the  vertical  stress  to  compute  the  Von  Mises  combined  stress 
at  the  center  of  the  seal  (experience  has  shown  that  the  Von  Mises 
combined  stress  is  usually  a  better  method  of  predicting  structural 
failure  than  by  relying  on  only  a  single  component  of  the  overall 
stress  state).  Libove's  analysis  appears  to  suggest  that  the  only 
significant  stress  at  the  outside  edge  of  the  seal  is  the  vertical 
tensile  stress;  therefore,  the  Von  Mises  combined  stress  for 
Libove's  analytical  stress  state  is  probably  very  close  to  the 
value  indicated  in  the  table  for  his  vertical  stress  component. 
Again,  however,  the  STARDYNE  Von  Mises  stress  should  probably  not 
be  expected  to  compare  with  the  Libove  stress  since  they  ere  not 
calculated  at  the  same  location  in  the  seal. 


LIBOVE'S  STRESS  MODEL: 
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The  next  and  final  analysis  which  was  performed  with 
regard  to  Critical  Problem  *1  concerned  a  model  of  a  single  section 
of  the  lid  seal.  The  purpose  of  this  final  analysis  was  to  obtain 
some  idea  of  the  vertical  stress  distribution  across  the  width  of 
the  seal  . 

6 .2 .3  Local  Seal  Model 

In  order  to  more  accurately  determine  the  vertical 
stress  distribution  across  the  seal  width,  we  decided  to  take  a 
closer  look  at  a  single  section  of  the  lid  seal.  Obviously,  from 
an  analytical  viewpoint,  it  would  have  been  desirable  to  simply 
modify  our  previous  work  by  incorporating  a  more  detailed  seal 
model  into  the  finite  element  model  of  the  flat-pack.  However,  the 
seal  had  to  be  modeled  with  several  elements  located  across  the 
seal  width  in  order  to  more  accurately  predict  the  stresses  along 
the  interior  and  exterior  wall  surfaces,  and  also  to  give  a  more 
meaningful  distribution  of  stress  values  across  the  width  of  the 
seal.  To  do  so  for  the  entire  flat-pack  model  would  have  resulted 
in  a  complex  and  relatively  expensive  finite  element  model.  There¬ 
fore,  it  was  decided  to  take  a  quick  and  easy  look  at  the  problem 
by  analyzing  only  a  small  portion  of  the  total  seal. 

In  effect,  a  single  “strip"  was  removed  from  the  most 
highly  stressed  seal  area  of  our  previous  "  ouarter-model  " ;  this 
strip  of  elements  included  lid  elements,  seal  elements,  and  wall 
elements.  It  was  then  attempted,  through  several  methods,  to  force 
this  strip  of  the  original  lid  into  the  deflection  shape  which  this 
particular  group  of  lid  elements  had  assumed  during  our  prior 
analysis.  By  doing  so,  we  hope  to  apply  the  same  forcing  function 
as  originally  seen  by  our  chosen  seal  elements  to  the  more  detailed 
model  of  this  section  of  the  seal.  Figure  41  describes  the  seal 
model  and  shows  its  relationship  to  the  previous  flat-pack  model. 
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As  shown  in  Figure  41,  the  seal  was  modeled  using  three 
rows  of  cube  elements  across  its  width.  To  better  determine  the 
vertical  stress  distribution  across  the  seal's  width,  four  thin 
(.0001")  membrane  plate  elements  were  added  to  the  vertical  faces 
of  the  three  seal  cube  elements  .  These  elements  were  located  at 
the  inner  and  outer  walls  of  the  seal  and  at  the  two  adjoining 
vertical  faces  of  the  interior  seal  cube  elements.  This  method,  in 
effect,  allowed  us  to  use  seven  elements  to  determine  the  vertical 
stress  distribution  across  the  seal's  width.  To  provide  the 
reouired  element  connectivity,  it  was  necessary  to  also  subdivide 
the  lid  and  wall  elements  in  the  immediate  vicinity  of  the  seal. 
With  the  above  modifications,  our  local  seal  model  consisted  of  48 
nodes,  13  cube  elements,  and  14  quadrilateral  plate  elements  (10  of 
these  plate  elements  were  added  to  the  exterior  surfaces  of  the 
lid,  seal,  and  wall  elements  so  that  pressure  could  be  applied  to 
our  cube-element  structure).  Also,  it  should  be  mentioned  that 
Kovar  was  considered  to  be  the  seal  material  for  this  analysis. 

To  recreate  the  state  of  stress  present  in  the  seal  of 
our  previous  model,  it  was  necessary  to  recreate  the  overall 
forcing  function  which  caused  the  original  stress  state.  Several 
potential  methods  were  considered,  but  time  did  not  allow  all 
methods  to  be  investigated.  Therefore,  two  loading  conditions  were 
considered  for  this  dislocated  portion  of  the  flat-pack  model. 

Both  of  the  loading  conditions  which  we  considered  were 
based  on  the  assumption  that  the  lid's  center  deflection  was  the 
prime  influence  on  the  stress  state  of  the  seal.  The  first  loading 
condition,  therefore,  was  to  simply  force  the  node  representing  the 
center  of  the  flat-pack  (node  41)  to  deflect  the  same  amount  as 
discovered  in  our  previous  lid  analysis  (0.00391").  We  felt  that 
the  rema i nder  of  the  nodes  would  then  deflect  in  a  manner  similar 
to  the  deflection  shape  previously  determined,  the  majority  of  any 
difference  probably  being  due  to  the  different  stiffness  formula¬ 
tion  of  the  seal  elements  with  the  new  gridwork  modification  in 
effect . 
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The  second  loading  condition  was  identical  to  the  -first 
except  that  the  30  psi  pressure  loading  was  applied  to  the  top 

surface  of  the  lid  and  to  the  exterior  surfaces  of  the  wall  and  the 
side.  It  was  hoped  that  the  addition  of  this  pressure  would  help 
recreate  the  original  deflection  shape  even  more  accurately  than 
for  the  first  loading  case. 

Another  possible,  seemingly  obvious,  loading  case 
should  be  discussed.  It  would  be  possible  to  force  every  node 
common  to  both  this  model  and  our  original  model  to  deflect  the 
amount  that  was  calculated  in  the  original  analysis  of  the  model 

with  the  seal.  This  would  seem  to  be  an  obvious  solution  method, 
except  for  one  consideration.  The  previously  calculated  nodal 
deflections  were  based  on  the  seal  stiffness  formulation  of  the 

original  model.  However,  by  changing  the  gridwork  of  the  seal 

elements,  the  effective  seal  stiffness  has  been  modified  somewhat. 
Therefore,  in  a  region  close  to  the  revised  seal  elements,  the 

previously  calc ul  a ted  deflections  might  not  truly  represent  the 

deflection  shape  which  the  lid  would  assume  with  the  modified  seal 

elements.  In  effect,  we  would  be  forcing  a  deflection  shape  which 

might  not  ever  be  realized  had  the  30  psi  pressure  been  applied  to 
a  complete  lid  model  with  our  modified  seal  elements.  This 
explains  the  reason  why  our  forcing  loading  conditions  were 
concerned  with  the  central  deflection  only:  we  assumed  that  the 

modified  seal  elements  would  not  have  significantly  affected  the 
central  deflection,  but  that  their  different  stiffness  formulation 
would  have  only  affected  nodes  in  the  local  region.  Therefore,  we 
forced  the  central  deflection  to  be  the  same  as  for  our  previous 
analysis,  but  we  allowed  the  deflections  of  nodes  near  the  seal  to 
assume  a  deflection  consistent  with  the  revised  seal  element 
stiffness . 
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6  .2  .3  .1 


Results  of  Local  Seal  Model 

The  stresses  which  were  calculated  for  the  elements 
across  the  width  of  the  seal  are  summarized  in  Table  29  for  both 
loading  cases.  As  shown  in  the  table,  the  vertical  seal  stresses 
vary  si gni f icantly  across  the  width  of  the  seal. 


Table  29.  Seal  Stress  Across  Seal  Width 

Seal  Stress  Location  Seal  Stress  ( ps  i  ) 

Von 


Loading  Description 

E 1 ement 

Location 

Vert . 

Long  . 

L at . 

M  i  s  e  s 

Forced  Center 

Plate  24 

Inside 

-4,488 

-1,170 

- 

4,191 

Di spl acement 

Cube  11 

-3,961 

-2  ,486 

-1,833 

4  ,471 

(load  case  1 ) 

Plate  25 

2,912 

898 

- 

2,591 

Cube  12 

Middle 

445 

1  ,612 

605 

2,155 

Plate  26 

-69 

-69 

- 

192 

Cube  13 

191 

633 

170 

1  ,273 

Plate  27 

Outs i de 

1  ,171 

245 

- 

1  ,095 

Forced  Center 

Plate  24 

Inside 

-8,207 

-2,125 

- 

7  ,462 

Di spl acement 

Cube  11 

-6,788 

-4,567 

-3,173 

7,778 

PI  us  Press  ure 

Plate  25 

4,890 

1,598 

- 

4,329 

(load  case  2) 

Cube  12 

Middle 

749 

2,680 

1,094 

3,600 

Plate  26 

-172 

-52 

- 

310 

Cube  13 

287 

1  ,049 

349 

2,159 

Plate  27 

Outs i de 

2  ,025 

503 

- 

1  ,849 

Fig u re  42  shows  the  vertical  stress  distribution  across 
the  seal  width.  Along  the  inside  edge  of  the  seal  there  are  large 
compressive  stresses;  along  the  outside  edge  of  the  seal,  the 
stresses  are  smaller  in  magnitude,  and  they  are  tensile  instead  of 
compressive.  For  both  load  cases,  it  is  obvious  that  the  vertical 
stress  distribution  is  not  linear.  To  better  understand  this 
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Figure  42.  Vertical  Stress  across  Seal  Width 
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phenomena,  one  must  consider  the  differences  between  Libove's 
closed-form  solution  and  the  finite  element  solution.  First, 
Libove  assumes  a  stress  distribution  as  shown  on  Page  62  of  his 
report  "881  .  Examination  of  this  distribution  and  Libove's 
formulas  for  the  maximum  seal  stress  appears  to  indicate  that  he 
assumes  ecual  maximum  tensile  and  compressive  stresses.  However, 
the  reactive  forces  applied  to  the  seal  appear  to  be  concentrated 
towards  the  inner  wall,  indicating  that  these  maximum  stresses 
should  really  not  be  ecual .  The  force  distribution  seems  to  imply 
that  the  inner  compressive  stress  should  be  greater  than  the  outer 
tensile  stresses. 

Second,  the  finite  element  solution  considered  the 
flexibility  of  the  wall  underneath  the  seal.  Libove's  solution 
appears  to  have  ignored  this  effect.  This  wall  flexibility  may 
have  contributed  somewhat  to  the  seal  stress  distribution 
determined  by  our  finite  element  model.  The  distribution  indicates 
that  the  effect  of  the  applied  moment  is  dominant,  and  is  reacted 
by  the  inner  half  of  the  seal.  Note  that  the  maximum  tensile  and 
compressive  stresses  in  the  inner  half  of  the  seal  width  are  on  the 
same  order  of  magnitude,  but  the  compressive  stresses  are  larger. 
Also,  note  that  there  is  a  tensile  lifting  of  the  seal  at  the  outer 
edge  of  the  lid.  Figure  42  also  shows  that  STARDYNE's  calculated 
stress  distribution  varies  depending  on  the  load  case.  It  is  not 
immediately  obvious  which  load  case  results  in  the  more  accurate 
stresses,  although  it  seems  reasonable  that  load  case  #2  might  be 
the  better  of  the  two. 

In  summary,  the  maximum  seal  stress  magnitudes  as 
predicted  by  Libove's  closed-form  solution  and  our  finite  element 
models  are  in  the  same  "ballpark."  However,  the  vertical  stress 
distributions  predicted  by  these  two  methods  are  markedly 
di  f f erent . 

Basically,  what  we  have  shown  with  this  local  seal 
analysis  is  that  the  vertical  seal  stress  distribution  is  not 
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likely  to  be  linear.  Howevc  ,  this  has  only  been  an  initial  look 
at  the  problem;  more  analysis  would  be  reouired  to  accurately 
determine  the  seal  stress  distribution. 

The  next  section  describes  the  second  critical  problem. 

6  .3  Critical  Probl  em  #2 

Critical  Problem  *2  was  concerned  with  a  hybrid  package 
consisting  of  an  Isotronics  base  and  a  Device  Closures  lid.  The 
dimensions  of  the  base  and  the  lid  are  shown  in  Figure  43. 

In  the  following  two  analyses,  the  hybrid  package  will 
be  analyzed  to  determine  its  response  to  a  pressure  load  of  30  psi, 
a  series  of  concentrated  loads  of  3  lb,  and  a  thermal  load  induced 
by  a  180°  F  temperature  differential.  The  first  analysis  will  be  a 
relatively  simple  analysis  of  the  hybrid  lid;  the  walls  will  be 
neglected  for  this  first  analysis.  Also,  this  analysis  will  take 
advantage  of  the  symmetry  of  the  lid  and  its  loads,  and  we  will, 
therefore,  use  a  "cuarter"  finite  element  model. 

The  second  analysis  of  this  package  will  use  a  "full" 
finite  element  model  of  the  hybrid's  base,  walls,  and  lid.  A  full 
model  was  used  so  that  we  could  determine  deflections  caused  by 
asymmetrical  point  loads.  These  calculated  deflections  were  then 
compared  to  the  actual  measured  deflections  induced  by  the  various 
point  loads  which  were  applied  to  the  hybrid  lid.  In  addition,  the 
analytical  thermal  strains  were  compared  to  the  thermal  strains 
measured  in  the  laboratory. 

6.3.1  "First-Cut"  Finite  Element  Model 

To  begin  the  analysis  of  Critical  Problem  *2,  we 
decided  to  perform  a  relatively  simple  analysis  of  the  hybrid  lid 
neglecting  the  effects  of  the  walls  and  the  hybrid  base.  Due  to 
our  use  of  a  "cuarter"  model  for  this  analysis,  we  were  limited  to 
only  investigating  symmetrical  load  cases. 

The  model  used  for  this  analysis  is  shown  in  Figure  44; 
as  shown,  it  consists  of  161  nodes,  72  cube  elements,  and  72 
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ISOMETRIC  VIEW 


elements,  and  the  model  included  the  s tepped- th i ck ness  feature  of 
the  actual  lid.  Also,  the  quadri 1 ateral  plate  elements  (0.0001" 
thick)  were  applied  to  the  top  and  bottom  surfaces  of  the  lid  to 
incorporate  the  effect  of  the  nickel  plating  used  on  the  real  lid. 
In  addition,  these  elements  allowed  us  to  apply  an  external 
pressure  to  the  upper  surface  of  the  lid  (remember  that  STARDYNE's 
cube  elements  may  not  be  directly  loaded  with  a  pressure). 

The  thin  outer  perimeter  of  the  lid  was  assumed  to  have 
a  thickness  of  0.005"  (average  of  the  dimensioned  values),  and  the 
central  portion  of  the  lid  was  assumed  to  be  0.015"  thick.  The  lid 
model  was  fully  restrained  at  the  nodes  along  the  bottom  surface  of 
the  lid  peri  meter;  in  a  more  complete  model  these  nodes  would  have 
connected  the  lid  to  the  hybrid  walls.  Also,  the  nodes  along  the 
axes  of  symmetry  were  restrained  consistent  with  the  necessary 
boundary  conditions  (note  that  it  was  not  necessary  to  restrain  the 
various  rotations  since  this  is  a  3-DOF  cube  element  model,  and  all 
nodal  rotations  are  identically  zero). 

It  should  be  mentioned  that  the  above  restraints  apply 
to  the  pressure  loading  case  and  the  point  loading  case  only.  For 
the  thermal  strain  analysis,  all  boundary  restraints  were  removed, 
one  node  was  full  restrained,  and  the  lid  was  allowed  to  expand 
freely  (the  base  is  also  Kovar,  and  will  expand  very  nearly  the 
same  as  the  lid). 

6. 3. 1.1  "First-Cut"  Results 

The  maximum  deflections  and  stresses  for  both  of  the 
force  load  cases  are  summarized  in  Table  30  below. 

Table  30.  "First-Cut"  Results  for  Critical  Problem  #2 
Load  Case  Center  Deflection  L i d  Stress  ( el ement ) 

30  psi  Pressure  0.02325"  28,700  psi  (cube  32) 

3  lb  Central  0.00513"  6,647  psi  (cube  72) 
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As  shown  in  Table  30,  the  30  psi  external  pressure 
produces  a  central  lid  deflection  of  0.02325",  and  a  maximum  Von 
Mises  combined  stress  of  28,700  psi  in  cube  32  (cube  32  is  located 
adjacent  to  the  midpoint  of  the  long  side  of  the  hybrid,  and  it  is 
one  of  the  0.005"  thick  perimeter  elements,  this  stress  is  less 
than  the  yield  strength  of  50,000  psi). 

The  equations  presented  in  Roark's  Formul as  for  Stress 
and  Strain  [85]  may  be  used  to  show  that  a  0.015"  thick,  clamped 
rectangular  Kovar  plate  (1.66"  x  1.06")  would  have  a  central 
deflection  of  0.01386"  under  a  30  psi  load.  This  is  less  than  our 
calculated  lid  deflection  because  Roark  does  not  consider  the 
effect  of  the  variable  thickness  lid  (refer  to  Appendix  A). 

Also,  the  calculated  central  lid  deflection  of  0.02325" 
due  to  the  pressure  loading  falls  into  the  range  of  deflections 
calc ul a  ted  by  libove;  Libove  suggests  that  the  deflection  will  fall 
between  0.0128"  and  0.0402"  depending  upon  the  extent  of  plastic 
flow  which  is  initiated  in  the  lid  (refer  to  Appendix  B,  and  Libove 
pp  .  23-26  )  . 

Table  30  also  shows  that  the  center  of  the  lid  will 
deflect  approximately  0.005"  when  subjected  to  a  3  1 b  central  load. 
The  correspond!'  ng  maximum  Von  Mises  stress  occurs  in  the  center  of 
lid,  and  has  a  value  of  about  6,650  psi.  Note  that  these  stress 

values  were  used  in  determining  maximum  loads  for  use  in  our 

laboratory  investigations.  (This  lid  deflection  may  be  compared 
with  actual  measured  deflection  discussed  later  in  this  report). 

The  final  portion  of  this  analysis  consisted  of  the 
180°F  differential  temperature  loading.  The  lid  was  assumed  to  be 
at  an  initial  temperature  of  77°F,  and  the  thermal  strains  were 
calculated  for  the  lid  after  it  was  raised  to  a  temperature  of 
257°  F  .  A  coefficient  of  thermal  expansion  ( CTE )  of  3.05  x  10‘6 

in./in.°F  was  used  for  the  Kovar  lid,  and  a  CTE  of  8.40  x  10'6 

in./in.°F  was  used  for  the  thin  nickel  plating.  The  resulting 
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in-plane  thermal  strain  was  calculated  to  be  5.70  x  10  in/in. 
Again,  this  strain  will  later  be  compared  to  experimental  data. 

6  .3  .2  Second  Finite  Element  Model 

To  continue  our  analysis  of  critical  problem  *2,  we 

decided  to  model  the  hybrid  with  a  somewhat  more  detailed  finite 

element  model  which  would  allow  us  to  calculate  deflections  for 

asymmetrical  loads.  This  model  included  elements  which  represented 
the  hybrid  lid,  walls,  and  base,  and  would  be  used  to  compare 
calculated  deflections  and  thermal  strains  to  laboratory  values. 

The  finite  element  model  is  shown  in  Figure  45.  As 

shown,  the  model  was  constructed  entirely  of  plate  elements  (122 
nodes,  116  cuad-i 1 ateral s ,  and  4  tri-plates).  For  the  pressure  and 
point  force  loads,  all  of  the  nodes  along  the  bottom  of  the  hybrid 
were  fully  restrained.  However,  for  the  thermal  load,  only  a 
single  corner  node  was  restrained  so  that  the  hybrid  would  be  free 
to  expand  as  dictated  by  its  CTE . 


6  .3  .2  .1 


Pressure  Loading  Results 


The  following  Table  31  summarizes  the  results  of  the 
pressure  loading  case. 

Table  31.  Pressure  Loading  of  Hybrid  Model 


Load 

Case 

30  psi 
Pressure 


F.E.M  Libove 

Lid  Deflection  Lid  Deflection 


Lid  Stress 
0.015"  0.005" 

Central  Edge 

Section  Section 


0  .03240" 


0.0138"  0.0402"  43,900  psi  210,300  psi 

(Quad  96)  (Quad  71) 


As  shown  in  Table  31,  the  calculated  lid  deflection  is 
approximately  0.0324"  which  is  greater  than  our  previous  calculated 
value  of  0.0233".  The  difference  is  due  to  the  different  stiffness 
formulations  of  the  cube  elements  from  the  previous  model  and  the 
plate  elements  of  the  present  model. 


Also,  from  the  table,  it  appears  that  the  maximum 
stress  in  the  thick  (0.015")  portion  of  the  lid  is  reasonable,  and 
this  stress  occurs  in  the  elements  adjacent  to  the  center  of  the 
lid.  However,  the  maximum  stress  calculated  for  the  thin  perimeter 
of  the  lid  (0.005")  is  not  reasonable.  This  stress  occurs  at  the 
midpoints  of  the  long  edges  of  the  lid,  but  its  magnitude  is 
adversely  affected  by  our  use  of  excessively  large  aspect  ratio 
quadrilateral  plates  in  this  region.  (Ideally,  aspect  ratios  for 
plate  would  be  kept  between  1  and  3  for  the  most  accurate  results.) 
If  we  desired  to  obtain  a  better  stress  in  this  area,  a  different 
gridwork  would  be  reouired  for  our  model. 

It  should  be  noted  that  no  laboratory  results  were 
available  to  correlate  with  our  calculated  deflections  for  the 
pressure  loading  case.  However,  for  both  the  point  load  cases  and 
the  thermal  load  case,  experimental  data  was  collected  and  used  for 
correlation  purposes. 

6  .3  .2  .2  Point  Loading  Results 

As  shown  in  Figure  46,  five  different  point  load  cases 
were  considered  in  the  laboratory;  in  each  case  a  load  of  3  lb  was 
applied  to  the  lid.  Figure  46  shows  the  points  of  application  of 
these  loads.  Appendix  C  presents  the  test  plan  which  was  followed 
during  the  laboratory  investigations. 

Table  32  indicates  both  our  measured  deflections  and 
our  calculated  deflections  for  each  point  of  load  application. 
Note  that  there  are  two  sets  of  calculated  deflections.  The  first 
corresponds  to  the  deflections  based  on  an  assumed  central  lid 
thickness  of  0.015"  and  an  assumed  edge  thickness  of  0.005". 
However,  after  the  deflection  tests  were  complete,  the  hybrid  lid 
was  sectioned,  and  actual  thickness  measurements  were  made.  It  was 
found  that  the  center  portion  of  the  lid  varied  between  0.0140" 
0.0145".  Therefore,  the  second  set  of  calculated  deflections 
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Figure  46.  Points  of  Concentrated  Load  Application 
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represents  the  deflections  after  the  above  thickness  changes  were 
included  in  the  lid  model. 


Point  of 

A  ppl i cati on 

Table  32  . 

Meas  ured 
Deflections 

Point  Loading 

Calculated 
Nominal  Thick 
for  Lid 

of  Hybrid  Model 

Deflections 
Actual  Thick 
for  Lid 

Of 

to 

Difference* 

A 

0 .00575" 

0  .00551" 

.00614" 

+6  .8% 

B 

0 .00625" 

0  .00551  " 

.00624" 

-0  .1% 

C 

0  .00456" 

0  .00444" 

.00493" 

+8.1% 

D 

0  .00775" 

0  .00693  " 

.00785" 

+  1  .3% 

E 

0  .00562 

0  .00444" 

.00518 

-7  .8% 

*Percent  difference  between  our  "best"  calculated  deflections  and 
the  measured  deflections 

As  shown  in  Table  32,  the  measured  lid  deflections 
correlate  very  well  (within  8%)  with  the  calculated  lid  deflections 
which  were  based  on  the  actual  lid  thickness.  The  error  would  have 
been  significantly  greater  had  we  not  reiterated  our  calculated  lid 
deflections  once  the  actual  lid  thicknesses  were  known .  This  is 
just  one  example  of  how  the  analyst  may  use  a  combination  of 
analysis  and  test  data  to  more  fully  understand  what  may  seem 
initially  to  be  significant  deviations  between  his  analysis  and  the 
test  res  ul ts  . 

6  .3  .2  .3  Thermal  Loading  Results 

As  previously  done  in  our  "first-cut"  analysis,  we 
thermally  loaded  the  hybrid  with  an  increase  in  temperature  of 
180°F,  and  calculated  the  resulting  thermal  strains.  Using  a  CTE 
of  3.05  x  10  6  in./in.°F  for  the  hybrid,  we  determined  in-plane 
strains  of  5.49  x  10"4  in/in.  However,  it  should  be  noted  that 
this  calculation  did  not  incl ude  the  effects  of  the  nickel  plating 
on  the  thermal  expansion  of  the  lid.  Therefore,  assuming  a  plating 
thickness  of  0.0001",  an  effective  CTE  of  3.16  in/in/°F  may  be 


found  from  our  previous  analysis  which  results  in  a  thermal  strain 

_4 

of  5.69  x  10  in/in  for  this  model. 

To  experimentally  determine  the  lid  strain  due  to  a 
180°F  temperature  change,  it  was  decided  to  instrument  the  hybrid 
lid  with  a  pair  of  strain  gages.  Unfortunately,  we  did  not  have 
available  the  required  gages  for  this  experiment  ( i  .e . ,  no  gages  of 
the  proper  size  and  with  a  CTE  of  3  x  10  ^  in./in.°F).  However,  we 
did  decide  to  go  ahead  with  the  test  using  the  gage  that  we  had 
available,  realizing  that  we  were  introducing  some  error  into  our 
meas  urements  . 

Basically,  we  were  required  to  adjust  our  measured 
strain  data  to  consider  the  following  four  aspects: 

1.  Mismatched  CTE  between  gage  and  lid 

2.  Apparent  strain  due  to  tempe r a t ure / r e s i s t i v i ty 
effects  (higher  temperature  increases  resistance  of 
gage) 

3.  Gage  factor  variation  with  large  temperature 
di f f erenti al s 

4.  Transverse  strain  effects  on  each  gage 

At  the  conclusion  of  our  strain  gage  testing,  and  after 
correcting  the  data  for  the  above  effects,  it  was  obvious  that  our 
data  was  not  sufficiently  accurate  to  correlate  with  our  analysis. 
Three  factors  had  combi ned  to  invalidate  our  results:  (1)  the 

relatively  large  temperature  differential  exaggerated  the  effect  of 
the  CTE  mismatch  between  the  gage  and  the  lid,  (2)  the  relatively 
high  test  temperatures  produced  apparent  strains  of  the  same  order 
of  magnitude  as  the  strains  we  were  trying  to  mr=»sure,  and  (3)  we 
relied  on  the  published  average  "correcti  on-data'1  for  this 
particular  type  of  gage;  for  better  accuracy,  the  gage  in  use  would 
need  to  be  tested  to  obtain  its  individual  correction  factors. 
(Time  constraints  did  not  allow  us  to  continue  with  this  effort.) 


Since  our  strain  gage  data  did  not  prove  to  be  us  able, 
we  dedicated  to  use  a  dilatometer  to  measure  the  thermal  strain  in 
the  hybrid  lid  (refer  to  Figure  47  for  a  sketch  of  a  dilatometer). 
The  following  data  was  collected  and  is  presented  as  Table  33. 


Table  33  . 

Di 1  a  tome  ter 

Strain  Measurements 

(zero  strain  0  74°F) 

T 

Dial  Gage 

Avg  .  Strain 

Effective  CTE 

183°  F 

+0  .0335  mm 

7  .60  x  10'4 

4.15  x  10'6/° 

93°F 

+0  .0180  mm 

4  .08  x  10"4 

4  .39  x  1  O' 6 / 0 

-42  0  F 

-0  .0080  mm 

-1  .82  x  10"4 

4.32  x  10_6/° 

- 142°  F 

-0  .0290  mm 

-6  .58  x  10-4 

4  .63  x  10‘6/° 

Average  CTE  = 

4.37  x  10'6/° 

As  shown  in  Table  33,  the  average  thermal  strain  due  to 
a  temperature  differential  of  180°F  is  measured  to  be  about  7.6  x 
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10  in/in.  This  does  not  correlate  well  with  our  calculated 

-4 

strain  of  5.69  x  10  in/in. 

Two  possible  explanations  of  this  discrepancy  are 
apparent.  First,  there  may  be  more  than  a  0.0001"  thick  layer  of 
nickel  plating  on  the  Kovar  lid  (0.0001"  was  only  an  assumed 
value).  If  this  is  true,  our  effective  CTE  would  increase  slightly 
from  the  calculated  value  of  3.16  x  10’6  in/in.  °F,  and  the 

calculated  thermal  strain  would  correspondingly  increase.  (Note 
that  the  measured  effective  CTE  of  this  hybrid  is  4.37  x  10"6 

in/in°F).  Figure  48  plots  the  variation  of  thermal  strain  with 

temperature,  and  shows  the  results  of  both  our  test  seouence,  and 
some  Westinghouse  test  data  [92].  Note  that  the  slope  of  the 

curves  indicates  the  CTE  for  Kovar.  Although  there  is  some 

deviation  between  the  two  sets  of  data,  it  is  not  excessive. 

Second,  we  are  relying  on  statistical  "average"  data 

for  our  values  of  the  CTE  for  Kovar.  This  does  not  consider  the 
variation  in  properties  from  lot  to  lot  of  the  material  ( 1 -e . .  it 

is  possible  that  our  hybrid  is  made  of  Kovar  which  may  be  several 
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standard  deviations  removed  from  the  "average"  Kovar  sample). 

Also,  the  Kovar  CTE  used  in  our  analysis  does  not  include  any 

effects  which  may  be  caused  by  the  different  working  processes  our 

sample  has  seen  during  the  production  of  the  hybrid. 

In  summary  of  the  thermal  strain  analysis,  it  is 

obvious  that  our  measured  and  calculated  strains  did  not  correlate 

well.  However,  this  appears  to  be  due  to  an  error  in  assuming  an 

effective  CTE  for  the  hybrid  model  (had  we  used  a  CTE  of  4.37  x 

1 0  "  6  in/in°F  in  our  analysis,  our  calculated  strain  for  a 

D  —  4 

temperature  differential  of  180  F  would  have  been  7.87  x  10  in/in 
--  much  closer  to  the  measured  strain). 

It  should  be  mentioned  that  two  control  tests  were  run 
with  the  dilatometer.  First,  the  empty  dilatometer  was  cycled  over 
temperature  extremes  to  measure  its  strain  over  this  temperature 
range  --  its  measured  strain  was  "zero".  This  means  that  the 

dilatometer  was  not  exaggerating  the  strain  of  any  sample  being 
tested.  Second,  a  sample  of  6061-T651  aluminum  was  tested;  its 

measured  CTE  was  approximately  7.5%  1  ower  J:han  published  values. 
Therefore,  we  concluded  that  the  dilatometer  was  not  producing 
consistently  higher  values  of  CTE  than  published. 

The  next  section  summarizes  the  results  of  the  critical 
problem  anal yses  . 

6  .4  General  Summary  of  Critical  Problem  Analysis  and 

Correlation  to  Test  Data  (Phase  3) 

The  following  is  a  summary  of  the  findings  of  the  two 
critical  problem  analyses. 


6.4.1  Critical  Problem  #1 

1.  Good  correlation  (3.7%)  was  found  between  the  lid 
deflections  calculated  by  Libove  and  our  "first- 
cut"  finite  element  model.  This  model  was 
purposely  kept  simple;  it  required  approximately 
one  hour  to  generate  the  model  and  to  obtain  the 


results.  However,  due  to  the  simplicity  of  this 
model,  no  elements  were  included  to  represent  the 
lid  seal,  and  therefore,  no  seal  stresses  were 

f 

obtained  with  this  model. 

2.  Good  correlation  was  also  found  between  the  lid 

deflections  calculated  by  Libove  and  those 
calculated  by  our  second,  more  detailed  finite 
element  model.  This  included  a  0.005"  thick  lid 

seal,  and  three  different  analyses  were  performed 
based  on  three  possible  seal  materials.  It  should 
be  noted  that  the  calculated  lid  deflections  and 
stresses  varied  only  slightly  with  the  seal  material 

3.  The  Von  Mises  seal  stresses  calculated  by  our 

second  finite  element  model  vary  significantly 
depending  upon  the  seal  material.  However,  it 
should  be  noted  that  the  vertical  stress  component 
at  the  center  of  the  seal  did  not  vary  much  with 
seal  material,  and  this  stress  is  low  (as  predicted 
by  Libove's  vertical  stress  distribution). 

4.  A  relatively  simple  finite  element  model  of  the 

most  highly  stressed  seal  area  was  generated  to 
better  determine  seal  stresses  and  the  seal  stress 
distribution  across  the  seal  width.  Results 
confirmed  that  the  seal  experienced  large 
compressive  stresses  nearest  the  package  cavity; 
however,  these  preliminary  results  also  indicated 
that  the  seal  stress  distribution  is  most  likely 
not  linear  across  the  seal  width.  A  more  detailed 
finite  element  analysis  would  be  required  to 
investigate  this  subject  in  greater  depth.  It  is 
recommended  that  an  investigation  of  the  seal 
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stress  distribution  be  included  in  the  work  of  a 
f uture  study  . 

Critical  Problem  #2 

1)  A  simple  "cuarter"  finite  element  model  of  the  lid 
may  be  used  to  approximate  lid  deflections  and 
stresses  for  symmetrical  loads.  This  model 
confirmed  that  a  stepped  lid  will  experience 
greater  deflections  for  any  given  load  than  would  a 
constant  thickness  lid.  The  maximum  lid  deflection 
calculated  by  this  model  fell  inside  the  deflection 
range  predicted  by  Libove. 

2)  A  more  complex  finite  element  model  was  generated 
in  order  to  handle  asymmetric  load  cases.  Calcu¬ 
lated  deflections  for  various  point  load  cases 
correlated  well  (within  8")  with  experimental  data 
once  the  actual  lid  thickness  had  been  measured. 
For  these  calculations,  it  was  not  sufficiently 
accurate  to  rely  on  the  nominal  lid  dimensions. 

3)  It  is  recommended  that  future  design  analyses 
consider  lid  thickness  tolerances  --  minimum  lid 
thickness  values  should  be  used  in  design  analyses. 

4)  The  tht'mal  strains  calculated  by  the  two  finite 
element  models  were  virtually  identical,  and  were 
constant  at  all  points  on  the  hybrid  lid.  The 
predicted  thermal  strains  occurring  over  a  wide 
temperature  range  did  correlate  well  with  our 
measured  strain  values.  However,  CTE  values  did 
not  correlate  well--this  may  be  caused  by  a 
difference  between  the  published  value  and  the  real 
value  of  the  effective  CTE  for  the  hybrid. 


7.0  ANALYSIS  GUIDELINES  (PHASE  4) 

7  .1  Introducti on 

"The  whole  concept  of  the  finite  element  technique 

can  be  compared  with  that  of  a  jigsaw  puzzle.  A  jigsaw  is 
usually  a  complete  picture  which  is  broken  down  into  small 
irregular-shaped  pieces.  The  problem  is  then  to  reconstruct  the 
picture  by  assembling  the  individual  parts.  Obviously,  the 

pieces  cannot  be  assembled  arbitrarily,  but  have  to  be  assembled 
according  to  certain  rules. "[89]  The  above  quote  is  an  accurate 
description  of  finite  el ement  modeling.  It  also  infers  that 
finite  element  modeling  can  sometimes  become  more  of  an  "art" 
than  a  "science "  . 

The  purpose  of  this  guidelines  section  is  to  provide 
the  mechanical  engineer  with  general  procedures  which  can  be 

applied  to  microelectronics  packaging  analysis.  Simply  stated, 
the  finite  element  process  can  be  divided  into  three  tasks:  (1) 
understanding  the  problem,  (2)  developing  the  model,  and  (3) 
interpreting  the  results.  These  three  tasks  will  now  be 

discussed  in  detail. 

7 .2  Understanding  the  Problem 

This  step  sets  the  tone  for  the  entire  analysis,  and 

is  oftentimes  not  considered  in  sufficient  detail.  It  is 

imperative  that  the  analyst  understand  the  problem  from  both  a 
physical  as  well  as  an  analytical  standpoint. 

A  complete  understanding  of  the  problem  requires  the 
analyst  to:  (1)  know  what's  required,  (2)  consider  program 
choices,  (3)  perform  simple  hand  calculations,  and  (4)  recognize 
typical  microelectronics  problems. 

a.  know  What's  Required 

When  an  analyst  is  assigned  to  a  problem,  he 
needs  to  understand  the  entire  scope  of  the 


problem.  Parameters  that  influence  the 
problem's  scope  are  those  engineering  and 
project  constraints  that  are  imposed  at  the 
outset.  For  microelectronics  packaging,  typical 
engineering  constraints  include  environmental 
definition,  system/subsystem  relationship, 
device  physical  description,  and  specific  design 
concerns  based  on  past  experience.  Project 
constraints  include  cost,  schedule,  workload, 
availability  of  analysts,  etc. 

Before  the  analyst  proceeds,  he  must  answer 
the  following  ouestions:  (1)  what  am  I 
analyzing,  (2)  why  am  I  analyzing  it,  (3)  who  am 
I  analyzing  it  for,  and  (4)  when  am  I  expected 
to  produce  results?  When  the  above  questions 
are  answered,  the  analyst  can  proceed  in 
selecting  a  suitable  finite  element  program. 

b.  Consider  Program  Choices 

The  topic  of  program  selection  has  been 
previously  discussed  in  Section  4.0  of  this 
report  -  Analysis  Methods.  For  microelectronics 
analysis,  the  choice  of  programs  is  usually 
dictated  by  the  imposed  environmental 
requirements.  For  instance,  if  random  vibration 
were  considered  the  most  important  environment, 
the  analyst  would  not  choose  a  program  such  as 
ANSYS,  since  it  cannot  provide  random  vibration 
analysis.  A  choice  such  as  STARDYNE  or  NISA 
would  be  more  appropriate. 

A  secondary  consideration  would  be  the 
availability  of  elements  necessary  to 
analytically  represent  the  device  being 
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modeled.  For  instance,  modeling  a  thin  membrane 
plate  structure  with  plate  elements  capable  of 
resisting  bending  forces  would  produce 
erroneously  large  stresses  and  lead  to  an 
overdesign.  However,  this  is  not  considered  a 
problem,  since  the  recommended  programs  have 
sufficiently  adequate  element  libraries  to 
describe  the  myriad  of  mi corel ectron ic  devices. 

Each  program  has  its  own  i d i osy nc r as i es  . 
This  was  illustrated  in  Sample  Problems  2 
through  4.  Sample  Problem  #2  showed  ANSYS  and 
ABAQUS  to  have  less  efficient  methods  for 
calculating  natural  frequencies  than  STARDYNE. 
Sample  Problem  #3  pointed  out  how  differently 
ABAQUS  and  ANSYS  treat  convection.  Finally, 
Sample  Problem  #4  showed  how  differently  ABAQUS, 
ANSYS,  and  STARDYNE  all  solved  the  same 
nonlinear  problem.  The  analyst  must  t.ry  to  be 
aware  of  program  idiosyncrasies  before  choosing 
a  program. 

Perform  Simple  Hand  Calculations 

Hand  calculations  can  provide  the  analyst 
with  a  "feel"  for  the  problem,  as  well  as  a 
check  on  results  obtained  from  the  finite 
element  model.  Depending  on  the  problem 
complexity,  this  step  may  or  may  not  be 
difficult.  Regardless  of  problem  complexity, 
the  analyst  must  make  some  assessment  of  the 
problem  before  he  starts  modeling. 

This  process  was  illustrated  in  Section  5.0 
of  this  report  -  Sample  Problems.  The  first 
three  sample  problems  showed  how  finite  element 
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results  can  be  compared  to  hand  calculations. 
However,  the  fourth  sample  problem  had  no 
closed-form  solution  available.  The  analyst  was 
then  recuired  to  interpret  results  based  on  a 
physical  interpretation  of  the  problem. 


Recognize  Typical  Microelectronics  Problems 

What  are  the  typical  microelectronics 
problems?  Based  on  this  writer's  experience, 
the  following  list  is  provided:  (1)  hybrids, 

DIPS,  flatpacks,  or  discrete  components  mounted 
on  printed  circuit  boards,  (2)  microware 
integrated  circuits  and  stripline  modules,  and 
(3)  printed  circuit  boards  in  general  including 
connectors,  solder  joints,  copper  runs, 
component  leads,  and  acceleration-sensitive 


components  . 

A  list  of  common  mi c roel ec tron i c  failures 
was  presented  in  Section  3.0  of  this  report  - 
Technical  Assessment.  This  list  is  reproduced 
here  to  help  the  analyst  identify  more 
specifically  microelectronic  problem  areas  that 
need  to  be  considered, 
o  Lack  of  hermeticity 
o  Die  bond  fail ure 
o  Broken  external  lead 


External  lead  corrosion 
Wirebond  and  wire  failures 
Microwave  package  failures 
Leadless  chip  carrier  solder  joint  failure 
Air  wound  inductor  breakage 
Connector  wear  or  lead  breakage 


Knowing  from  past  experience  what  typically 
causes  real-world  failures  will  help  the  analyst 
to  concentrate  on  solving  meaningful  problems, 
which  will  help  to  identify  potential  failures 
and  increase  device  reliability. 

Once  the  analyst  has  completed  the  above  steps,  he 
is  ready  to  develop  the  finite  element  model. 

7 .3  Model  Development 

Since  the  finite  element  model  is  a  mathematical 
idealization  of  the  actual  structure,  the  modeling  process  is 
critical  to  the  analysis'  success.  Once  the  model  is  defined, 
it  will  form  the  basis  for  analyses  which  may  consider  a  dozen 
environmental  or  loading  conditions.  The  finite  element 
modeling  procedure  involves  the  following:  (1)  grid 

optimization,  (2)  proper  element  selection,  (3)  restraints,  (4) 
material  properties,  (5)  types  of  analyses,  (6)  recuired  input 
data,  and  (7)  submitting  the  job  to  the  computer.  These  seven 
steps  will  now  be  described  in  detail. 

a.  Grid  Optimization 

The  process  of  defining  a  grid  (nodes)  is 
that  part  of  the  finite  element  process  which 
can  be  described  as  an  "art".  The  general 
rule-of-thumb  is  to  use  as  many  nodes  as  you  can 
afford.  Generally  speaking,  the  more  nodes  a 
model  has,  the  more  accurate  it  becomes. 
However,  the  strength  of  the  finite  element 
method  is  that  the  model  will  converge  an 
acceptable  answer  quickly  when  on  an 
" appropri ate"  number  of  nodes  is  used. 
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There  are  two  physical  constraints  that 
require  nodes.  The  first  is  physical 
interfaces.  Examples  include  where  two  planes 
of  a  component  intersect,  where  a  hybrid  lead 
attaches  to  a  circuit  board,  or  where  a  ciruit 
board  attaches  to  an  electronic  chassis.  The 
second  constraint  is  boundary  conditions. 
Examples  include  the  edges  of  a  circuit  board, 
the  ends  of  a  component  lead,  or  the  attachment 
surface  of  an  MIC  module. 

The  real  "art"  of  finite  element  modeling  is 
the  spacing  of  nodes  when  their  presence  is  not 
required  by  physical  constraints.  In  static 
analysis,  highly  stressed  regions  will  have  a 
higher  concentration  of  nodes  than  elsewhere. 
Dynamic  analysis  requires  that  regions  where 
higher  order  vibration  modes  need  to  be 
determined  will  also  have  a  higher  concentration 
of  nodes  than  elsewhere.  The  key  to  this  step 
is  the  analyst's  preparation  -  in  particular, 
his  hand  calculations  and  experience  will 
describe  where  regions  of  high  stress  or  higher 
order  vibration  modes  are  expected. 

When  static  loading  is  symmetrical,  device 
symmetry  can  be  used  to,  in  effect,  decrease  the 
model  complexity  by  only  modeling  a  portion  of 
the  device  that  exhibits  symmetry.  For 
instance,  Sample  Problem  #1  used  a  quarter 
model,  as  did  Critical  Problem  #1  -  Models  1  and 
2  and  Critical  Problem  #2  -  Model  1.  The 

analyst  is  cautioned  that  symmetrical  models 
should  only  be  used  when  the  modeled  device 
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possesses  symmetry  and  the  loading  is  static  and 
symmetr i cal. 

When  the  static  loading  is  asymmetrical  or 
vibration  modes  are  needed,  the  entire  device 
should  be  modeled.  This  techni cue  is  especially 
important  in  dynamic  modeling  since  few,  if  any, 
vibration  modes  are  symmetric.  Sample  Problem 
#2  used  this  techt.i  que  for  dynamics  analysis, 
whereas  Critical  Problem  #2  -  Model  2  needed 

this  modeling  method,  due  to  asymmetrical  static 
loading. 

Another  aspect  of  grid  optimization  is  the 
effect  of  a  program  which  generates  several 
nodes  and/or  elements  at  one  time.  The  analyst 
must  be  constantly  aware  of  this  program 
feature,  since  it  can  save  him  a  lot  of  time.  A 
discussion  of  this  techni oue  is  presented  in 
Sampl e  P robl em  #1 . 

A  commonly  used  technique  is  that  of 
creating  a  "first-cut"  model.  This  model  is 
designed  to  be  simple  and  fast,  but  still 

provides  the  analyst  with  an  "idea"  of  what  to 

expect  with  a  more  complicated  model.  This 
technique  was  used  in  "Critical  Problem  #1  - 

Model  1.  The  purpose  of  this  model  was  to 
obtain  hybrid  lid  deflections,  and  good 

correlation  with  closed-form  results  was 
obtai ned . 

As  the  analysis  progresses,  the  analyst 
learns  about  device  behavior  and  then 

concentrates  on  areas  that  may  indicate 
potential  failures.  this  was  illustrated  in 
Critical  Problem  #2  -  Model  2,  when  the  seal  was 
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introduced  into  the  model  and  in  Critical 
Problem  #2  -  Model  3,  where  a  highly  stressed 
region  of  the  seal  was  analyzed  with  a  local 
model  . 

Proper  Element  Selection 

The  most  common  finite  elements  are  shown  in 
Figures  49  and  50.  Choosing  the  proper  elements 
requires  the  analyst  to  have  some  preconception 
as  to  the  expected  structural  behavior  of  the 
device  being  modeled.  For  instance,  beam/bar 
elements  are  typically  used  to  model  component 
leads,  plate  elements  are  used  to  model  circuit 
boards,  and  solid  elements  are  used  to  model 
component  bodi es  . 

Many  elements  are  able  to  resist  different 
types  of  loading.  The  most  general  beam  element 
can  resist  all  three  types  of  loading  shown  in 
Figure  49.  However,  if  only  axial  loads  were 
expected,  then  the  general  element  should  be 
limited  to  resisting  only  axial  bar  loads.  If 
in-plane  membrane  forces  are  expected,  then  the 
membrane  plate  element  should  be  used,  etc. 
Refer  to  Appendix  D  for  a  glossary  of  finite 
element  terms  . 

Optimal  element  usage  can  be  enhanced  by 
understanding  element  behavior  and  element 
limitations.  Standard  tests  [57]  have  been 
developed  and  limitations  of  elements  [57,90] 
have  been  published  in  the  literature.  Since 
each  element  is  a  stand-alone  structure,  the 
analyst  needs  to  understand  the  element's 
theoretical  makeup.  When  this  is  done,  optimal 
use  of  the  elements  will  be  achieved.  No 
element  is  theoretically  perfect  for  all  types 
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Figure  49.  Beam  and  Plate  Elements 
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FI  g  ure 


50.  Plate  and  Solid  Elements 


of  usage.  However,  if  the  theory  is  understood, 
then  an  element's  strengths  can  be  best  utilized 
in  model  development. 

It  is  common  to  mix  element  types  within  a 
finite  element  model.  For  example,  a  model  of  a 
hybrid  attached  to  a  circuit  board  would  include 
beam,  plate,  and  solid  elements.  All  three  of 
these  elements  resist  different  types  of  forces. 
The  analyst  must  exercise  care  when  connecting 
elements  of  different  types,  since  force- 
compatibility  does  not  necessarily  exist.  The 
analyst's  knowledge  of  what's  expected 
physically  at  these  interfaces  determines  that 
restraints  are  to  be  applied  to  assure  transfer 
of  load. 

For  example,  consider  a  beam  element 
attached  at  one  common  node  to  a  plate  element. 
Assume  that  the  beam's  axis  is  perpendicular  to 
the  plane  of  the  plate.  The  beam  element  can 
resist  forces  in  the  three  translational  and 
three  rotational  directions  (six  degrees  of 
freedom  (DOF)).  However,  the  plate  element  can 
resist  all  forces  except  rotation  about  an  axis 
perpendicul ar  to  its  plane  (5  DOF).  Assume, 
now,  that  a  twisting  load  about  the  beam's  axis 
is  applied  at  the  beam's  end  opposite  where  it 
is  attached  to  the  plate.  This  load  will  get 
transferred  to  the  beam-plate  common  node  where 
the  plate  has  no  stiffness  to  resist  this  load. 
The  result  will  be  a  theoretically  infinite 
rotation  about  the  beam's  axis  and  the  program 
will  produce  a  warning  message.  This  problem 
can  be  alleviated  by  restraining  the  rotation 
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about  the  beam's  axis  at  the  common  beam-plate 
node,  and  is  discussed  in  detail  in  Sample 
Problem  #1.  Restraints  will  be  further 
discussed  in  the  next  portion  of  this  section. 

Throughout  the  critical  problem  analysis 
cubes  and/or  plates  were  used  to  model  the 
hybrid  walls,  lids,  and  seals.  When  considering 
plate  bending,  either  element  type  may  be  used. 
Plate  elements  are  more  commonly  used  because 
they  have  typically  half  as  many  nodes  (four, 
compared  to  eight)  as  a  cube  element.  However, 
plate  elements  do  have  more  DOF  per  node  (five 
compared  to  three)  than  cube  elements.  There¬ 
fore,  the  total  DOF  per  element  would  be  20  for 
a  plate  and  24  for  a  cube.  The  general  rule- 
of-thumb  is  to  use  plate  elements.  Cubes  are  to 
be  used  when  they  are  needed  to  physically 
describe  the  device's  geometry.  For  example, 
plates  were  used  in  Critical  Problem  #1  -  Model 
1  and  Critical  Problem  #2  -  Model  2.  Cubes  were 
used  in  Critical  Problem  #1  -  Models  2  and  3  and 
Critical  Problem  #2  -  Model  1. 

The  critical  problem  analysis  also  used 
cubes  and  membrane  plates.  The  STARDYNE  program 
does  not  allow  direct  application  of  pressure 
loading  to  its  cube  elements,  and  cube  stresses 
are  only  available  at  the  element’s  C.G.  There¬ 
fore,  attaching  thin  (.0001  in)  membrane  plates 
to  these  cube  elements  allowed  the  analyst  to 
apply  pressure  loading  and  to  obtain  stresses  at 
the  cube ' s  faces  . 
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c.  Restraints 

Restraints  are  one  area  of  modeling  that  is 
critical  and  not  necessarily  strai ght-f orward . 
Basically,  a  restraint  is  a  boundary  condition 
that  is  an  imposed  zero  displacement  at  the 
node.  Since  each  node  has  six  degrees  of 
freedom  (three  translations  and  three 
rotations),  a  maximum  of  six  displacements  can 
be  restrained  at  any  node.  Once  again, 
knowledge  of  what  happens  physically  determines 
what  restraints  are  to  be  applied. 

For  example,  connectors  on  circuit  boards 
are  usually  assumed  clamped  (all  six  displace¬ 
ments  restrained),  bolted  electronic  chassis 
attachment  points  are  usually  assumed 
ball-and-socket  joints  (all  three  translations 
restrained,  all  three  rotations  unrestrained), 
and  circuit  board  edges  supported  by  rubber  are 
assumed  simply  supported  (two  translations 
restrained,  and  two  rotations  restrained).  Test 
data  can  be  very  beneficial  in  determining 
difficult- to -describe  restraints  . 

The  .ample  problem  and  critical  problem 
analyses  provides  two  additional  restraint 
examples.  Sample  Problem  #1  illustrates  the 
procedure  for  applying  restraints  along  an  axis 
of  symmetry  when  using  a  symmetrical  model. 
Critical  Problem  #2  -  Model  1  shows  how  only  one 
node  need  be  fully  restrained  for  f ree-expans  i  on 
thermal  stress  models. 
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d.  Material  Properties 

Typical  material  properties  recuired  for  a 
dynamic  analysis  include  Young's  modulus, 
Poisson's  ratio  and  density.  Static  analyses 
require  Young's  modulus,  Poisson's  ratio, 
coefficient  of  expansion  (for  thermal  stresses), 
and  density  (if  weight  loading  is  involved). 
Depending  on  the  program,  properties  may  be 
linear,  nonlinear,  isotropic,  anisotropic,  or 
orthotropi c  . 

When  modeling  a  circuit  board,  the  copper 
circuitry  stiffens  the  board  and  its  effect  must 
be  included  in  the  Young's  modulus  value  used  in 
the  model.  An  "equivalent"  modulus  is  calcula¬ 
ted  in  this  case,  which  now  represents  the 
composite  effect  of  the  copper  and  board 
material.  This  calculation  is  aptly  described 
in  Steinberg's  book,  "Vibration  Analysis  for 
Electronics  Equipment"  [91].  Sample  Problem  #2 
also  illustrates  how  the  total  suspended  weight 
of  a  circuit  Sard  is  used  in  computing  an 
appropriate  density.  This  problem  also  pointed 
out  the  program  idiosyncrasies  that  STARDYNE 
recuired  a  weight  density  input,  whereas  ABAQUS 
and  ANSYS  required  a  mass  density  input. 
Finally,  Critical  Problem  #2  -  Model  1  showed 
how  an  equivalent  CTE  can  be  calculated.  This 
shows  how  a  simple  model  can  be  used  to 
determine  properties  necessary  for  a  more 
detai 1 ed  model  . 
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Types  of  Analyses 

The  two  types  of  mi croel ectroni c  analyses 
are  static  and  dynamic.  The  purpose  of  static 
analysis  is  to  apply  a  load  and  obtain  deflec¬ 
tions  and  stresses.  Typical  loads  include  point 
loads,  distributed  loads,  applied  temperatures 
(thermal  stress  analysis),  applied  displace¬ 
ments,  and  accelerations.  Dynamic  analysis 
first  involves  obtaining  a  number  of  natural 
freouencies  (eigenvalues)  and  mode  shapes 
(eigenvectors).  Then,  the  response  to  a 
particular  dynamic  environment  is  computed. 
Dynamic  environments  commonly  found  in 
microelectronic  analysis  include  sinusoidal 
vibration,  random  vibration,  shock  pulse,  shock 
spectrum,  and  transient  waveforms.  Finally,  the 
dynamic  response  can  be  used  to  calculate 
dynamic  stresses.  These  dynamic  stresses  can, 
in  turn,  be  used  to  assess  a  device's  fatigue 
life. 

Referring  to  the  typical  causes  of 
microel ectroni c  failures,  the  following  list  is 
provided  to  show  which  types  of  analysis  can  be 
performed  to  pre.ent  these  types  of  failures, 
o  Lack  of  hermeticity  -  static,  thermal 

stress ,  dynamic 

o  Die  bond  failure  -  static,  thermal  stress, 

dynami c 

o  Broken  external  lead  -  static,  dynamic 
o  External  lead  corrosion  -  usually  prevented 
by  proper  plating 

o  Wirebond  and  wire  failures  -  static,  dynamic 
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o  Microwave  package  failures  -  static,  thermal 
stress 

o  Leadless  chip  carrier  solder  joint  failure  - 
nonlinear  statics,  thermal  stress 
o  Air  wound  inductor  breakage  -  dynamic 
o  Conductor  wear  or  lead  breakage  -  thermal 
stress,  dynamic 

Recuired  Input  Data 

The  required  input  information  is  dependent 
on  the  type  of  analysis  and  the  program  being 
used.  Each  program  has  a  specific  format  for 
inputting  data.  The  user's  manual  provides  this 
format  which  must  be  followed.  The  input  data 
is  usually  recorded  on  computer  coding  sheets 
which  are  used  to  create  a  deck  of  punched 
cards.  Examples  of  this  are  shown  in  the  four 
sample  problems.  However,  recent  advances  i*n 
graphics  terminal  capabilities  have  improved 
these  procedures  . 

There  are  an  excellent  number  of 
preprocessors  available  which  allow  the  analyst 
to  create  and  check  a  model  graphically. 
Examples  of  the  procedures  are  documented  in  the 
literature  and  substantial  savings  can  be 
realized  by  using  a  preprocessor  [57],  The  basic 
strength  of  a  preprocessor  is  that  the  analyst 
can  see  his  model  displayed,  correct  any  errors  , 
and  submit  it  to  the  computer  error-free.  This 
is  shown  in  the  critical  problems  analyses  where 
the  Unistruc  preprocessor  was  used  to  generate 
isometric  views  of  the  finite  element  models. 
The  advent  of  color  graphics  [57]  and  a 
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seemingly  endless  number  of  improvements  makes 
this  topic  one  of  great  promise  and  should  form 
the  basis  for  future  work.  Future  work  would 
include  a  study  comparing  various  pre-post 
processors,  their  capabilities,  speed,  cost,  and 
acc  uracy  . 

g.  Submitting  the  Job  to  the  Computer 

If  a  preprocessor  were  used,  then  the  model 
will  be  debugged.  All  that's  required  is  to 

then  submit  the  job  to  the  computer,  calling  on 
the  appropriate  finite  element  program  to 
perform  the  analysis.  If  a  preprocessor  were 
not  used,  then  a  number  of  geometry  runs  will  be 
needed  until  the  model  is  debugged.  Then,  the 
analysis  can  be  accomplished  using  the 

appropriate  program. 

7  .4  Interpreting  Results 

For  a  static  analysis,  expected  results  include 
deflection  shapes  and  stresses.  Dynamic  analysis  results 

include  natural  frequencies,  mode  shapes,  plots  of  dynamic 

response,  and  dynamic  stresses.  A  graphics  postprocessor  is  an 
invaluable  tool  for  interpreting  deflection  shapes,  stress 
contours,  mode  shapes,  and  dynamic  responses.  Since  much  of 
finite  element  analysis  results  are  graphical,  displayed  results 
are  nearly  a  necessity. 

The  analyst  must  ask  the  question,  "Are  the  results 
reasonable?"  Simply  stated,  do  the  results  reflect  what's 
expected  and  reasonable?  Are  the  deflections,  stresses,  natural 
frequencies,  or  mode  shapes  consistent  with  those  found  by  hand 
calculations  or  experience?  Is  there  continuity  of  deflections 
and  stresses?  Do  the  applied  loads  e  qua!  the  re  actions?  All  of 
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the  above  cuestions  must  be  considered  by  the  analyst  once 
results  are  obtained. 

The  following  is  a  synopsis  of  how  the  results  of  the 
sample  problem  and  critical  problem  analyses  were  interpreted. 

o  Sample  Problem  *1  -  The  maximum  stresses  and 

deflections  compared  very  well  with  the  hand 
calculations,  considering  the  coarseness  of  the 
model  '  s  grid. 

o  Sample  Problem  #2  -  The  natural  frequencies 

compared  well  with  the  hand  calculations.  The 
mode  shapes  also  compared  well,  but  care  was 
needed  when  interpreting  (visually)  the  mode 
shapes  . 

o  Sample  Problem  #3  -  The  ANSYS  results  agreed 

closely  with  the  closed-form  solution.  The 
ANSYS  and  ABAQUS  programs  had  different  ways  of 
handling  convection. 

o  Sample  Problem  #4  -  No  closed-form  solutions 

were  available.  Therefore,  STARDYNE's  solution 
method  and  corresponding  answers  appeared  most 
valid  when  considering  a  physical  assessment  of 
the  problem. 

o  Critical  Problem  #1  -  Model  1  -  Good  deflection 
correlation  was  obtained  with  a  simple  "first- 
cut"  model.  Poor  stress  correlation  resulted 
because  of  a  coarse  grid  occurring  in  the  region 
of  maximum  stress  . 

o  Critical  Problem  #1  -  Model  2  -  Different  seal 
materials  resulted  in  varying  seal  stresses.  A 
more  detailed  model  of  the  seal  is  needed  to 
assess  the  vertical  stress  distribution  across 
the  seal's  width. 


182 


o  Critical  Problem  #1  -  Model  3  -  The  s  r  a  1 

vertical  stress  distribution  appears  reasonable 
when  considering  the  physical  aspects  of  the 
problem.  The  resulting  moment  appears  to  be 
reacted  by  the  inner  half  of  the  hybrid  wall.  A 
more  complicated  model  is  needed  to  further 
assess  seal  stresses  . 

o  Critical  Problem  #2  -  Model  1  -  The  quarter 

model  provided  stresses  for  establishing  static 
loading  for  hybrid  testing.  Nickel  plating  has 
little  effect  on  lid  CTE . 

o  Critical  Problem  #2  -  Model  2  -  Good  correlation 
was  obtained  between  analysis  and  test  lid 
deflections  under  static  loading.  Test  and 
analysis  CTE  at  room  temperatures  did  not 
compare  well,  but  thermal  strains  did  compare 
reasonably  over  a  range  of  temperatures.  Closer 
scrutiny  of  published  manufacturer's  CTE  values 
is  needed.  Lid  thickness  tolerance  effects  are 
sizeable  when  considering  lid  deflections. 
Worst  case  (thinnest)  lid  dimensions  should  be 
used  in  design  analysis. 

8.0  CONCLUSIONS  AND  RECOMMENDATIONS 

8 .1  S  ummary 

The  conclusions  and  recommendations  will  be  presented 
according  to  each  section  in  this  report.  The  conclusions 
represent  a  summary  of  the  information  gained  through  the  technical 
assessment  and  the  various  analyses  performed  on  microelectronic 
packages.  The  conclusions  presented  here  are  abbreviated  versions 
of  those  presented  at  the  end  of  each  section.  Refer  to  each 
section  for  a  complete  wording  of  appropriate  conclusions. 
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The  recommendations  will  be  directed  toward  the  need 
for  additional  work.  This  work  must  be  accomplished  with  respect 
to  information  gained  through  the  technical  assessment,  more 
detailed  analysis,  analysis  substitution  for  testing,  and  an 
indepth  study  of  pre  and  postprocessor  programs  in  conjunction  with 
CAD/CAM  acti vi ties  . 

8 .2  Cone! us i ons 

8.2.1  Technical  Assessment  (Phase  1) 

o  Although  FEA  is  primarily  a  tool  for  designers 
of  large  complex  structures  a  great  potential 
exists  if  applied  to  mi c roel ectroni c  packages, 
o  By  using  FEA  for  early  warnings  of  potential 

failures  the  mechanical  engineer  can  play  a 
significant  role  in  assuring  or  improving  the 
reliability  of  microelectronic  systems, 
o  The  survey  showed  that  the  areas  giving 

microelectronic  packages  the  most  problems  were 
hermeticity,  broken  or  corroded  internal  wires, 
and  broken  wirebonds. 

o  FEA  can  provide  accurate  answers  if  proper 

assumptions  are  made  and  if  one  understands  the 
finite  element  theory. 

o  Commercial  hybrid  houses  claim  lower  cost  and 

higher  reliability  than  similar  products 
produced  at  military  hybrid  houses, 
o  Most  companies  do  not  use  FEA  to  assess  the 

structural  and  thermal  integrity  of 
microelectronic  packages. 

o  The  "full  blown"  screening  in  MIL-STD-883  is  not 
recommended  by  most  of  the  companies  polled. 
Instead,  an  abbreviated  version  that  works  for 
them  is  utilized. 
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o  Closed-form  solutions  can  result  in  accurate, 

cost-effective  solutions. 

o  Considerable  cost  savings  can  be  realized  by 

substituting  FEA  for  MIL-STD-883  screening 
tests.  An  example,  showed  a  cost  savings  from 
$25,000  to  $50,200  for  a  lot  of  1,000  hybrids 
with  11  hybrid  types. 

8.2.2  Analysis  Methods  (Phase  2) 


o  Five  programs-ADINA,  FINITE,  MARC-CDC,  NASTRAN, 
and  SAP  IV  were  rated  “very  good"  and  have 
specific  solution  methods  that  were  rated 
“excellent".  However,  these  programs  should 
only  be  used  for  analyzing  unioue,  specific 
microelectronic  problems  that  address  the 
program's  strengths. 

o  The  four  programs  recommended  for* 
m  i  c roel  ec t r on i c  analysis  are  ABAQUS,  ANSYS, 
NISA,  and  STARDYNE  .  These  programs  will  best  be 
utilized  when  analyzing  the  following  types  of 
microelectronic  problems: 

a.  ABAQUS  -  nonlinear  statics  and  dynamics 

b.  ANSYS  -  heat  transfer  and  thermal  stress 

c.  NISA  -  composite  and  sandwich  structures 

d.  STARDYNE  -  linear  statics  and  dynamics 

o  STARDYNE  is  the  best  program  for  analyzing  the 
majority  of  microelectronic  problems. 


185 


8  .2  .3  Sample  Problems  Analysis  (Phase  2) 

8. 2. 3.1  linear  Statics  Problem 

o  ABAQUS,  ANSYS,  and  STARDYNE  had  good  correlation 
with  theoretically  determined  maximum  lid 
deflections  (within  3%)  and  stresses  (within  8%) 
with  a  relatively  coarse  grid. 

o  ANSYS  and  STARDYNE  deflection  shapes  compared 
well  with  each  other  -  ABAQUS  results  were 
slightly  different. 

8  .2  .3  .2  Linear  Dynamics  Problem 

o  STARDYNE  has  the  most  accurate,  user-oriented 
methods  for  obtaining  natural  frecuencies  and 
mode  shapes  . 

o  ANSYS  obtained  accurate  val ues  of  natural 

frequencies  and  mode  shapes--however ,  its  method 
is  not  efficient  and  is  therefore  more  expensive 
than  STARDYNE's  . 

o  ABAQUS  provided  inaccurate  values  of  natural 

frequencies  and  mode  shapes  -  its  method  is  not 
user-ori ented . 

8  .2 .3 .3  Heat  Transfer  Problem 

o  ANSYS  provided  an  accurate,  straight-forward 

solution  to  this  problem. 

o  ABAQUS  did  not  have  an  obvious  method  for 

handling  convection  for  this  simplified  problem. 

o  STARDYNE  does  not  have  heat  transfer  analysis 

capabi 1 i ty  . 
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8  .2  .3  .4 


Nonlinear  Statics  Problem 


8  .2  .4 


o  STARDYNE  provided  the  most  us e r- o r i en ted , 
acc urate  sol uti on  . 

o  ANSYS  provided  a  user-oriented,  yet  somewhat 
inaccurate  solution. 

o  ABAQUS  could  not  solve  this  problem  with  its 
four  node  plate  element. 


Critical  Problem  Analysis  and  Correlation  to  Test 


Data  ( Phase  3) 


8  .2. 4.1  Critical  Problem  #1 


o  Lid  deflections  obtained  with  a  simple  finite 

element  model  correlated  well  (within  3.7%)  with 
Libove's  closed-form  solution, 
o  Different  seal  materials  result  in  different 

seal  stresses  . 

o  The  vertical  stress  distribution  across  the  seal 
width  is  probably  not  linear  and  needs  a  more 
in-depth  model  to  better  determine  this 
di stri b  uti on  . 


8  .2  .4  .2  Critical  Problem  #2 

o  Lid  deflections  computed  with  a  simple  quarter 
model  fell  within  the  range  predicted  by 
Libove's  closed-form  solution, 
o  Good  correlation  (within  8%)  between  finite 
element  model  results  and  experimentally 
determined  hybrid  lid  deflections  under  point 
loading  was  achieved. 
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o  Design  analyses  should  use  minimum  hybrid  lid 

thickness  val ues  . 

o  Good  correlation  was  not  obtained  between 

predicted  and  measured  hybrid  CTE  at  room 
temperature-- however ,  reasonable  correlation  of 
thermal  strains  was  observed  over  a  wide 
temperature  range. 

8 .2 .5  Analysis  Guidelines  (Phase  4) 

The  guidelines  section  was,  in  itself,  a  set  of 
"conclusions  and  recommendations".  Therefore,  the  end  of  the 
guidelines  section  did  not  include  a  set  of  conclusions  and 
recommendations.  The  following  are  provided  as  general 
recommended  procedures  for  performing  mi c roel ec tron i c  finite 
el emen t  analysis  . 

1.  The  finite  element  process  involves 
understanding  the  problem,  developing  the  model, 
and  interpreting  results. 

2.  Understanding  the  problem  requires  the  analyst 
to: 

o  Know  what's  reouired 

o  Consider  program  choices 
o  Perform  simple  hand  calculations 
o  Recognize  typical  microelectronic  problems 

3.  Developing  the  model  recuires  the  analysis  to 
understand: 

o  Grid  optimization 

o  Proper  element  selection 
o  Restraints 

o  Material  properties 

o  Types  of  analysis 

o  Required  input  data 

o  Submitting  the  job  to  the  computer 
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4.  Interpreting  results  requires  the  analyst  to 
answer  the  following  Questions: 

o  Do  the  results  reflect  what's  reasonable  and 
expected? 

o  Are  the  deflections,  stresses,  natural 

frecuencies,  and  mode  shapes  consistent  with 
those  found  by  hand  calculations? 
o  Is  there  continuity  of  stresses  and 

deflections? 

o  Do  the  applied  loads  ecual  the  reactions? 

8  .3  Recommendati ons 

As  in  any  study,  many  interesting  problems  were 
discovered  that  were  beyond  the  scope  of  this  effort  but  should 
form  the  basis  for  future  studies.  Below  is  a  list  of 
recommendations  for  future  work: 

1.  Sealing  of  large  custom  packages  should  be 
investigated  to  gain  more  knowledge  on  fracture 
failures  of  brittle  materials  used  in  packages 
such  as  glasses  and  ceramics. 

2.  There  is  a  need  to  develop  a  criteria  for 
inspecting  mounted  hermetic  chip  carriers 
(leaded  and  leadless). 

3.  Military  specifications  should  include 
provisions  for  selective  cualification  by 
analysis.  To  fully  utilize  this  approach  more 
work  is  needed  in  the  tes t- anal ys i s  correlation 
area . 

4.  The  area  of  nonlinear  analysis  -  particularly 
solder  material  effects  -  is  a  promising  area 
for  future  work  and  should  be  investigated. 
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5.  More  work  is  reouired  in  the  area  of  hybrid  lid 
seal  analysis  to  more  accurately  determine  the 
stress  distribution  -  this  will  provide  more 
insight  to  design  criteria  for  sealing  large 
custom  packages  . 

6.  Pre  and  postprocessor  computer  program  studies 
should  be  performed,  similar  to  the  present  studies 
on  FEA,  to  determine  the  optimum  programs  for 
microelectronic  package  analysis.  Combined  with 
CAD/CAM  activities,  this  will  enable  designers 

to  improve  their  producti vi ty  . 
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APPENDIX  A 

HAND  CALCULATIONS  -  ROARK  -  CRITICAL  PROBLEM  #2 


From  Roark  [85],  page  392  ,  Case  8,  the  expression 
for  maximum  deflection  is: 

.4 

v  =  "Ob 

Jmax  £.3 


where : 


y 

■'max 

= 

maximum  lid  deflection 

a 

= 

.0247  for  a/b  =  1.66/1.06  =  1.566 

P 

= 

pressure  =  30  psi 

b 

= 

pi  ate  wi dth  =  1  .06 " 

E 

= 

Young's  modulus  =  20  x  106  psi 

t 

= 

lid  thickness  =  0.015" 

^max 

= 

0  .01386" 
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APPENDIX  B 

HAND  CALCULATIONS  -  LIBOVE  -  CRITICAL  PROBLEM  #2 


Referring  to  Libove  [88],  page  23,  ecuation  34,  the 
lower  bound  for  the  maximum  lid  deflection  is: 


^maxl 

= 

10.92  £  (f)3  an4 

where:  6  , 

maxi 

= 

maximum  lid  deflection 

P 

= 

pressure  =  30  psi 

E 

= 

Young's  modulus  *  20  x 

a 

= 

1  id  width  =  1 .06" 

t 

= 

lid  thickness  =  .015" 

n4 

= 

.00225 

6 _ 

- 

0  .01379" 
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APPENDIX  B  (CONTINUED) 


Referring  to  Libove  [85],  page  25,  ecuation  38,  the 
upper  bound  for  the  maximum  lid  deflection  is: 


^  max2 

= 

-IJU.l-v2-)..2  ,Pa2  „4  (0) 

EtJ 

where:  ^3x2 

= 

maximum  lid  deflection 

V 

= 

Poisson's  ratio  =  0.317 

P 

= 

Pressure  =  30  psi 

a 

= 

lid  width  =  1 .06 " 

E 

= 

Young's  modulus  =  20  x  10 

t 

= 

lid  thi ckness  =  .015” 

n4(0) 

= 

.008 

me 

= 

.446 

n9 

= 

.103 

^  max2 

= 

0  .04020" 

i 


1 


i 

i 
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APPENDIX  C 
TEST  PLAN 
FOR 

HYBRID  LID  TEST 
(CRITICAL  PROBLEM  #2) 

1.0  INTRODUCTION 

This  test  is  necessary  in  order  to  obtain 
correlation  between  the  finite  element  analysis  on  a  hybrid  lid 
and  test  data.  Two  tests  will  be  performed:  1)  A  concentrated 
load  with  deflection  measurements,  and  2)  a  high  and  low 
temperature  test  with  strain  gages  mounted  to  the  top  of  the  lid 
for  correlation  to  thermal  stress  analysis. 

2.0  TEST  PROCEDURE 

2  .1  Test  Sample 

The  test  sample  is  a  large  hybrid  (1.74“  x  1.14") 

shown  in  Figure  2.1-1.  The  lid,  shown  in  Figure  2.1-2,  is  a 
stepped  lid,  .015  inches  thick.  The  lid  is  welded  to  the  hybrid 
package  . 

2 .2  Test  Conditions  and  Measurements 

The  first  test  to  be  performed  shall  be  a 
concentrated  load  at  the  locations  shown  in  Figure  2.2-1  with 
deflection  measurements  being  made  at  the  point  of  load 

application.  The  base  of  the  hybrid  shall  be  supported  with  a 
solid  piece  of  aluminum  so  that  a  fixed  support  can  be  assumed 
in  the  analysis.  Loads  shall  be  applied  continuously  from  0.0 

to  3.0  pounds  at  the  defined  locations.  Deflections  shall  be 

recorded  for  the  maximum  load  as  shown  in  Table  2.2-1. 
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t  .003 


.  01  0 
MAX  TYP 


.035 

*«*-  DIA  TYP 


H  :i 

:!  i! 

'J  J  v 


u:  y 


.040  ,TF5 


i 


-  .002 


rr 


.018  +  .002-*JU- 

PIN  TYP 


1 


250  ♦  .01 
.042  DIA  MIN  TYP 


SECTION  B-B 

4  ground  pins  marked  "A"  to  be  flush,  or  underflush  with  inside 
surface  of  case. 


OUTLINE  AND  DIMENSIONS  (-401) 
Figure  2  .1-1 


Table  2.2-1.  Deflection  Measurements 


! 

t 

Deflection  measurement  for  points  A  -  E,  inches 
Load ,  1b.  ABODE 

3  .0 


Table  2.2-2. 


Strain  gage  location 

- sr  - 

S2 


Strain  Gage  Measurements* ,  PSI 

.  .  Temperature 

25°C  -55°C  2  5  °C  125° 


25° 


Measurements  taken  after  30  minutes  at  each  temperature  for 
stabil itation  . 
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The  second  test  will  be  a  high  and  low  temperature 
test  with  strain  gages  mounted  as  shown  in  Figure  2.2-2.  Strain 
gage  measurements  shall  be  taken  at  -55°C,  25°C  and  125°C  as 
shown  In  Table  2  .2-2  . 

3.0  RESULTS 

A  short  report  will  be  generated  to  document  this  test, 
the  results  and  equipment  used  for  correlation  to  the  finite  element 
analysis  . 
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APPENDIX  D 

GLOSSARY  OF  FINITE  ELEMENT  TERMS 


1.  Linear  Statics  -  statics  problems  that  can  be  solved  without 
an  iterative  process 

2.  Thermal  Stress  -  stresses  caused  by  temperature-only  loading 

3.  Nonlinear  Statics  -  statics  problems  involving  an  iterative 
sol uti on 

4.  Geometric  nonlinearity  -  nonlinearity  caused  by  gaps, 
bottoming  out  of  adjacent  structures,  one-way  springs, 
tension-only  and  compres s i on-onl y  members 

5.  Material  Nonlinearity  -  nonlinearity  caused  by  materials  which 
exhibit  nonlinear  behavior  such  as  elastic-plastic, 
viscoelastic,  creep,  or  swelling 

6-  Large  Deflection  Nonlinearity  -  nonlinearity  which  exceeds  the 
bounds  of  small  deflection  theory  engineering  mechanics 

7.  Linear  Dynamics  -  dynamics  not  involving  geometric,  material, 
or  large  deflection  nonlinearities 

8.  Modal  Extraction  -  the  process  of  obtaining  eigenvalues 
(natural  frequencies)  and  eigenvectors  (mode  shapes) 

9.  Transient  Response  -  dynamic  response  to  transients  such  as 
shock  pulses,  step  changes,  etc. 

10.  Harmonic  Response  -  dynamic  response  to  sinusoidal  vibration 
exci tati on 

11.  Random  Vibration  Response  -  dynamic  response  to  random 
vibration  excitation 

12.  Shock  Spectra  Response  -  dynamic  response  to  shock  spectrum 
exci tati on 

13.  Nonlinear  Dynamics  -  dynamics  involving  geometric,  material, 
or  large  deflection  nonlinearities 

14.  Linear  Heat  Transfer  -  heat  transfer  involving  only  linear 
phenomena  such  as  conduction 
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15.  Nonlinear  Heat  Transfer  -  heat  transfer  involving  nonlinear 
phenomena  such  as  convection  or  radiation 


16.  Bar  Element  -  element  capable  of  resisting  axial-only  forces 

17.  Beam  Element  -  element  capable  of  resisting  axial,  torsional, 
and  bending  forces 


18  . 

Membrane  Plate  Element  - 
forces 

el ement 

which 

can  resist 

i n- pi ane 

19  . 

Bending  Plate  Element  - 
forces 

el ement 

which 

can  resist 

bending 

20  . 

Thin  Plate/Shell  Element  - 
via  thin  plate/shell  theory 

el ement 
only 

that 

conforms  to 

loading 

21  . 

Thick  Plate/Shell  Element 

element  that  is  sensitive  to 

variations  in  stress  through  its  thickness 

22.  Isoparametric  Solid  Element  -  element  whose  displacement 
function  is  identical  to  its  shape  function  which  produces 
better  accuracy  within  the  element  and,  therefore,  requires 
less  elements  for  modeling 

23.  Axisymmetric  Element  -  an  element  used  to  model  structures 
which  have  circular  symmetry  about  an  axis  such  as  a  cylinder 

24.  Pipe  Element  -  a  beam  element  capable  of  resisting  fluid  flow 
loading 

25.  Gap  Element  -  one  dimens1'  nal  element  used  for  sensing  gaps  or 
bottoming  out  of  adjacen  structures 

26.  Friction  Element  -  one  dimensional  element  used  to  introduce 
f ri cti on 

27.  Spring,  Mass,  Damper  Elements  -  one  dimensional  dynamics 
elements  used  to  introduce  stiffness,  mass,  or  damping  at  a 
node 

28.  Node  -  a  point  within  a  model  which  defines  an  element's 
boundaries  or  a  point  where  element(s)  intersect 

29.  Node  Generation  -  the  process  of  automatically  generating  a 
series  of  e '"ally-spaced  nodes 

30.  Element  -  a  finite  continuum  defined  by  nodes 
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31.  Element  Generation  -  the  process  of  automatically  generating  a 
group  of  elements 

32.  Load  Generation  -  the  process  of  automatically  generating  a 
series  of  loads 

33.  Restart  Capability  -  the  ability  to  use  stored  model 
information  for  multiple  load  cases 

34.  Free  Format  Input  -  input  that  has  no  rigid  format,  such  as 
information  needing  to  be  only  input  in  certain  columns  of  a 
field 

35.  Data  Input  Check  -  a  capability  checking  for  correctness  of 
input  data  prior  to  making  a  complete  computer  run  with  loads, 
etc . 

36.  Substructuring  -  the  process  of  dividing  an  element  into  a 
detailed  model  within  itself  (a  substructure)  and  then 
assembling  the  substructures  into  a  large  model 

37.  Equivalent  Stiffness  Properties  -  replacing  the  composite 
stiffness  of  a  layered  section  with  an  equivalently  stiff 
uniform  section 

38.  Isotropic  -  having  properties  which  are  equal  in  all 
directions 

39.  Anisotropic  -  having  unequal  properties  in  all  directions 

40.  Orothotropic  -  having  material  properties  which  vary  in  three 
orthogonal  directions 

41.  Visoelastic  -  having  properties  that  cause  behavior  like  a 
thick  fluid  and/or  a  solid  material 

42.  Composite  -  a  layered  structure  that  has  predominately  planar 
properties  which  can  vary  according  to  direction 

43.  Sandwich  -  a  layered  structure  with  two  outer  bending  layers 
and  one  inner  core  which  transmits  shear  and  is  very  strong 
transversely 

44.  Deflection  Coupling  -  the  process  of  forcing  deflections  from 
two  adjacent  elements  to  be  equal  at  a  common  node 
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45.  Degree  of  Freedom  Releasing  -  the  process  of  not  restraining  a 
specific  degree  of  freedom 

46.  Closed  Form  Solution  -  an  analytical  approximation  (usually  in 
the  form  of  equations)  to  the  solution  of  a  physical  problem. 
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MISSION 

of 

Rome  Air  Development  Center 

RAPC  plant  and  executes  res  earch,  development,  tut  and 
t ejected  acquisition  programs  In  support  o f  Command,  Control 
Communications  and  Intelligence  (C3 I)  activities  .  Technical 
and  engineering  support  within  areas  of  technical  competence 
Is  provided  to  BSD  Program  Calces  IPOs]  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


